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Abstract 

Electronic  and  optoelectronic  devices  having  exceptional  perfonnance  capabilities 
are  being  fabricated  from  wide  band  gap  semiconductor  materials  such  as  gallium  nitride 
(GaN)  and  its  ternary  alloy  aluminum  gallium  nitride  (AlGaN).  These  devices  are  also 
attractive  candidates  for  integration  into  future  Air  Force  communication  and  sensor 
platforms,  including  those  that  must  operate  in  the  space  radiation  environment. 

However,  radiation-induced  perfonnance  degradation  in  GaN  and  AlGaN  devices  has 
been  observed,  and  this  has  been  attributed  to  the  creation  of  microscopic  point  defects 
such  as  vacancies,  interstitials,  and  related  defect  complexes,  which  have  been  associated 
with  energy  levels  deep  in  the  semiconductor  band  gap.  In  recent  years,  several  of  these 
defect  energy  levels  in  GaN  have  been  experimentally  observed  and  characterized,  but 
very  few  studies  have  been  performed  on  radiation-induced  defect  energy  levels  in 
AlGaN. 

In  this  dissertation  research,  the  electrical  and  optical  properties  of  high  energy 
electron-irradiated  n-AlxGai_xN  epitaxial  films  are  characterized  for  aluminum  mole 
fractions  from  x  =  0.0  to  0.3  using  deep  level  transient  spectroscopy  (DLTS), 
temperature-dependent  Hall,  and  cathodoluminescence  (CL)  measurements,  and  the 
results  are  compared  with  those  of  the  as-grown  n-AlxGai_xN.  The  DLTS  measurements 
reveal  the  presence  of  electron  traps  in  as-grown  GaN  that  are  characteristic  of  those 
designated  A,  B,  C,  D,  and  E  in  the  literature.  Also,  three  electron  traps,  labeled  D,  PI, 
and  P2  are  observed  in  the  as-grown  AlGaN.  Following  1.0  MeV  electron  irradiation  of 


IV 


the  AlGaN,  it  has  been  found  that  four  additional  electron  traps  labeled  Rl,  R2,  R3,  and 
R4  are  created  and  their  properties  are  characterized  for  the  first  time,  to  the  best  of  our 
knowledge.  Three  of  these  traps,  R1-R3,  correspond  to  radiation-induced  traps  reported 
in  GaN.  The  newly  discovered  fourth  trap,  R4,  appears  to  be  unique  to  AlGaN,  and  may 
thus  be  related  to  aluminum  displacement.  Furthermore,  this  new  electron  trap  R4  is  the 
most  prominent  radiation-induced  defect  in  the  DLTS  spectra.  Although  the  other 
radiation-induced  traps  begin  to  anneal  under  and  applied  voltage  bias  at  or  below  430  K 
in  AlGaN  Schottky  barrier  diodes,  this  R4  trap  is  thermally  stable  up  to  450  K,  and  could 
be  a  contributor  to  the  deleterious  radiation  effects  reported  for  AlGaN/GaN  high  electron 
mobility  transistors  (HEMTs).  It  has  also  been  found  that  trap  levels  deepen  significantly 
in  the  energy  band  gap  with  increase  in  aluminum  mole  fraction. 

The  room  temperature  concentration  of  electrons  in  the  conduction  band, 
determined  by  Hall  measurement,  is  found  to  decrease  significantly  in  both  GaN  and 
AlGaN  following  1.0  MeV  electron  irradiation  at  a  fluence  of  lxlO17  cm"2.  The  dose- 
averaged  carrier  removal  rate,  q,  was  found  to  depend  foremost  on  the  initial  carrier 
concentration,  n0,  regardless  of  the  aluminum  mole  fraction.  For  6.5xl016<  nQ<  8.2xl017 
cm"  ,  the  carrier  removal  rate  shows  a  linear  dependence  on  nQ  given  by  q  =  (3.96x10"  n0 
-  0. 15)  cm'1.  This  relationship  is  attributed,  in  part,  to  a  process  whereby  nitrogen 
interstitials  passivate  shallow  silicon  donors  by  forming  Si-N;  complexes. 

Following  1.0  MeV  electron  irradiation  at  a  fluence  of  lxlO17  cm"2,  the  peak  CL 
intensities  of  the  GaN  and  AlGaN  samples  were  reduced,  on  average,  by  50%.  This 
reduction  is  attributed  to  the  creation  of  non-radiative  recombination  centers. 
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For  all  of  the  GaN  and  AlGaN  samples,  the  observed  radiation  effects  of  carrier 
removal  and  luminescence  degradation  are  stable  at  room  temperature.  For  comparison 
of  radiation  tolerance,  a  sample  of  a  commonly  used  semiconductor  material,  gallium 
arsenide  (GaAs),  was  irradiated  and  characterized  along  with  the  GaN  and  AlGaN 
samples.  It  is  shown  that  GaAs  has  a  carrier  removal  rate  2  times  higher  than  the  GaN 
and  AlGaN,  and  the  luminescence  intensity  degrades  to  less  than  1%  of  the  original  value 
after  irradiation  with  the  1.0  MeV  electron  dose  of  lxlO17  cm"2. 
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DEEP  LEVEL  DEFECTS  IN  ELECTRON-IRRADIATED  ALUMINUM  GALLIUM 
NITRIDE  GROWN  BY  MOLECULAR  BEAM  EPITAXY 


I.  Introduction 

Throughout  the  last  decade,  GaN  and  related  Ill-nitride  alloys  have  emerged  as 
the  leading  materials  for  many  optoelectronic  devices  operating  in  the  blue  to  ultra-violet 
(UV)  spectrum.  More  recently,  GaN-based  electronic  devices,  such  as  AlGaN/GaN  high 
electron  mobility  transistors  (HEMTs),  have  received  much  attention  for  their  superior 
high-power,  high-frequency  perfonnance  characteristics  [1,2].  The  GaN  and  AlGaN 
device  technologies  are  still  immature  relative  to  those  of  GaAs  and  AlGaAs,  but  a  great 
deal  of  scientific  and  engineering  effort  has  been  focused  on  the  advancement  of  the 
GaN-based  device  technologies.  The  last  five  years  have  witnessed  a  steep  increase  in 
the  number  of  defect  studies  involving  thin  epitaxial  films  of  AlxGai_xN,  particularly 
epitaxial  GaN  (x  =  0).  These  studies  have  been  driven  by  several  different  technological 
challenges.  One  pressing  challenge  has  been  to  reduce  the  concentration  of  interface- 
related  defects  such  as  dislocations  and  stacking  faults  [3].  These  can  generally  be 
attributed  to  the  large  lattice  mismatch  and  difference  in  thermal  expansion  coefficient 
between  AlxGai.xN  and  available  substrates  such  as  ALO3  and  SiC.  Another  pressing 
challenge  is  to  achieve  efficient  p-type  doping.  For  example,  as-grown  native  defects 
create  problems  by  compensating  the  impurity  dopant.  Even  after  these  fundamental 
growth  problems  have  been  solved,  optimization  of  device  performance  requires 
characterization  of  the  remaining  defects  that  have  energy  levels  in  the  band  gap. 
Radiation-induced  defects  must  be  well  characterized  to  ensure  reliable  device  operation 
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in  a  space  or  nuclear  environment.  Processing  steps,  such  as  electrical  isolation  via  high- 
energy  ion  bombardment  [4,  5,  6]  and  dopant  activation  via  low-energy  electron 
irradiation  [7]  can  also  be  improved  when  the  relevant  defect  characteristics  are 
understood. 

Native  defects  have  been  studied  extensively  for  epitaxially  grown  GaN  (x  =  0), 
and  some  useful  conclusions  have  been  drawn.  For  example,  as-grown  undoped  GaN 
always  has  a  net  background  concentration  (typically  1016  -  10 19  cm'3,  depending  on  the 
growth  technique)  of  uncompensated  shallow  donors.  For  three  decades,  this  was 
attributed  to  the  nitrogen  vacancy,  Vn,  until  the  mid- 1 990 ’s  when  first-principles 
theoretical  calculations  indicated  that  the  formation  energy  for  a  Vn  is  too  high  to  occur 
so  readily  during  growth  [8,  9,  10].  To  date,  the  experimental  consensus  is  that  the 
dominant  residual  shallow  donors  are,  in  most  cases,  oxygen  or  silicon  impurities,  or  both 
[8,  9,  10].  Furthermore,  theoretical  studies  have  explained  why  native  donor  defects  are 
more  readily  formed  as  GaN  material  is  made  more  p-type,  and  why  native  acceptor 
defects  are  more  readily  formed  as  GaN  material  is  made  more  n-type.  That  is,  the  most 
likely  mechanisms  responsible  for  the  efficient  self-compensation  observed  in  GaN  are 
now  well  understood  [8].  Although  there  are  still  fundamental  questions  that  remain 
unanswered,  much  has  also  been  reported  recently  about  radiation-induced  defects  in 
GaN.  There  exists,  for  example,  some  consensus  about  the  primary  radiation-induced 
defect  energy  levels  in  GaN.  However,  no  such  consensus  exists  for  AlGaN  [11], 
because  only  a  couple  of  experimental  results  are  available  for  irradiated  AlGaN.  In  fact, 
a  clear  picture  has  still  not  emerged  regarding  the  dominant  as-grown  defects  in  AlGaN 
because  DLTS  measurements  have  not  been  reported  for  varying  aluminum  mole 


2 


fractions.  Specifically,  it  is  unclear  how  any  given  deep  defect  level  in  GaN  would 
change  with  aluminum  mole  fraction  in  AlGaN.  Consequently,  the  understanding  of 
performance-limiting  radiation  effects  in  AlGaN  electronic  devices  is  usually  based  on 
hypothesized  defect  energy  levels  that  have  yet  to  be  verified  [12].  Additionally,  there  is 
a  lack  of  dose-dependent  data  on  carrier  concentration  n  and  mobility  //  in  AlGaN.  Of 
particular  motivation  to  this  dissertation  research  is  the  fact  that  no  single  study  has 
correlated  measurements  of  radiation-induced  carrier  removal  with  aluminum  mole 
fraction  and  silicon  dopant  concentration.  These  and  other  such  gaps  in  knowledge  are 
compounded  by  the  different  growth  techniques  in  common  use.  GaN  and  AlGaN  layers 
grown  by  vapor  phase  epitaxy  (VPE)  techniques  such  as  metal-organic  vapor  phase 
epitaxy  (MO VPE)  and  hydride  vapor  phase  epitaxy  (HVPE)  have  received  the  most 
comprehensive  treatment  to  date.  Epitaxial  layers  of  AlxGai_xN  grown  by  MBE  generally 
contain  fewer  impurities  and  more  dislocation-related  defects.  These  different  initial 
conditions  for  the  thin  films  on  substrates  may  affect  the  radiation  response  of  the 
material,  and  thus  there  exists  the  need  for  a  comprehensive  study  of  radiation  effects  in 
MBE-grown  AlxGai.xN. 

In  this  dissertation  research,  the  electrical  and  optical  properties  of  high  energy 
electron-irradiated  n-AlxGai.xN  epitaxial  films  are  characterized  for  aluminum  mole 
fractions  from  x  =  0.0  to  0.3  using  deep  level  transient  spectroscopy  (DLTS), 
temperature-dependent  Hall,  and  cathodoluminescence  (CL)  measurements,  and  the 
results  are  compared  with  those  of  the  as-grown  n-AlxGai_xN.  The  main  purpose  of  the 
DLTS  measurements  was  to  obtain  accurate  DLTS  signatures  (total  activation  energies 
and  capture  cross-sections)  of  as-grown  and  radiation-induced  electron  traps.  Hall  effect 
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measurements  were  performed  in  order  to  characterize  room  temperature  carrier  removal 
due  to  1 .0  MeV  electron  radiation.  Additionally,  the  temperature  dependent  carrier 
concentration  was  fit  with  a  simple  theoretical  model  in  order  to  extract  radiation-induced 
defect  concentrations  and  activation  energies.  The  CL  characterization  was  perfonned 
before  irradiation  and  after  each  of  two  consecutive  1 .0  MeV  electron  radiation.  The 
scope  of  the  CL  analysis  was  limited  to  looking  for  significant  radiation-induced  spectral 
variation  and  quantifying  the  radiation-induced  degradation  of  the  dominant 
luminescence  peaks. 

This  dissertation  research  resulted  in  the  first  observation  and  characterization  of 
the  main  radiation-induced  electron  traps  in  AlGaN.  Also,  the  energy  levels  of  radiation- 
induced  and  as-grown  electron  traps  were  correlated  with  aluminum  mole  fraction  for  the 
first  time.  Furthermore,  the  degradation  of  carrier  concentration  and  luminescence 
intensity  in  AlGaN  was  correlated  with  aluminum  mole  fraction  and  initial  carrier 
concentration  for  the  first  time.  The  new  findings  enable  detailed  modeling  that  can  be 
used  to  engineer  radiation-hardened  device  designs.  The  findings  also  validate  some 
previous  speculations  and  invalidate  others  about  radiation  effects  in  AlGaN  compared  to 
GaN. 

The  following  chapters  begin  by  introducing  semiconductor  device  physics  and 
AlxGai_xN  properties  relevant  to  this  dissertation  research.  Then,  the  necessary 
background  of  defect  properties  and  radiation  effects  is  presented  with  review  of  related 
AlxGai_xN  literature.  The  experimental  and  analytical  procedures  used  to  obtain  the  data 
in  this  dissertation  research  are  presented  next.  Then,  the  experimental  results  are 


4 


presented  and  discussed.  Finally,  the  major  conclusions  resulting  from  this  dissertation 
research  and  recommendations  for  future  research  are  presented. 
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II.  Background 


2. 1  Physics  of  Semiconductor  Devices 


2.1.1  Semiconductor  Properties 


In  a  crystalline  material,  the  quantum  mechanical  interaction  of  conduction 
electrons  with  the  periodic  potential  of  the  host  atoms  gives  rise  to  allowed  energy  bands 
in  which  the  electrons  may  exist.  For  semiconductors,  this  phenomenon  gives  rise  to  the 
existence  of  an  energy  gap  between  the  band  of  electrons  tightly  bound  in  atomic  orbitals 
(valence  band)  and  the  energy  band  of  electrons  that  can  move  freely  throughout  the 
crystal  (conduction  band).  The  valence  band  is  normally  filled  and  the  conduction  band 
is  normally  empty  with  the  exception  of  a  small  percentage  of  electrons  that  are  thermally 
excited  from  the  valence  band  to  the  conduction  band.  This  concentration  is  given  as 


n  =  Nc  exp 


kT  / 


II- 1 


In  this  expression,  Ec  is  the  energy  at  the  conduction  band  minimum,  Ef  is  the  Fermi 
energy,  and  Nc  is  the  conduction  band  density  of  states  given  by 


where  Me  is  the  number  of  equivalent  conduction  band  minima,  and  m*  is  the  electron 
effective  mass.  Because  the  electrons  in  the  conduction  band  are  thermally  excited  from 
the  otherwise  full  valence  band,  the  intrinsic  concentration  of  holes  is  equal  to  the 
intrinsic  concentration  of  electrons  given  by  equation  II- 1 .  Holes  have  charge  +e  (e  = 
1.6022  x  10'19  Coulombs)  and  can  move  freely  throughout  the  valence  band.  The  Fenni 
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energy  of  an  intrinsic  semiconductor  is  located  near  the  middle  of  the  energy  band  gap, 
Eg.  Many  commonly  used  semiconductors  have  Eg  on  the  order  of  1  eV,  and  in  these 
semiconductors  the  intrinsic  carrier  concentration  at  room  temperature  ( kT ~  0.026  eV)  is 
many  orders  of  magnitude  less  than  Nc ■  For  many  applications,  the  intrinsic  carrier 
concentration  is  too  low  and  thus  extrinsically  raising  the  concentration  of  conduction 
band  electrons  or  valence  band  holes  via  doping  is  desirable.  In  fact,  extrinsic  doping 
takes  place  to  some  degree  in  all  semiconductors  due  to  unintentional  impurities  and 
defects  in  the  crystal  lattice.  A  donor  level  is  introduced  when  an  impurity  dopant  has 
one  more  valence  electron  than  the  atom  that  it  replaces  in  the  crystal  lattice,  and  an 
acceptor  level  is  introduced  when  the  impurity  dopant  has  one  less  valence  electron  than 
the  atom  that  it  replaces  in  the  crystal  lattice.  Defects  in  the  crystal  lattice  can  also 
introduce  centers  in  the  band  gap  that  behave  as  donors  or  acceptors.  A  donor  transition 
involving  one  of  these  centers  is  said  to  occur  when  a  neutral  center  loses  an  electron,  or 
when  a  positively  charged  center  gains  or  loses  an  electron  [13:109].  An  acceptor 
transition  is  said  to  occur  when  a  neutral  center  gains  an  electron  or  a  negatively  charged 
center  gains  or  loses  an  electron  [13:109].  Usually  the  tenn  “level”  is  used  to  indicate 
the  energy  of  a  particular  transition  [13:109].  The  occupation  of  a  donor  (acceptor)  level 
is  detennined  by  the  electron  distribution  function  as 


RED(A),T)  = 


1  +  —  expl 
Si 


JD(A) 

kT 


-E  ^ 


II-3 


where  ED(a)  is  the  donor  (acceptor)  energy  with  respect  to  the  conduction  band  minimum 
(valence  band  maximum),  and  the  constant  ( gjgi )  accounts  for  the  degeneracy  of  the 
level  with  (gi)  and  without  (g0)  the  trapped  electron.  The  appropriate  value  for  gQ/gi  is 
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often  1/2  for  donor  levels  and  4  for  acceptor  levels  [13:108].  With  a  donor  concentration. 
No,  the  concentration  of  donor  electrons  thermally  excited  to  the  conduction  band  is 

n  =  ND(l-f(ED,T)).  H-4 

In  the  presence  of  other  discrete  energy  levels  in  the  band  gap,  an  analogous 
expression  may  be  determined  from  the  charge  neutrality  condition  expressed  as  [13] 

«+ZA'i=/>  +  XA,i,-  n-5 

k  m 

This  expression  simply  states  that  at  thermodynamic  equilibrium,  the  sum  of  conduction 
band  electrons  and  ionized  acceptors  at  k  different  levels  must  be  equal  to  the  sum  of 
valence  band  holes  and  ionized  donors  at  m  different  levels.  In  extrinsic  semiconductors, 
either  n  or  p  in  equation  II-5  will  become  vanishingly  small,  with  the  remaining  carrier 
concentration  dictating  whether  the  semiconductor  is  n-type  or  p-type.  Considering  n- 
type  material,  for  which  the  hole  concentration  is  vanishingly  small,  equation  II-5  may  be 
written  as 

”  =  n-6 

m  k 

In  order  for  a  semiconductor  to  conduct  electricity,  it  must  posses  partially  filled 
electron  energy  bands  and  the  carriers  must  have  a  net  direction  of  motion  [14].  The  first 
condition  requires  thennal  ionization  processes  that  place  electrons  in  the  conduction 
band  or  holes  in  the  valence  band  as  discussed  above.  The  second  condition  requires  the 
action  of  an  electric  field  or  the  existence  of  a  carrier  concentration  gradient  in  the 
material.  In  the  presence  of  an  electric  field,  E  ,  the  coulombic  force,  F  =  qE  ,  causes 
holes  ( q  =  +e)  to  accelerate  in  the  direction  of  the  field  and  electrons  (q  =  -  e)  to 
accelerate  in  the  direction  opposite  the  field.  This  acceleration  is  interrupted  by  various 
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temperature-dependent  scattering  mechanisms  in  the  semiconductor  material.  There 
exists  a  mean  time,  r,  between  these  scattering  events  called  the  momentum  relaxation 
time  because  it  detennines  the  average  velocity  of  carriers  under  the  influence  of  a  given 
electric  field.  Using  the  notation  pertaining  to  electron  conduction  alone,  this  relationship 
between  the  average  electron  velocity  (called  drift  velocity)  and  electric  field  is  given  as 


m 


where  n„  ~  exjm  is  the  electron  mobility.  This  mobility  is  related  to  electron  current 
density  by  the  expression 

Jndrift  =  neMnE.  II-8 

The  total  drift  current  density  due  to  electron  and  hole  conduction  is  given  by  the  sum  of 
electron  and  hole  current  densities  as 

J  dnf,  =  e(n/un  +  p/up  )E  =  oE,  II-9 

where  a  =  e(np„  +  pup)  is  defined  as  the  conductivity  of  the  given  material.  The  same 
carrier  mobilities  are  also  important  to  diffusion  currents  that  result  from  carrier 
concentration  gradients.  Due  to  random  motion  of  carriers  alone,  there  will  be  a  net 
movement  of  carriers  from  regions  of  high  density  to  regions  of  low  density,  reflected  by 
a  net  current  density.  The  diffusion  current  due  to  electrons  is  expressed  as 

Jdijf  =  P,J<TVn  =  eDS/n,  II- 10 

where  Dn  =  pnkT/e  is  defined  as  the  diffusion  coefficient  of  the  given  material.  If  an 
electric  field  is  present  as  well  as  a  concentration  gradient,  the  current  densities  of 
electron  and  holes  are  given  by  [14]. 

Jn  =  -enjunE  +  eDVn,  II- 1 1 
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and 


JP  =  ePM„E  ~  eDnVn. 


11-12 


If  conduction  by  both  electrons  and  holes  is  present,  the  total  current  density  is  then 

J  =  Jn+Jp.  H-13 

The  concepts  outlined  above  lay  the  groundwork  for  understanding  electrical 
characterization  of  semiconductors  and  electronic  devices.  The  same  concepts  are  also 
relevant  to  the  characterization  of  optoelectronic  semiconductor  devices,  as  these  devices 
are  nearly  all  based  on  the  creation  or  annihilation  of  electron-hole  pairs  [14].  An 
electron-hole  pair  is  created  when  an  electron  is  excited  from  the  valence  band  to  the 
conduction  band.  As  mentioned  above,  this  excitation  can  take  place  spontaneously  with 
enough  thennal  energy,  or  it  can  take  place  upon  interaction  with  various  types  of 
radiation  that  have  energy  greater  than  the  band  gap.  In  either  case,  the  excited  electron 
will  eventually  recombine  with  the  hole  and  lose  its  excess  energy  in  the  process.  This 
energy  loss  can  be  radiative  or  non-radiative,  or  it  can  be  a  combination  of  the  two.  In  a 
radiative  recombination,  this  excess  energy  is  emitted  as  a  photon.  In  a  non-radiative 
recombination,  this  excess  energy  is  emitted  as  a  phonon.  Radiative  recombination  is  the 
basis  of  light  emitting  devices  and  luminescence  characterization  techniques.  Figure  II- 1 
shows  different  ways  in  which  radiative  recombination  transitions  frequently  take  place. 
In  figure  II- 1 ,  transition  (a)  is  a  band-to-band  transition  producing  a  photon  of  energy 

fico  =  E  g  +  kT  /  2,  11-14 

where  the  second  tenn  on  the  right  hand  side  of  the  equation  accounts  for  the  peak  in  the 
density  distribution  of  electrons  above  the  conduction  band  edge  [14].  Transition  (b)  is 
an  excitonic  transition,  of  which  there  are  two  types.  A  free  exciton  is  formed  when  a 
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Figure  II-l.  Various  radiative  transitions  observed  in  semiconductor  luminescence 
spectra. 


conduction  band  electron  and  valence  band  hole  are  bound  by  Coulombic  attraction.  The 
energy  of  the  electron-hole  system  is  reduced  by  a  small  free  exciton  binding  energy. 
When  the  electron  and  hole  recombine,  a  photon  is  produced  with  energy  [15] 


ha>  =  Eg 


*  4 

mre 
2 K2h2  ’ 


11-15 


where  k  =  m*r  is  the  reduced  effective  mass  of  the  exciton.  The  exciton  thus 

formed  can  also  form  a  bound  state  with  a  donor  or  acceptor,  and  this  is  referred  to  as  a 
bound  exciton.  When  a  bound  exciton  recombines,  the  emitted  photon  energy  is  reduced 
further  from  that  of  equation  11-15  by  an  additional  small  binding  energy  that  is  typically 
— 1/1 0th  the  donor  or  acceptor  energy  level  [14].  Transitions  (c)  and  (d)  are 
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donor/acceptor-band  transitions.  The  energy  of  the  photon  emitted  from  such  a  transition 
is  given  by  [15] 

ha>  =  Eg  ~ED{A)  +~kT.  11-16 

Here,  Ed(a)  is  the  ionization  energy  corresponding  of  the  donor  (acceptor)  level. 
Transition  (e)  is  a  donor  acceptor  pair  (DAP)  transition.  The  photon  energy  released  in  a 
DAP  transition  is  given  by  [15] 

2 

fi  co  =  E  —Ed  -  E  4  +  —,  11-17 

KT 

where  r  is  the  separation  distance  of  the  donor-acceptor  pair.  Phonons  can  be  emitted  in 
any  of  the  above  transitions,  and  whenever  that  is  the  case,  the  photon  energy  is  reduced 
by  the  energy  of  the  phonon(s)  involved.  In  the  case  of  sharp  peaks,  these  phonon 
replicas  can  be  clearly  identified  because  their  energies  are  shifted  down  from  the 
fundamental  peak  by  integer  multiples  of  a  known  phonon  energy. 

2.1.2  Metal-Semiconductor  Contacts 

In  order  to  fabricate  a  useful  electronic  or  optoelectronic  semiconductor  device,  it 
is  necessary  to  make  electrical  contact  with  the  semiconductor.  This  is  usually 
accomplished  by  depositing  metal  onto  a  clean  semiconductor  surface  in  a  vacuum. 
Figure  II-2(a)  shows  the  energy  band  diagram  of  an  isolated  metal  adjacent  to  an  isolated 
n-type  semiconductor  having  shallow  donor  concentration  No.  The  Fermi  level  Ef  of 
both  the  metal  and  semiconductor  are  labeled  as  well  as  the  energy  difference  between  Ef 
and  the  vacuum  level.  This  energy  difference,  called  the  work  function,  is  labeled  ecpm 
for  the  metal  and  ecps  for  the  semiconductor.  In  the  case  considered  here,  (pm  >  (ps.  The 
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Vacuum  level 


Interface 


Figure  II-2.  (a)  Energy  band  diagram  of  an  isolated  metal  adjacent  to  an  isolated  n-type 
semiconductor,  (b)  Energy  band  diagram  of  a  metal-semiconductor  contact  in  thermal 
equilibrium  (after  Sze  [16]). 
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energy  difference  ey„  between  the  conduction  band  minimum  and  the  vacuum  level  is 
called  the  electron  affinity.  When  the  metal  comes  into  contact  with  the  semiconductor 
as  depicted  in  figure  II-2(b),  the  Fermi  levels  of  the  two  materials  must  be  equal  at 
thennal  equilibrium.  The  vacuum  level  must  also  be  continuous.  These  two 
requirements  dictate  that  the  semiconductor  energy  bands  will  bend  such  that  E(  -  of  the 
semiconductor  is  eVbi  (=  e<ps  -  expm)  higher  at  the  interface  than  in  the  neutral  region  of  the 
semiconductor  [16],  The  quantity  Vhl  is  referred  to  as  the  built-in  voltage.  The  physical 
process  that  corresponds  to  this  energy  band  diagram  is  diffusion  of  donor  electrons  from 
the  semiconductor  to  the  metal.  This  diffusion  process  leaves  a  depletion  region  of 
ionized  donor  atoms  which  results  in  an  electric  field  directed  from  the  semiconductor 
toward  the  metal.  This  field  opposes  the  diffusion  current  of  electrons,  and  the  depletion 
region  increases  to  the  depth  W  required  to  exactly  counteract  the  diffusion  current. 

Thus,  as  required  by  definition,  no  net  current  flows  across  the  interface  at  thermal 
equilibrium.  Applying  the  depletion  approximation  that  ps  =  +eND  for  x  <  W  and  ps  =  0 
for  x  >  W,  solution  of  Poisson’s  equation  allows  W  to  be  expressed  as  [16] 

w  =  11-18 

\eND 

where  es  =  s,-e0  is  the  static  dielectric  constant  of  the  semiconductor.  If  the  potential  of  the 
neutral  semiconductor  region  can  be  controlled,  the  depletion  width  can  be  increased  or 
decreased  by  the  application  of  a  bias  voltage  V  across  the  interface.  With  the  applied 
bias  voltage,  the  expression  for  W  becomes 
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11-19 


The  applied  bias  voltage  V  in  equation  11-19  is  a  positive  quantity  when  the  voltage  on  the 
metal  is  made  positive  with  respect  to  the  n-type  semiconductor.  As  mentioned  above, 
the  electric  field  in  the  depletion  region  points  in  the  x-direction  (figure  II-2(b))  and  the 
field  strength  F  is  given  as  a  linear  function  of  x  by 

F(x)  =  eh^{W  -x).  11-20 

•A 

A  capacitance  is  associated  with  the  fixed  charge  in  the  depletion  region.  This  “space 
charge”  capacitance  per  unit  area  is  given  by  [16:164] 

Qsc  =  eN DW  =  p essND{Vb,  -V),  11-2 1 


and  the  capacitance  C  is  given  by 


C  =  A 


=  A 


e£sND  _  M 

2(V„,-V)  W 


11-22 


where  A  is  the  metal-semiconductor  contact  area.  The  metal-semiconductor  contact 
described  above  is  called  a  Schottky  contact. 

In  the  case  that  (ps  >  <pm  in  figure  II-2(a),  the  n-type  semiconductor  energy  bands 
bend  downward  upon  contact  with  the  metal,  and  there  is  an  accumulation  rather  than 
depletion  of  electrons  in  the  interface  region  of  the  semiconductor.  In  this  case,  the 
application  of  bias  V  causes  current  to  flow  through  the  contact  according  to  Ohm’s  law 
which  can  be  expressed  as 

V  =  IRC.  11-23 


This  expression  simply  states  that  the  current  /  is  directly  proportional  to  the  applied  bias 
voltage  Vby  the  contact  resistance  Rc  of  the  metal-semiconductor  contact.  A  contact 
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which  exhibits  this  behavior  is  called  an  Ohmic  contact.  For  semiconductor  devices,  it  is 
also  desirable  to  make  Ohmic  contacts  having  very  small  contact  resistances  such  that 
current  flows  nearly  unimpeded  through  the  metal-semiconductor  interface.  For  this 
reason,  an  Ohmic  contact  is  commonly  defined  as  one  for  which  Rc  ~  0.  Practically, 
metals  and  semiconductors  for  which  <ps  >  (pm  are  not  usually  available,  so  Ohmic  contacts 
must  be  fabricated  by  an  alternative  method  to  that  described  above  [14].  Alternatives 
include  making  the  semiconductor  interface  region  highly  n-type  by  implantation, 
diffusion,  or  alloying  [14].  Implantation  and  diffusion  are  both  ways  of  obtaining  a  large 
concentration  of  chosen  shallow  donors  in  the  interface  region  of  the  semiconductor. 
Alloying  is  usually  accomplished  by  depositing  a  properly  selected  material  on  the 
semiconductor  surface  and  annealing  at  temperatures  which  can  range  from  400  to  1000 
°C.  As  the  name  implies,  Ohmic  alloys  of  the  contact  material  and  the  semiconductor  are 
formed. 

Schottky  and  Ohmic  contacts  are  essential  to  semiconductor  device  fabrication.  A 
simple  and  important  device  called  the  Schottky  Barrier  Diode  (SBD)  consists  of  a 
Schottky  contact  as  described  above  with  an  Ohmic  contact  placed  elsewhere  on  the 
semiconductor.  The  Ohmic  contact  provides  control  of  the  semiconductor  potential  and 
thus  pennits  the  bias  V  to  be  applied  across  the  Schottky  contact.  A  diagram  of  a  SBD  is 
shown  in  figure  II-3.  In  figure  II-3,  the  Ohmic  contact  is  placed  on  the  back  surface  of  a 
bulk  semiconductor,  but  in  the  case  of  thin  epitaxial  semiconductor  layers  it  is  commonly 
placed  on  the  front  surface  next  to  the  Schottky  contact.  The  energy  band  diagram  of  the 
SBD  is  shown  in  figure  II-4  under  applied  bias  voltage  conditions  of  a)  no  bias  voltage; 
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semiconductor 


Figure  II-3.  Diagram  of  a  Schottky  Barrier  Diode  (SBD). 

b)  forward  bias  voltage;  and  c)  reverse  bias  voltage.  In  the  figure,  it  is  shown  that  a 
forward  bias  voltage  +Vf  causes  the  built-in  voltage  barrier  to  decrease  by  an  energy  eVf. 
Also,  it  is  shown  in  figure  II-4  that  a  reverse  bias  voltage  -  Vr  increases  the  built-in  voltage 
bias  by  an  energy  eVr.  Figure  II-4  also  shows  that  an  applied  bias  voltage  causes  the 
Fenni  levels  of  the  metal  and  semiconductor  to  misalign,  indicating  a  non-equilibrium 
condition  whereby  a  net  current  may  flow  across  the  metal-semiconductor  interface. 
Under  all  three  applied  bias  conditions,  it  is  seen  that  the  Schottky  barrier  energy  expn 
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faced  by  electrons  approaching  the  semiconductor  from  the  metal  side  remains  the  same. 
Decreasing  the  built-in  energy  barrier  eVbi  relative  to  expn  allows  more  conduction  band 
electrons  on  the  upper  end  of  the  Boltzmann  energy  distribution  to  overcome  the  barrier 
and  create  a  net  electron  current  across  the  interface.  The  transport  mechanism  thus 
described  is  called  thermionic  emission.  At  reverse  bias,  eVbi  is  increased  relative  to  etps, 
so  a  net  electron  current  flows  from  the  metal  to  the  semiconductor.  However,  etps  is 
relatively  large  in  a  good  SBD,  so  this  net  current  due  to  thermionic  emission  is 
approximately  zero  at  reverse  bias.  Thus,  the  SBD  is  a  rectifier  with  electrical  current 
density  given  by  [14] 


J  =  A**T 2  exp 


e<Ps 

kT 


exp 


\jifkT  j 


-1 


11-24 


where  A  is  the  Richardson  constant,  a  material  specific  quantity;  and  n/ -is  the  ideality 
factor  that  introduces  a  correction  to  account  for  other  transport  mechanisms  besides 
thennionic  emission.  Various  tunneling  mechanisms  and  carrier  generation  in  the 
depletion  region  are  examples  of  non-ideal  carrier  transport.  At  reverse  bias,  tunneling 
through  the  Schottky  barrier  is  an  important  current  transport  mechanism  that  causes 
undesirable  leakage  currents. 


2.2  Properties  of  AlxGa^xN  Material  and  Devices 

The  wurtzite  crystal  structure  of  GaN  is  shown  in  figure  II-5.  Though  a  zinc- 
blende  crystal  phase  of  GaN  can  be  grown  epitaxially  on  properly  chosen  substrates,  the 
wurtzite  phase  is  the  most  widely  grown  and  studied  [17].  Only  the  wurtzite  crystal 
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Figure  II-5.  Atomic  locations  for  wurtzite  crystal  structure  of  GaN.  The  graph  on  the 
left  illustrates  the  orientation  of  the  (0001)  direction,  with  Ga  represented  by  black  and  N 
by  white  balls.  The  graph  on  the  right  is  presented  looking  down  the  (000  1)  direction 
for  a  Ga  face  surface  [17]. 

phase  of  GaN  and  AlGaN  will  be  discussed  in  this  dissertation.  In  figure  II-5,  the  gallium 
atom  positions  are  represented  by  black  balls  and  the  nitrogen  atom  positions  are 
represented  by  white  balls.  The  axis  units  are  in  angstroms  and  the  lattice  constants  are  a 
=  3.189  A  and  c  =  5.185  A.  The  energy  band  structure  of  GaN  is  shown  in  figure  II-6.  It 
is  seen  in  this  figure  that  GaN  is  a  direct  band  gap  semiconductor.  That  is,  the 
conduction  band  minimum  and  valence  band  maximum  are  located  at  the  same  point  (T) 
in  momentum  space.  The  band  gap  energy  is  Eg  =  3.39  eV  at  T  =  300  K.  This  wide  band 
gap  affects  nearly  all  of  the  other  physical  properties  presented  here.  These  properties 
enable  the  desirable  perfonnance  capabilities  of  GaN-based  electronic  and  optoelectronic 
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Figure  II-6.  Wurtzite  GaN  band-structure  as  computed  by  Suzuki  [18],  et  al.  via  a  plane- 
wave  pseudopotential  calculation  method.  All  energy  values  are  at  300  K. 


devices.  AIN  has  the  same  wurtzite  crystal  structure  as  GaN,  but  aluminum  atoms  have 
lower  covalent  radii  than  gallium  atoms  [19:13].  The  crystal  lattice  constants  in  AIN  are 
smaller  than  those  in  GaN,  and  the  tighter  packing  of  atoms  contributes  to  larger  band 
gap  energy.  The  band  gap  energy  in  AIN  is  6.2  eV  at  T  =  300  K.  Table  II- 1  contains 
some  selected  physical  properties  of  GaN  and  AIN  that  are  directly  relevant  to  the 
analyses  and  discussions  in  this  dissertation. 

In  the  case  of  AlxGai.xN  where  aluminum  substitutes  for  gallium  on  x  %  of  the 
cation  lattice  sites,  the  physical  parameter  values  are  somewhere  in  between  those  of  GaN 
and  AIN.  For  many  physical  parameters,  Vegard's  law  is  a  very  good  way  to  determine 
the  parameter  value  in  AlxGai_xN  for  a  given  aluminum  mole  fraction  x.  Vegard's  law, 
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Table  II-l.  Selected  physical  properties  of  GaN  and  AIN  at  T  =  300  K  (from  [20]  and 
references  contained  therein). _ 


Property  (units) 

Symbol 

GaN 

Value 

AIN 

Value 

Energy  Gap  (eV) 

E* 

3.39 

6.2 

Density  (g/crn3) 

P 

6.15 

3.23 

Thermal  Conductivity 
(W/crnK) 

K 

1.3 

2.85 

Static  Dielectric 

Constant  (8.854  x  10"14 
F/cm) 

Sr 

8.9 

8.5 

Lattice  Constant  (A) 

a 

3.189 

3.11 

Lattice  Constant  (A) 

c 

5.185 

4.98 

Electron  Effective  Mass 
(m0) 

meff 

0.20 

0.48 

Electron  mobility,  bulk 
(cm2/V-sec) 

Pn 

1000 

135 

Hole  mobility,  bulk 
(cm2/V-sec) 

Mp 

30 

14 

Saturatio  n  Ve  loc  ity 
(cm/s) 

vsat 

2.5xl07 

1.4xl07 

Breakdo  wn  F  ie  Id 
(V/cm) 

Fh 

>  2.5xl06 

applied  to  GaN  and  AIN  may  be  expressed  as  [14] 

r(AlxGai.xN)  =  x  T(AIN)  +  ( 1  -x)  T(GaN)  ,  11-25 

where  T  represents  any  physical  parameter  for  which  the  relation  holds.  One  noteworthy 
physical  parameter  for  which  equation  11-25  may  not  be  adequate  is  band  gap  energy  Eg. 
When  calculating  Eg  of  a  ternary  alloy  it  is  often  necessary  to  introduce  a  quadratic  term 
to  equation  11-25  to  obtain  the  expression 

r(AlxGai.xN)  =  x  T(AIN)  +  (1-x)  T(GaN)  -  b  x(l-x),  11-26 

where  b  has  the  same  units  as  T  and  is  called  the  bowing  parameter.  In  the  case  of  Eg  in 
AlxGai.xN,  reported  bowing  parameters  have  ranged  from  b  =  0  [21]  (no  bowing)  to  b  = 
1.3  [22],  so  equation  11-26  can  not  presently  be  applied  with  much  certainty.  The 
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temperature  dependence  of  Eg  is  given  by  the  Varshni  equation  as  [17] 

Ji  2 

E(T)  =  E(  0) - .  11-27 

g  (T  +  p) 

In  GaN,  a  =  -7.7x1  O'4  eV/K,  =  600  K,  and  E(0)  =  3.49  eV.  In  AIN,  a  =  -1.8x1  O'3  eV/K, 
ft  =  1462  K,  and  E(0)  =  6.28  eV.  Figure  II- 7  shows  the  room  temperature  band  gap  of 
AlxGai_xN  calculated  with  equation  11-26  and  various  choices  of  b.  In  figure  II-8,  the 
temperature  dependence  of  the  b  =  0.35  curve  is  calculated  with  equation  11-27. 

Upon  application  of  equation  II- 1  and  II-2  it  is  found  that  the  intrinsic  carrier 
concentration  n,  of  AlxGai_xN  (0  <  x  <  1)  is  essentially  zero  at  T  =  300  K,  and  n,  remains 
small  enough  so  as  to  have  a  negligible  effect  on  the  operation  of  most  devices  until  T  ~ 
1000  K.  This  property  makes  these  and  other  wide  band  gap  semiconductors  well  suited 
to  use  in  high  temperature  environments.  The  wide  band  gap  of  AlxGai_xN  also  makes 
the  material  resistant  to  electrical  breakdown  in  high  electric  fields.  The  breakdown  field 
strength  of  GaN  is  approximately  10  times  higher  than  that  of  GaAs  and  Si.  GaN  also 
has  an  electron  saturation  velocity  more  than  2  times  greater  than  that  of  GaAs  and  Si. 
The  high  breakdown  field  and  electron  velocity  make  GaN  an  excellent  material  for  high 
power  and  high  frequency  electronic  devices.  As  are  all  Ill-nitride  semiconductors 
AlxGai_xN  is  highly  polar.  The  tensile  strain  caused  by  the  growth  of  AlGaN  on  GaN 
gives  rise  to  a  piezoelectric  polarization  that  combines  with  a  spontaneous  polarization  to 
create  a  net  positive  charge  at  the  AlGaN/GaN  interface  [2],  The  polarization  charge 
induces  a  2-dimensional  "gas"  of  electrons  (2DEG)  at  the  GaN  side  of  the  interface.  This 
2DEG  forms  the  basis  for  a  very  important  GaN-based  electronic  device  known  as  a  high 
electron  mobility  transistor  (HEMT).  A  schematic  cross  section  of  such  an  AlGaN/GaN 


23 


Al  mole  fraction  x 


Figure  II-7.  Calculated  room  temperature  band  gap  energy  of  AlxGai_xN  versus 
aluminum  mole  fraction  for  different  bowing  parameters. 


Al  mole  fraction  x 


Figure  II-8.  Calculated  band  gap  energy  of  AlxGai_xN  versus  aluminum  mole  fraction  at 
three  different  temperatures.  A  bowing  parameter  of  b=0.35  is  assumed. 
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HEMT  is  shown  in  figure  II-9.  In  figure  II-9,  a  dashed  line  indicates  the  location  of  the 
2DEG.  The  metal-AlGaN  contact  labeled  "gate"  is  a  Schottky  contact.  The  contact 
regions  labeled  "source"  and  "drain"  are  Ohmic  contacts.  When  a  bias  voltage  is  applied 
to  the  source  and  drain  contacts  current  flows  readily  in  the  2DEG  channel.  This  channel 
has  a  very  high  mobility,  especially  when  the  GaN  is  undoped.  Upon  application  of  a 
negative  voltage  to  the  gate  contact,  the  2DEG  diminishes  in  proportion  to  the  magnitude 
of  the  applied  gate  voltage  until  a  voltage  is  reached  at  which  current  is  no  longer  carried 
between  the  biased  source  and  drain.  The  speed  with  which  channel  conduction  can  be 
switched  on  and  off  and  the  high  source-drain  bias  voltages  that  can  be  supported  make 
AlGaN/GaN  HEMTs  well  suited  to  advance  the  state-of-the-art  in  high  frequency  power 
applications.  The  demonstrated  power  perfonnance  of  AlGaN/GaN  HEMTs  is  generally 
6-10  times  better  than  that  of  GaAs  or  InP  up  through  20  GHz  [23]. 
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Figure  II-9.  Schematic  of  an  AlGaN/GaN  HEMT. 


The  continuous  choice  of  direct  band  gap  energies  from  3.4  to  6.2  eV  make 
AlxGai_xN  a  very  attractive  material  for  optoelectronic  device  applications  in  the 
spectrum  ranging  from  blue  to  ultraviolet.  Functionally,  such  devices  can  be  categorized 
as  those  which  detect  light  and  those  which  emit  light.  The  metal-AlxGai.xN  SBD,  when 
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properly  designed,  functions  as  a  solar-blind  detector  for  x  ~  0.35.  This  Alo.35Gao.65N 
SBD  operates  at  reverse  bias.  When  a  photon  having  energy  greater  than  the  band  gap 
enters  the  depletion  region,  an  electron-hole  pair  is  created.  The  depletion  region  field 
causes  the  electron  and  hole  to  move  in  opposite  directions  out  of  the  depletion  region.  A 
current  thus  flows,  and  it  can  be  measured  in  external  circuitry.  The  tenn  "solar-blind" 
indicates  that  the  detector  is  blind  to  the  spectrum  of  solar  photons  found  at  the  surface  of 
the  earth.  These  solar  photons  have  wavelengths  X  longer  than  ~  285  mn,  and  they 
accordingly  can  not  excite  electrons  across  a  band  gap  greater  than  4.35  eV.  In  the  light 
emitter  category,  GaN-based  light  emitting  diode  (LED)  and  laser  diode  technology  have 
been  developed  to  the  point  of  significant  commercial  success  in  recent  years.  The 
physics  of  various  LED  and  laser  diodes  is  discussed  in  several  good  texts  [14,  16]  and 
will  not  be  presented  here. 

2.3  Defect  Characteristics 

There  are  many  types  of  native  point  defects  in  a  ternary  alloy  such  as  AlxGai_xN. 
Example  point  defects  include  vacancies  (VN,  Voa,  VAi),  interstitials  (Ni,  Gai,  Ah),  and 
antisites  (Noa,  GaN,  A1N).  Two  or  more  of  these  point  defects  can  also  couple  with  each 
other  (or  with  impurity  atoms)  to  form  electrically  active  defect  complexes.  Additionally, 
each  defect  type  can  have  more  than  one  possible  charge  state.  In  ternary  alloys  like 
AlGaN,  there  exist  additional  variations  due  to  different  cation  nearest-neighbor 
configurations  (e.g.  the  number  of  gallium  versus  aluminum  atoms  surrounding  a 
nitrogen  vacancy).  Which  of  these  defects  are  of  interest  depends  on  the  intended 
AlxGai_xN  device  application  and  on  the  relative  concentrations  of  the  various  defect 
species.  The  relative  concentration  of  a  given  point  defect  in  as-grown  semiconductor 
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material  can  depend  on  the  formation  energy  and  thennal  stability  of  the  defect  species 
[24],  the  stoichiometry  and  purity  of  feed  material  [8],  and  the  Fermi  level  [24].  In 
irradiated  semiconductor  material,  the  concentration  of  a  given  defect  can  depend  on  the 
type,  energy,  and  flux  of  radiation;  the  threshold  energy  for  atomic  displacements;  the 
thennal  stability  of  the  defect  species;  the  existence  of  defect  precursors;  and  the  Fenni 
level  [25,  26,  27]. 

An  electrically  active  defect  (that  is,  a  donor  or  acceptor)  is  normally  categorized 
as  being  either  shallow  or  deep.  This  tenninology  refers  to  the  magnitude  of  the  defect’s 
ionization  energy,  but  there  are  more  fundamental  characteristics  that  distinguish  the  two 
categories.  A  shallow  defect  level  tends  to  behave  as  a  “hydrogenic”  donor  or  acceptor 
with  ionization  energy  defined  by  the  canier  effective  mass  in  the  conduction  (valence) 
band.  A  carrier  occupying  one  of  these  states  tends  to  be  delocalized  in  space,  and  thus 
have  well  defined  momentum  according  to  the  Heisenberg  uncertainty  principle. 
Conversely,  the  large  ionization  energy  of  a  deep  level  suggests  a  more  tightly  bound 
electron  orbit  that  is  localized  in  space.  Accordingly,  deep  levels  tend  to  be  delocalized 
in  momentum  space.  This  results  in  two  significant  properties  of  deep  levels  [28],  First, 
deep  level  defects  are  efficient  recombination  centers  because  the  large  extent  in 
momentum  space  allows  the  level  to  couple  with  a  variety  of  phonons.  Nonradiative 
recombinations  involving  deep  states  are  then  much  more  likely  than  radiative 
recombinations,  with  phonon  emission  replacing  photon  emission  for  energy 
conservation.  Second,  the  entire  band  structure  must  be  involved  in  a  theoretical 
description  of  a  deep  defect,  making  its  energy  level  difficult  to  calculate  from  first 
principles.  The  nitrogen  vacancy  is  observed  in  both  as-grown  and  irradiated  GaN  layers. 
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First-principles  calculations  predict  that  Vn  is  a  donor  with  three  possible  charge  states 
(0,  1+,  and  3+)  in  AlxGai_xN  [29,  9,  30].  The  2+  state  is  unstable.  The  T2  state  of  the 
neutral  nitrogen  vacancy  is  resonant  with  the  conduction  band  in  GaN.  The  single  donor 
electron  in  that  state  automatically  ionizes  to  the  bottom  of  the  conduction  band  where  it 
forms  an  effective  mass  state  bound  by  the  Coulomb  tail  of  the  vacancy  potential  [24].  A 
central  cell  correction  potential  (the  short  range  component  due  to  local  chemical 
structure)  apparently  causes  this  level  to  be  40-50  meV  deeper  than  a  pure  effective  mass 
donor,  leading  to  the  measured  thermal  activation  energy  of  60-70  meV  [8,  31].  Thus, 

Vn  behaves  somewhat  like  a  shallow  single-electron  donor.  Under  the  effective  mass 
donor  model,  this  shallow  level  is  expected  to  become  deeper  with  increased  A1  mole 
fraction.  To  date,  these  energy  levels  have  not  been  systematically  measured  (i.e. 
temperature  dependent  Hall  (TDH)  and  deep  level  transient  spectroscopy  (DLTS)  for 
incremental  steps  in  x).  Around  x  =  0.45,  first  principles  calculations  have  indicated  that 
a  fundamental  change  in  the  electronic  structure  takes  place  [30],  in  that  the  Vn  neutral 
state  enters  the  band  gap  and  continues  deeper  into  the  band  gap  for  increasing  A1  mole 
fraction.  Zhou,  et  al.  [32]  observed  this  behavior  using  visible  absorption  spectroscopy 
of  proton-irradiated  AlxGai_xN.  For  x  <  0.45,  one  might  expect  the  Vn  donor  in  AlGaN 
to  be  similar  in  nature  to  Vn  in  GaN,  though  it  is  unclear  how  the  hypothesized  central 
cell  potential  that  is  so  significant  in  GaN  would  change  in  AlGaN.  The  triply  ionized 
state  is  expected  to  have  an  energy  level  about  300  meV  above  the  VBM  in  GaN.  This 
charge  state  is  metastable  in  n-type  GaN.  The  ubiquitous  electron  trap  in  GaN  that  is 
often  associated  with  VN  would  correspond  to  the  singly  ionized  state.  Look  et  al.  [33] 
reported  this  donor  level  to  be  64  ±10  meV  below  the  CBM,  and  almost  all  other  defect 
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studies  report  a  similar  level.  This  donor  energy  level  was  detennined  by  temperature 
dependent  Hall  effect  measurements,  and  it  has  been  observed  as  an  electron  trap  in  many 
different  DLTS  studies  as  well.  Most  of  these  DLTS  studies  reported  a  trap  level  around 
200  meV,  but  it  was  shown  that  this  is  probably  a  multi-component  manifestation  of  the 
same  Vn  level  having  different  capture  cross-section  and  capture  barrier  heights  [31]. 

One  plausible  explanation  for  these  different  components  is  variation  in  the  spacing  of 
Vn-Ni  Frenkel  pairs  [8].  Several  considerations  detailed  by  Look  [8]  make  it  appear  that 
Vn  (either  isolated  or  complexed  in  a  way  that  changes  the  energy  slightly)  is  the 
dominant  radiation-induced  donor-type  point  defect.  The  dominant  acceptor-type  defect 
induced  by  radiation  is  less  clear.  Theory  [24,  29]  points  unambiguously  to  VGa,  a  triple 
acceptor,  as  the  dominant  acceptor-type  defect  in  as-grown  n-type  material.  In  irradiated 
material,  Ni  point  defects  are  also  created  which  act  as  deep-level  single  acceptors.  Look 
[8]  makes  the  case  that  Ni  is  the  dominant  acceptor  in  irradiated  GaN.  The  as-grown  and 
radiation-induced  defect  levels  reported  in  GaN  and  AlGaN  to-date  are  summarized  in 
table  II-2  and  table  II-3.  The  format  of  the  tables  is  similar  to  that  presented  by  Look  [8]. 

The  annealing  behavior  of  the  above  point  defects  helps  correlate  different 
characterization  techniques  such  as  DLTS  and  cathodoluminescence.  The  typical 
annealing  trend  in  III-V  materials  is  for  interstitials  and  cation  vacancies  (Ni,  Gai,  and 
VGa)  to  anneal  out  at  temperatures  below  150  °C.  Anion  vacancies  (Vn)  tend  to  anneal 
out  in  the  temperature  range  of  200-400  °C,  and  antisites  (NGa  and  GaN)  tend  to  anneal 
out  at  temperatures  higher  than  600  °C  [32],  Look  et  al.  [33]  applied  a  first-order 
analysis  to  GaN  isochronal  (fixed  time  and  varied  temperature)  annealing  data  following 
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Table  II-2  Reported  defect  energy  levels  in  GaN. 


Category 

Experimental 
Energy 
as-grown  (eV) 

Introduced/ 
Increased  by 
Radiation 
(eV) 

Thermal 

Energy 

(eV) 

Type  of 
Radiation 

Observation 

Technique 

Material 

I.D. 

Ref. 

Donors 

0.06 

0.018 

0.06 

0.06 

e'  (1  MeV) 

TDH 

undoped  HVPE 
n~1017 

vN 

SlGa 

33 

Acceptors 

deep 

e'  (1  MeV) 

TDH 

II 

Ni 

33 

deep 

e'  (2  MeV) 

Positron 

Annihil. 

Semi-insulating  bulk 
[MgMOlMO20 

vGa 

44 

Electron 

Traps 

0.21,0.27, 

0.45,  0.61 

0.22 

0.22  peak  fit: 
0.06,  0.10, 

0.20 

e'  (90Sr) 

DLTS 

undoped  MOVPE  ELO 
n~1016 

vN 

34 

0.21,0.27, 

0.45,  0.61 

0.13,0.16, 

0.20 

H+  (2  MeV) 

DLTS 

undoped  MOVPE  ELO 
n~1016 

VN 

35 

0.21,0.27, 

0.45,0.61 

0.20,  0.78 

0.90 

He+ 

(5.4  MeV) 

DLTS 

undoped  MOVPE  ELO 
n~1016 

VN 

Spike 

36 

0.20 

0.12,  0.16, 
0.52,  0.59, 

0.90 

0.08 

H+  (24  GeV) 

DLTS 

undoped  HVPE 
n~1017 

vN? 

Ni  or 

Voa'Ni 

37 

0.27,  0.36,  0.58 

0.09,  0.13 

y-ray  (60Co) 

DLTS 

Si-doped  MOCVD 
n~1017 

vN 

38 

0.24,  0.45,  0.62 

0.18,  0.85 

0.06 

e'  (1  MeV) 

DLTS 

undoped  MOCVD 
n~1016 

vN 

39 

0.2,  0.55,  0.65, 
0.9 

DLTS 

undoped  freestanding 
HVPE  n~1016 

40 

Hole  Traps 

Ev  +  0.6,  0.9 

ODTLS 

undoped  freestanding 
HVPE  n~1016 

40 
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Table  II-3  Reported  defect  energy  levels  in  AlxGai_xN  (x  <  0.5) 


Category 

Experimental 
Energy 
as-grown  (eV) 

Introduced/ 
Increased  by 
Radiation 
(eV) 

Thermal 

Energy 

(eV) 

Type  of 
Radiation 

Observation 

Technique 

Material 

I.D. 

Ref. 

Electron 

Traps 

0.61 

x  =  0.12 
[Si]  ~  3  xlO17 
MOVPE  on  SiC 

similar  to 

0.5  eV  GaN 
level 

41 

0.77,  0.83,  1.01 

0.274 

0.187 

He+2 

(1.8  MeV) 

DLTS 

TSCAP 

x  =  0.12  and  0.42 
[Si]  ~  5  x  1017 

MBE  from  SVTA 

Similar  to 
0.234  eV 

GaN  level 

42 

0.567 

Similar  to 

0.61  eV 
level  in  [4 1  ] 
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0.7  MeV  electron  irradiation.  The  equation  that  describes  first-order  annealing  of 
mobility  is 

Mi1  =  M*  +  (/f  i  -  ) exp(-vt exp[-£ / kTt]),  11-28 

where  i  denotes  the  annealing  temperature  step,  visa  frequency  factor  (commonly  set 
equal  to  10  s'  ),  and  E  is  the  activation  energy  for  defect  migration.  Look  et  al.  used 
this  model  to  fit  the  data  precisely.  In  order  to  do  so,  they  increased  E  linearly  with 
temperature  in  the  range  250-400  °C.  This  would  be  expected  if  the  mobility  degradation 
is  caused  by  Frenkel-pairs  with  different  separations  [33].  Achieving  an  accurate  first- 
order  fit  to  the  data  in  this  way  suggests  that  vacancies  are  recombining  with  their  initial 
interstitials,  i.e.  Frenkel-pair  annealing.  Gallium  interstitials  and  related  complexes  were 
tentatively  identified  by  optically  detected  magnetic  resonance  (ODMR)  [43].  The 
doubly-charged  gallium  interstitial,  Gaf  ,  apparently  anneals  near  room  temperature  and 
complexes  to  form  a  different  paramagnetic  center  [43,44].  A  0.85  eV 
photoluminescence  (PL)  band  appeared  along  with  this  complex.  Gallium  vacancies  are 
detectable  by  positron  annihilation  spectroscopy.  Saarinen  et  al.  [44]  reported  that  the 
Voa  created  by  2  MeV  electron  irradiation  anneal  in  a  long-range  migration  process  at 
200-300  °C.  On  the  other  hand,  this  group  reported  that  the  VGa  in  as-grown  GaN  is 
stable  up  to  1000-1200  °C,  and  thus  posited  that  VGa  forms  complexes  with  oxygen 
impurities. 

Two  other  general  defect  categories  besides  point  defects  are  line  and  areal 
defects.  The  important  areal  defects  in  as-grown  GaN  are  stacking  faults,  and  these  are 
especially  concentrated  in  a  0.2  pm  thick  region  at  the  GaN/sapphire  interface.  This 
region  often  has  a  degenerate  donor  defect  concentration  which  may  be  related  to 
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stacking  faults  [45],  It  has  also  been  proposed,  however,  that  this  interface  conduction 
layer  is  caused  by  an  impurity  band  of  oxygen  [46].  Of  certain  consequence  to  the 
electrical  measurements  in  this  study  are  line  defects;  specifically  threading  edge  and 
screw  dislocations,  which  are  believed  to  be  electrically  active  in  GaN  [47].  These 
defects,  with  typical  densities  of  the  108-10l()  cm 2  in  epitaxial  GaN  grown  on  sapphire, 
are  a  consequence  of  the  14%  lattice  size  mismatch  and  34%  thermal  expansion 
coefficient  mismatch  between  the  GaN  and  the  sapphire  substrate  [47].  Some  theoretical 
calculations  suggest  that  the  core  of  such  defects  is  composed  of  Voa  [48].  Being 
acceptor-like  and  deep  in  the  band  gap,  these  ionized  defects  would  cause  a  significant 
reduction  in  mobility.  Indeed,  an  inverse  relationship  between  dislocation  concentration 
and  mobility  has  been  noted  by  many  observers  [47].  It  has  been  suggested  that  a  variety 
of  point  defect  species  segregate  around  core  dislocations  in  GaN  [49]. 


2.4  Radiation  Damage  Mechanisms 

2.4.1  Interactions 

To  first-order,  degradation  of  electrical  and  optical  characteristics  of  a 
semiconductor  usually  scale  with  non-ionizing  energy  loss  (NIEL).  NIEL  is  a  measure  of 
the  fraction  of  energy  that  is  transferred  from  the  bombarding  particle  to  primary  knock- 
on  atoms  (PKAs)  as  opposed  to  ionization  and  excitation  in  the  crystal  lattice.  This 
concept  is  very  useful  in  that  it  allows  a  direct  comparison  of  the  dose  imparted  from 
different  radiation  types  and  energies.  One  can  thus  analytically  detennine  the  fluence  of 
2  MeV  electrons  necessary  to  produce  the  same  displacement  damage  as  a  1012  cm'2 
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fluence  of  14  MeV  neutrons.  Additionally,  the  popular  Monte-Carlo  program  "Transport 
and  Range  of  Ions  in  Matter"  (TRIM)  [50]  may  be  used  to  calculate  the  actual  defect 
concentration  due  to  collision  chains.  Whereas  NIEL  is  an  analytically  detennined  value 
proportional  to  the  product  of  PKA  collision  cross-section  and  average  recoil  energy, 
TRIM  uses  the  PKA  recoil  energy  spectrum  to  incorporate  the  range  of  possible  collision 
chains  [25: 18],  The  TRIM  output  gives  the  concentration  of  displacements  generated  at  a 
given  material  depth,  but  this  says  nothing  about  the  thennal  stability  and 
electrical/optical  characteristics  of  the  resultant  defects. 

A  Frenkel-pair  is  generated  by  any  radiation  that  can  impart  the  necessary  displacement 
energy.  In  GaAs,  the  average  displacement  energy  is  about  9.8  eV  for  both  Ga  and  As 
atoms  [25].  The  precise  value  can  range  from  7-11  eV  depending  on  crystal  orientation 
[25:1 10].  For  a  relativistic  electron  to  impart  this  much  energy  in  an  elastic  collision,  it 
must  have  energy  greater  than  or  equal  to  the  threshold  energy,  Ethr.  The  maximum 
imparted  energy  is  given  in  units  of  MeV  by  the  following  equation  [8]: 

„  2E(E  +  2m0c2)  0.002 147£(£  + 1.022) 

^  max  —  ,  r  2  —  a  ?  H-29 

Me  A 

where  c  is  the  speed  of  light,  m0  is  the  rest  mass  of  an  electron,  M  is  the  mass  of  the 
target  atom,  A  is  the  atomic  weight  of  the  target  atom  in  units  of  amu,  and  E  is  the  energy 
of  the  bombarding  electron  in  MeV.  The  Emax  value  of  each  constituent  atom  in  AlGaN 
is  listed  in  table  II-4  for  the  given  electron  energies.  By  substituting  Ed  for  Emax,  equation 
11-29  may  be  solved  to  obtain  E,  which  in  this  case  is  equivalent  to  Ethr.  It  is  thereby 
calculated  that  the  threshold  electron  energy  for  Ga  and  As  atoms  in  GaAs  is  0.25  and 
0.27  MeV,  respectively.  Additional  transferred  energy  will  produce  Frenkel-pairs  with 
larger  separations  and  secondary  displacement  due  to  the  recoiling  target  atom. 
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Electrons  in  the  laboratory  testing  environment  and  in  the  near-earth  space  environment 
typically  have  energies  less  than  5  MeV  [27].  At  these  energies,  electrons  create  simple 
Frenkel  pairs  uniformly  in  semiconductor  material.  For  example,  1  MeV  electrons  have  a 
range  of  about  0.2  cm  in  GaN  [51].  The  1  MeV  electrons  traverse  a  1  pm  epitaxial  layer 
with 


Table  II-4.  Maximum  energy  transferred  to  AlGaN  atomic  constituents  for  given 
incident  electron  energies. _ 


Incident 
Energy  (MeV) 

Max  Energy 
Transferred  to  Ga 
Atom  (eV) 

Max  Energy 
Transferred  to  A1 
Atom  (eV) 

Max  Energy 
Transferred  to  N 
Atom  (eV) 

0.42 

18.7 

48.2 

92.9 

0.62 

31.4 

81.0 

156 

1.0 

62.3 

161 

310 

1.2 

82.1 

212 

409 

equal  displacement  probability  per  unit  depth  and,  following  a  displacement  collision,  the 
energy  of  the  electron  drops  below  Et/ir.  The  displaced  atom  becomes  an  interstitial 
relatively  close  to  the  remaining  vacancy. 

If  Ed  is  known,  the  Rutherford  cross  section  for  displacements  by  relativistic 
electrons  is  given  in  units  of  cm'  by  [33], 
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2  1/2 

where  Z  is  the  target  nucleus  charge,  e  is  the  electron  charge,  [3  =  [(>'"  -  1)  ]/y  ,  y  =  1  + 
E/me c2,  and  E  is  the  electron  beam  energy.  This  expression  includes  a  quantum 
mechanical  correction  formulated  by  McKinley  and  Feshbach  [52].  The  displacement 
cross  section  thus  calculated  is  used  to  predict  the  displacement  introduction  rate,  r\d,  of 
certain  atomic  displacements  during  electron  irradiation  by 
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Vd  =°dNti 


11-31 


where  Nt  is  the  target  concentration  of  the  particular  atomic  species.  The  displacement 
introduction  rate  has  units  of  cm’1  gives  the  concentration  of  the  displaced  atoms,  Nj, 
resulting  from  a  given  fluence,  ®,  by 

Nd=Vd  <*>•  H-32 

In  this  way,  an  introduction  rate  is  analogous  to  a  macroscopic  cross  section  in  nuclear 
engineering.  Figure  II- 10  shows  the  calculated  r)d  in  GaN  and  Al0.30Ga0.70N  for  a  range 
of  assumed  displacement  energies  (Ed)  based  on  the  displacement  cross  section  given  by 
equation  11-30.  The  displacement  energies  of  the  atomic  components  are  not  well  known, 
but  various  theoretical  considerations  and  experimental  results  [53,  44,  33,  54]  limit  the 
possible  range  of  values  to  10  <Ed<  50  meV  for  Ga  and  N  in  GaN  with  perhaps  the 
strongest  consensus  being  that  20  <Ed<  25  meV  [33,  54], 

Another  factor  that  can  conceivably  produce  different  damage  from  two 
irradiations  having  the  same  NIEL  is  beam  flux.  Increased  flux  causes  an  increase  in 
ionization  throughout  the  material  during  irradiation.  This  can  affect  the  mobility  of 
certain  defects  during  irradiation  via  the  increase  in  defect  charge  states  and  non-radiative 
recombination  energy  (phonons).  Additionally,  the  overlap  of  damage  trails  within  a 
short  time  period  can  cause  changes  in  the  stable  defect  concentrations.  In  the  case  of 
high-energy  electron  irradiation,  the  most  important  consideration  concerning  flux  is  the 
rate  at  which  heat  is  generated  in  the  sample.  If  the  sample  is  not  cooled  adequately  for  a 
given  flux,  the  sample  temperature  will  rise  and  increased  dynamic  annealing  will  take 
place. 
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Figure  11-10.  Calculated  displacement  introduction  rates  in  GaN  and  Al0.30Ga0.70N  for  a 
range  of  assumed  displacement  energies  (Ed)  based  on  the  displacement  cross  section 
given  by  equation  11-30. 
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2.4.2  Carrier  Removal 


Carrier  removal  is  a  general  term  referring  to  the  decrease  in  majority  carrier 
concentration  following  irradiation.  The  two  established  mechanisms  of  carrier  removal 
are  carrier  compensation  and  donor  passivation.  In  an  n-type  semiconductor,  the 
existence  of  any  acceptor-type  traps  below  the  Fenni  level  will  lower  the  free  carrier 
concentration.  The  acceptor-type  traps  that  are  more  than  a  few  kT  below  the  Fenni  level 
must  be  totally  (practically  speaking)  filled  in  thermal  equilibrium,  and  those  electrons 
are  supplied  by  the  donor-type  impurities  and/or  traps  above  and  near  the  Fermi  level. 
This  phenomenon  is  called  carrier  compensation,  and  it  occurs,  to  some  degree,  in  every 
real  semiconductor.  Donor- type  traps  that  are  well  below  the  Fenni  level  are  filled  and 
are  thus  unable  to  contribute  to  the  free  carrier  concentration.  Radiation  generally  creates 
both  donor  and  acceptor-type  traps,  but  the  net  effect  is  to  decrease  the  majority  carrier 
concentration,  as  described  above.  The  stability  of  certain  traps  depends  on  their  charge 
states  (i.e.  the  location  of  the  Fermi  level).  So,  the  relative  concentrations  of  stable 
donors  and  acceptors  also  vary  based  on  whether  the  material  is  n-type  or  p-type.  This 
additional  mechanism  is  apparently  responsible  for  type  conversion  in  n-type  GaAs  [55] 
(that  is,  the  n-type  material  not  only  becomes  more  intrinsic,  but  proceeds  to  become  p- 
type  under  irradiation). 

At  low  doses,  the  carrier  concentration  n  is  reduced  as  a  function  of  fluence  O 
according  to 

n(O)  =  n0  (l  -  aO)  ,  11-33 

where  the  constant  a  has  units  of  cm"  and  is  applicable  only  to  a  given  type  and  energy  of 
radiation.  The  constant  nQ  is  the  initial  carrier  concentration,  and  <I>  is  the  fluence  (units 
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2 

of  cm'  )  of  the  given  type  of  radiation.  The  carrier  removal  rate,  defined  as  //  =  dn/dO,  is 
equivalent  to  n0 a  in  equation  11-33  and  has  units  of  cm"1.  Though  a  simple 
proportionality  of  removal  rate  to  initial  free  carrier  concentration  is  implied  by  11-12,  the 
constant  a  can  in  principle  be  proportional  to  the  inverse  of  n0  or  a  fractional  power 
thereof.  For  energetic  ion  and  neutron  irradiation  of  GaAs,  a  is  inversely  proportional  to 
a  fractional  power  of  the  initial  carrier  concentration,  as  will  be  discussed  below. 

When  the  fluence  is  large  enough  to  cause  a  significant  decrease  in  majority 
carrier  concentration,  equation  11-33  begins  to  become  less  reliable.  The  main  reason  for 
the  deviation  from  the  linear  model  is  that  trap  filling  is  a  non-linear  function  of  the 
existing  concentration  of  majority  carriers,  which  is  decreasing  significantly  with  fluence 
[].  In  other  words,  at  high  doses  the  Fermi  level  begins  to  move  towards  the  middle  of 
the  band  gap  and  trap  filling  in  the  upper  half  of  the  band  gap  decreases  accordingly. 
Dose-dependent  carrier  removal  in  GaAs  has  been  studied  extensively  in  the  past  decade 
for  the  purpose  of  electrical  isolation  via  ion  implantation  [56,  57].  Prior  to  that,  a 
number  of  neutron  [58,  59,  60]  and  electron  [61,  55]  irradiation  studies  were  performed 
in  order  to  characterize  intrinsic  defects  in  GaAs  and  understand  radiation  effects  in 
GaAs  devices. 

For  GaAs,  the  linear  model  of  equation  11-33  is  only  applicable  up  to  doses  where 
carrier  concentration  decreases  by  25%  [58].  At  higher  doses,  the  decrease  is  seen  to  be 
approximately  exponential,  corresponding  to  a  linear  shift  in  the  Fenni  level  with  dose. 
That  is, 

n(O)  =  n0  exp(-aO)  ,  11-34 

where  the  constant  a  takes  on  a  value  that  is  different  than  in  equation  11-33. 
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An  additional  consideration  that  is  certain  to  apply  to  GaN  is  the  efficient 
dynamic  annealing  that  takes  place  during  room  temperature  irradiation.  That  is,  the 
introduction  rate  of  stable  and  observable  defects  is  may  be  much  less  than  the 
introduction  rate  of  primary  point  defects  (vacancies,  interstitials,  and  antisites).  This 
behavior  has  been  studied  in  connection  with  electrical  isolation  via  ion  implantation 
damage  in  GaN.  In  one  study  [62],  point  defects  were  found  to  anneal  efficiently  during 
irradiation  at  both  77  K  and  room  temperature.  Stable  defect  complexes  were  fonned  in 
the  process,  but  no  amorphization  was  observed  by  Rutherford  backscattering  (RBS)  and 
transmission  electron  microscopy  (TEM)  until  a  threshold  dose  was  reached.  At  that 
threshold  dose,  an  amorphization  layer  formed  on  the  surface  and  grew  toward  the  bulk. 
At  this  point,  it  is  important  to  clarify  that  the  doses  required  to  achieve  electrical 
isolation  are  several  orders  of  magnitude  lower  than  those  required  to  achieve 
amorphization.  In  the  fonner  case,  the  threshold  dose  is  that  required  to  produce  defect 
concentrations  on  the  order  of  the  number  of  initial  free  carriers  (~10  cm"  ).  In  the  latter 
case,  the  threshold  dose  is  that  required  to  produce  defect  concentrations  on  the  order  of 
the  crystal  atomic  concentration  (--TO”  cm'  ).  What  the  RBS  study  makes  clear  is  that 
some  defects  anneal  in  GaN  at  temperatures  as  low  as  77  K. 

Dose-dependent  studies  of  carrier  removal  are  not  widely  reported  for  either  GaN 
or  AlGaN.  Irradiation  with  1  MeV  electrons  was  reported  to  introduce  compensating 
acceptors  at  a  rate  of  approximately  1  cm"  in  GaN  having  n0  =  2.3x10  cm"  [31]. 

Green  [17]  reported  a  carrier  removal  rate  of  0.14  cm'1  for  1.0  MeV  electron  irradiation 
of  unintentionally  doped  freestanding  GaN  layers  grown  by  HVPE.  For  1  MeV  proton 
irradiation  of  unintentionally  doped  GaN  ( nQ  =  1.2x10  cm"  ),  Auret  et  al.  [63]  found  a 
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carrier  removal  rate  of  3880  ±  380  cm’1  using  C-V  measurements.  They  also  used  DLTS 
to  measure  the  introduction  rate  of  the  ER3  defect  as  being  290  ±  50  cm"1. 

Some  groups  have  measured  dose-dependent  sheet  resistance  in-situ.  Boudinov  el 
al.  [64]  reported  that  sheet  resistance  in  n-type  GaN  (n0  =  3x10  cm'  )  increases  by  50% 
when  exposed  to  a  0.6  MeV  proton  fluence  of  5xl013  cm’2.  The  threshold  dose  (the  dose 
at  which  resistivity  reaches  a  plateau  following  several  orders  of  magnitude  increase)  is 
5x10  cm’  .  It  was  also  reported  that  the  threshold  dose  is  directly  proportional  to  the 
initial  carrier  concentration  in  silicon  doped  GaN  [64].  This  is  in  qualitative  agreement 
with  the  linear  carrier  removal  model  of  equation  11-33.  More  recently,  TDH 
measurements  were  reported  for  the  same  type  of  GaN  sample  discussed  above  [62]. 
Following  a  single  dose  of  2.4xl014  cm’2,  the  sheet  carrier  density  dropped  by  about  5 
orders  of  magnitude  whereas  the  effective  Hall  mobility  dropped  by  about  1  order  of 
magnitude.  A  recently  proposed  model  [65]  with  a  simple  theoretical  basis  is  consistent 
with  the  above  observations  (specifically,  the  resistivity  vs.  dose  curve  and  the  threshold 
dose  dependence  on  n0 ).  The  model  is  based  on  the  transformation  of  shallow  dopant 
levels  into  deep  levels  via  complexes  that  are  formed  between  point  defects  and  the 
dopants. 

It  is,  at  present,  unclear  how  the  aluminum  mole  fraction  in  AlxGai_xN  affects 
dose-dependent  carrier  removal.  A  carrier  removal  rate  (applicable  at  low  doses)  has 
been  reported  for  Alo.12Gao.88N  [42]  and  a  threshold  dose  for  isolation  has  also  been 
reported  [66].  Legodi  et  al.  [42]  concluded  that  the  carrier  removal  rate  in  AlxGai_xN  (x  = 
0.12  or  0.41)  due  to  1.8  MeV  He“  ions  is  80  times  higher  than  in  similarly  irradiated 
GaN.  This  may  be  the  only  published  statement  on  the  subject,  thus  far,  and  the  result  is 
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not  at  all  conclusive.  The  GaN  carrier  removal  study  used  for  comparison  is  not  cited. 
One  of  the  few  published  reports  of  carrier  removal  rate  in  GaN  is  that  done  by  two  of 
Legodi’s  co-authors,  Auret  and  Goodman  [36]  for  5.4  MeV  He2+  irradiated  GaN.  The 
dose  deposited  through  a  gold  contact  in  both  studies  can  be  compared  using  the  TRIM 
program  [50],  Such  a  calculation  reveals  that  the  damage  created  per  unit  fluence  in 
Auret’s  study  is  1.35  times  greater  than  in  Legodi’s  study.  Thus,  the  irradiation 
conditions  of  these  two  studies  are  indeed  similar.  The  GaN  carrier  removal  rate  reported 
by  Auret  is  6200  ±  300  cm'1.  After  accounting  for  the  ion  energy  difference,  this  is  only 
5  times  lower  (not  80  times!)  than  the  carrier  removal  rate  reported  for  AlGaN  by  Legodi. 
Additionally,  the  GaN  sample  in  Auret’s  study  was  unintentionally  doped  with  n0  = 
3xl016  cm'3.  The  AlGaN  sample  in  Legodi’s  study  was  silicon-doped  with  n0  =  5-9  xlO17 
cm'  .  Thus,  it  is  unclear  from  these  studies  how  carrier  removal  would  compare  in  GaN 
and  AlGaN  with  similar  doping  levels. 

Structural  disorder  was  measured  in  proton  and  heavy  ion-irradiated  AlxGai_xN 
[62].  Increasing  the  A1  mole  fraction  (from  x  =  0.05  to  x  =  0.60)  was  shown  to  increase 
the  dynamic  annealing  at  room  temperature  [32],  Thus,  the  samples  with  the  highest  A1 
mole  fraction  were  the  most  resistant  to  amorphization.  In  contrast  with  GaN, 
preferential  surface  disordering  was  not  observed  in  AlxGai_xN  [32], 

2.4.3  Carrier  Trapping,  Generation,  and  Recombination 

The  occupation  of  the  defect-related  energy  levels  mentioned  above  was 
discussed  in  the  context  of  the  semiconductor  being  in  thermodynamic  equilibrium.  That 
is,  the  occupation  probability  of  each  level  was  defined  by  its  energetic  proximity  to  the 
Fenni  level.  When  a  non-equilibrium  situation  is  created  by  a  stimulus,  such  as  the 
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application  of  an  external  electric  field  or  exposure  to  ionizing  photons,  a  non¬ 
equilibrium  distribution  of  electrons  is  created.  When  the  stimulus  is  removed,  the 
carrier  distribution  at  each  level  will  return  to  thennodynamic  equilibrium  with  an 
emission  rate  defined  by  its  energetic  proximity  to  the  conduction  band  minimum  and  the 
temperature  of  the  material.  The  situation  described  above  is  referred  to  as  carrier 
trapping.  Specifically,  an  electron  trap  is  defined  as  a  defect  for  which  the  electron 
capture  rate  constant,  cn,  is  much  larger  than  the  hole  capture  rate  constant,  cp  [28],  A 
microscopic  capture  cross  section  <r,hp,  having  units  of  cm',  may  be  defined  as 


cn  =  crn<vn>n 

11-35 

and 

Cp  =  <7p<Vp>p 

11-36 

for  electron  and  hole  capture,  respectively.  The  symbol  <v,hp>  denotes  the  root  mean 
square  (rms)  value  of  the  carrier  thermal  velocity  given  as 


<  V„,P  >= 


3  kT 


m 
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where  m  is  the  effective  mass  of  the  applicable  carrier  type.  The  emission  rate  constant, 
en,  of  a  trapped  electron  to  the  conduction  band  has  the  form  [28] 


e„  =  A,  exP 


~E,h 

kT 
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where  Et/,  in  the  Boltzmann  factor  is  the  energy  difference  between  the  conduction  band 
minimum  and  the  trap  energy  level  (Eth  =  Ec  -  Et),  and  the  constant  An  is  a  property  of 
the  specific  defect.  Applied  to  an  electron  trap,  the  statistical  principle  of  detailed 
balance  states  that  in  a  neutral  material  at  thermal  equilibrium  the  capture  and  emission 
rates  must  be  equal.  This  principle  is  expressed  as  [72] 


enP  =  cn(\-P), 


11-39 
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where  P  is  defined  as  the  probability  that  the  trap  is  filled  by  an  electron.  This 
probability  is  given  by  the  Fenni-Dirac  distribution  function  of  equation  II-3, 


P  =  f(ET)  = 


1  +  —  exp 
gi 


Et  -  Ef 
kT 
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The  coefficient^,,  in  equation  11-38  may  be  detennined  by  inserting  equations  11-35  and 
11-40  into  equation  11-39.  Equation  11-38  can  then  be  expressed  as 


=  <?n  <  v„  >  n  ■ 
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The  explicit  dependence  on  the  Fenni  level  can  be  removed  from  this  expression  by 
substituting  for  n  from  equation  II- 1  to  obtain 


e„  =  <  v„  >  Nc  ■ 
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It  is  often  observed  that  the  capture  cross  section  of  a  particular  trap  has  a  temperature 
dependence  of  the  form  [28] 

r-E.  ^ 


o-n(T)  =  aii(c°)-ex  p 


V  kT  j 
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where  Ea  is  effectively  a  capture  energy  barrier  and  cr„( qo)  is  the  high  temperature  limit  of 
the  capture  cross  section.  Showing  this  explicit  temperature  dependence  of  an,  and  that  of 
the  Nc  and  <v„>,  equation  11-42  becomes 


=  7 '  crn  (co)  ■ 
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where  ;/  =  4 V 6 7Z2  M  cmeffk2h  3 .  The  temperature  dependence  of  one  more  tenn  is  still 

implicit  in  equation  11-44.  The  temperature  dependence  of  E,h  is  given  by 

Eth  =  Et0  -  aT.  11-45 

Thermodynamically,  Eth,  ET0,  and  a  are  the  respective  changes  in  Gibbs  free  energy, 
enthalpy,  and  entropy  due  to  the  change  in  the  occupation  of  the  trap  level  [28].  Upon 
substitution  of  equation  11-45  into  equation  11-44,  and  some  rearrangement,  the  explicit 
temperature  dependence  of  en  is  given  by 

e„ ( T)  =  r  ■  <Tn («0  ■{—)■  e^T2  ■  expf  ~(^° +^g)|  11-46 

\  kT  ) 

This  expression  is  critical  to  the  analysis  of  the  DLTS  data.  Both  capture  and  emission 
rate  constants  (often  simply  referred  to  as  capture  and  emission  rates)  have  units  of  s'1, 
and  the  specific  ways  in  which  they  describe  carrier  trapping  and  detrapping  in  the 
depletion  region  of  a  Schottky  diode  will  be  discussed  in  chapter  IV. 

Deep  traps  are  known  to  be  responsible  for  causing  significant  degradation  of 
electrical  device  perfonnance.  In  GaAs-based  electrical  devices,  for  instance,  deep  traps 
have  been  shown  to  be  responsible  for  I-V  hysteresis  loops,  “backgating  phenomena” 
(where  a  voltage  on  one  device  affects  the  current  on  another  device),  and  low  frequency 
oscillations  [13:187]. 

A  recombination  center  is  a  defect  for  which  both  cn  and  cp  are  large.  In  this  case, 
a  trapped  electron  is  likely  to  recombine  with  a  hole  before  it  has  the  chance  to  be  emitted 
back  into  the  conduction  band.  In  generation,  an  electron  is  trapped  from  the  valence 
band  and  subsequently  emitted  into  the  conduction  band.  A  trap  level  must  be  located 
near  the  middle  of  the  band  gap  to  be  an  efficient  generation  center. 
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A  specific  example  of  an  electron  trap  observed  by  DLTS  is  the  nitrogen  vacancy 
single  donor  level  located  60-80  meV  below  the  conduction  band  minimum  of  GaN. 
Polenta  et  al.  [31]  proposed  that  the  EDI  and  ED2  trap  levels  originate  from  the  auto- 
ionized  T2  state  of  the  nitrogen  vacancy.  Equation  11-46  was  applied  in  an  attempt  to 
model  these  two  overlapping  peaks,  and  the  fitting  suggested  ED  1  had  an  activation 
energy  of  about  60  meV.  Thus,  it  was  concluded  that  ED  1  corresponds  to  the  isolated 
Vn,  whereas  ED2  was  suggested  to  be  related  to  Vn  but  with  some  difference  in 
microscopic  configuration  causing  it  to  have  a  different  capture  cross  section  and 
activation  energy  [31]. 

2.4.4  Mobility  Degradation 

In  many  cases  of  interest,  the  degradation  of  mobility  can  be  described  as  a  linear 
function  of  fluence  by  the  equation  [25:138], 

//  1  =//01(l  +  6O),  11-47 

where  b  has  units  of  cm-.  In  analogy  to  the  carrier  removal  rate,  a  damage  constant  Kfl  = 
bl/j-o  may  be  defined  and  used  to  determine  the  expected  increase  in  // 1  with  fluence. 

The  degradation  in  mobility  due  to  radiation-induced  defects  is  expected  to  be  greatest  at 
low  temperature,  because  ionized  impurity  (or  defect)  scattering  is  the  dominant 
scattering  mechanism  in  semiconductors  at  low  temperature. 

2.4.5  Luminescence  Degradation 

Radiation  generally  causes  the  luminescence  intensity  of  a  semiconductor  to 
decrease  by  creating  defects  that  act  as  non-radiative  recombination  centers.  Previous 
investigations  of  luminescence  degradation  in  GaN  have  placed  more  emphasis  on  GaN- 
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based  LEDs  [67,  53]  than  on  epitaxial  layers.  The  radiation  response  of  GaN-based 
LEDs  is  of  practical  importance  and  the  luminescence  can  be  measured  with  less 
experimental  uncertainty  than  is  found  in  photoluminescence  and  cathodoluminescence 
measurements,  however,  LEDs  have  additional  degradation  mechanisms  than  those 
measured  directly  in  epitaxial  GaN  layers.  Several  different  expressions  for  damage 
constants  have  been  proposed  depending  on  the  device  design,  as  changes  in  series 
resistance  and  other  factors  can  affect  the  electroluminescence  of  LEDs  [25].  Gaudreau 
et  al.  [67]  reported  that  a  GaN-based  LED  was  about  two  orders  of  magnitude  more 
resistant  to  2  MeV  proton  radiation  than  GaAs  LEDs.  The  same  type  of  GaN-based  LED 
was  irradiated  with  2.5  MeV  electrons  and  the  damage  was  observed  to  disappear  after  16 
hours  of  room  temperature  annealing  [53].  The  photoluminescence  degradation  in  thin 
epitaxial  films  of  GaN  grown  by  MBE  was  characterized  with  2  MeV  proton  radiation  by 
Khanna  et  al.  [68].  Following  previous  work  in  GaAs,  the  dose  dependence  of  the 
photoluminescence  intensity  /  was  expressed  as 


-1  =  AO, 


11-48 


9 

where  Io  is  the  PL  intensity  before  irradiation,  O  is  the  fluence  in  units  of  cm'  ,  and  K  is 
the  damage  constant  in  units  of  cm  .  Equation  11-48  was  not  applicable  until  O  ~  10 
cm'”  For  lesser  values  of  O  some  luminescence  degradation  was  apparent,  but  the  change 
in  (A//  -  1 )  was  considered  to  be  within  the  margin  of  error  of  the  pre-irradiated 
measurement.  The  damage  constant  obtained  for  O  >10  cm'  was  K  =1.4x10'  cm”. 
This  damage  constant  was  detennined  for  the  dominant  donor  bound  exciton  transition  in 
the  PL  spectrum.  Comparing  with  a  damage  constant  pertaining  to  the  dominant 
excitonic  peak  in  neutron- irradiated  GaAs,  Khanna  et  al.  concluded  that  GaN  is  ~  100 
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times  more  tolerant  to  radiation  than  GaAs  with  regard  to  luminescence  degradation.  It 
does  not  appear  that  any  comparable  study  has  been  perfonned  on  AlGaN  to  date. 

2.5  Radiation  Effects  on  AixGai.xN  Devices 

The  radiation  damage  mechanisms  outlined  above  can  work  independently  or  in 
concert  to  affect  the  functionality  of  a  particular  semiconductor  device.  Because 
AlGaN/GaN  high  electron  mobility  transistors  (HEMTs)  have  shown  great  promise  for 
military  and  space  applications,  several  researchers  have  focused  on  the  radiation 
response  of  these  devices,  and  their  findings  are  discussed  below.  The  Schottky  barrier 
diode  (SBD)  forms  an  integral  part  of  the  HEMT  structure.  A  stand  alone  SBD  is  also 
useful  as  a  photodetector  or  an  electrical  rectifier.  An  overview  of  radiation  effects  on 
both  of  these  devices  is  presented  below. 

2.5.1  SBD  Response 

As  majority  carrier  devices,  SBDs  are  inherently  more  radiation  tolerant  than  p-n 
junctions  [27].  The  most  radiation  sensitive  parameter  in  SBDs  is  usually  reverse  leakage 
current,  caused  by  trap-assisted  carrier  generation  in  the  depletion  region  [69].  Radiation 
can  also  causes  a  decrease  in  Schottky  barrier  height  due  to  carrier  removal.  However, 
the  Schottky  barrier  height  of  real  devices  is  often  determined  by  surface  defects  [14,  16]. 
These  can  be  present  in  large  densities  and  relatively  insensitive  to  carrier  removal.  The 
SBD  ideality  factor,  «/,  also  increases  with  increases  in  radiation  dose  due  to  the  increase 
in  non-thennionic  currents.  Analytical  expressions  have  been  formulated  for  defect 
related  contributions  to  the  SBD  forward  and  reverse  currents  [70,  71]  and  defect 
concentrations  and  activation  energies  may  be  extracted  from  temperature  dependent  I-V 
data  [72].  This  technique  is  generally  not  as  reliable  as  DLTS,  however. 
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2.5.2  HEMT  Response 

The  primary  DC  performance  parameters  of  standard  GaN  field  effect  transistors 
(FETs)  that  decrease  with  radiation  displacement  damage  are  maximum  transconductance 
and  drain  saturation  current.  Additionally,  low-frequency  and  generation-recombination 
noise  have  both  been  observed  to  increase  with  increase  in  defect  concentration  [79].  In 
GaN-based  HEMTs  specifically,  radiation  damage  in  the  thin  n-doped  AlGaN  layer 
causes  a  decrease  in  the  two-dimensional  electron  gas  (2DEG)  concentration.  This  could 
be  caused  by  compensation  from  acceptor-type  defects  in  the  AlGaN  layer,  but  the  effect 
may  in  fact  be  dominated  by  a  radiation-induced  reduction  in  the  piezoelectric  field 
within  the  AlGaN  layer.  Bradley  et  al.  [73]  correlated  the  reduction  of  as-grown  2DEG 
carrier  density  in  certain  MBE-grown  HEMTs  with  an  Alo.3Gao.7N  defect  luminescence 
peak  at  2.34  eV.  Cai  et  al.  [74]  reported  a  68%  and  62%  reduction  in  the 
transconductance  and  drain  saturation  current  of  an  MBE-grown  Alo.15Gao.85N/GaN 
HEMT  following  irradiation  with  1.0x10  protons/cnr  at  1.8  MeV.  After  annealing  for 
40  seconds  at  800  °C,  the  transconductance  and  drain  saturation  current  returned  to  70% 
and  85%  of  the  respective  pre-irradiation  values.  Hall  measurements  revealed  that  an 
increase  in  electron  carrier  concentration  was  responsible  for  the  recovery,  whereas 
mobility  decreased  after  the  annealing.  Luo  et  al.  [75]  performed  a  radiation  study  on 
similar  HEMTs  of  varying  channel  geometries  which  were  fabricated  by  SVT  Associates 
[76].  The  devices  were  irradiated  with  40  MeV  protons  from  5xl09  to  5xl010  cm'2, 
followed  by  isothennal  annealing  at  300  °C.  Transconductance  decreased  by  less  than 
15%  for  a  dose  of  5xl09  cm'2.  The  authors  attributed  this  to  an  increase  in  channel 
resistance.  Additionally,  Luo  et  al.  found  that  reverse  bias  breakdown  voltage  increased 
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in  magnitude  with  dose,  suggesting  that  carrier  removal  dominated  the  rise  in  resistivity. 
The  decrease  in  saturation  current  was  30-50%  for  the  highest  dose  of  5xl010  cm'2.  They 
reported  a  significant  increase  in  both  transconductance  and  saturation  current  following 
60  seconds  of  annealing  at  300  °C.  There  was  no  systematic  effect  of  gate  width  or 
length  on  the  HEMT  parameters.  Gaudreau  et  al.  [67]  recently  reported  temperature 
dependent  Hall  effect  measurements  on  2  MeV  proton-irradiated  MBE-grown 
Alo.3Gao.7N/GaN  2DEG  test  structures.  Between  a  fluence  of  3xl014  and  3xl015  cm'2,  the 
structures  changed  from  being  a  conductor  to  being  an  insulator.  The  sheet  charge 
density  dropped  to  about  50%  of  its  original  value  and  the  mobility  dropped  to  about  20% 
of  its  original  value.  A  further  increase  in  resistivity  beyond  a  fluence  of  3x10 15  cm'2  was 
due  to  a  large  drop  in  mobility,  whereas  sheet  charge  density  remained  mostly 
unchanged. 


50 


III.  Experimental  Procedures 


This  chapter  describes  the  procedures  used  to  obtain  the  experimental  data 
presented  in  chapter  V.  The  reasons  for  choosing  the  given  test  materials  are  discussed 
first,  followed  by  a  description  of  the  preparation  of  test  samples.  Next  are  described  the 
procedures  by  which  electron  radiation  doses  were  applied  to  the  samples.  Finally,  the 
three  main  characterization  techniques;  DLTS,  Hall  Effect,  and  CL  are  introduced  along 
with  specific  procedures  followed  when  applying  these  characterization  techniques. 

3.1  Sample  Selection  and  Preparation 

The  AlxGai_xN  material  used  in  this  study  was  purchased  from  SVT  Associates 
[76].  Previous  experience  proved  that  this  company  was  capable  of  producing  good 
quality  GaN  for  electrical  and  optical  characterization  as  well  as  good  quality  AlxGai_xN 
(x  <  0.3)  for  optical  characterization.  The  1pm  thick  AlxGai_xN  layers  were  grown  by 
plasma-assisted  MBE  on  the  c-plane  of  0.5  mm  thick  polished  sapphire  substrates  with 
low-temperature  grown  AIN  buffer  layers  between  the  two.  Each  wafer  was  custom 
ordered  with  a  requested  aluminum  mole  fraction  and  silicon  doping  level.  Before 
choosing  the  requested  silicon  doping  levels,  a  simple  estimate  of  room  temperature 
carrier  concentration  was  made  by  solving  the  charge  neutrality  equation  with  the 
assumption  of  a  single  hydrogenic  donor  level  and  a  given  background  concentration  of 
acceptors  [77], 
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Here,  n  is  the  carrier  concentration,  ND  is  the  intentionally  doped  donor  concentration,  Na 
is  the  unintentionally  doped  acceptor  concentration,  and  g  is  the  degeneracy,  which  is 
equal  to  2.  ED  is  the  donor  ionization  energy  given  by 


where  m*  and  er  are  the  electron  effective  mass  and  static  dielectric  constant  of  AlxGai_xN 
for  the  given  A1  mole  fraction.  From  equation  III-2,  the  er  values  in  table  II- 1,  and  the 
effective  mass  m  =  0.22  mQ  the  donor  ionization  energies  are  38,  48,  and  53  meV  for  x  = 
0,  0.2,  and  0.3,  respectively.  Theory  indicates  that  the  Vc,a  acceptor  fonnation  energy 
decreases  as  increases  in  GaN  [29].  Thus,  more  VGa  acceptors  are  formed  during 
growth  of  GaN  that  is  more  n-type.  Another  source  of  compensation  is  Si  atoms  that 
reside  on  N  lattice  sites,  acting  as  acceptors  instead  of  donors.  Thus,  Na  is  not 
independent  of  No,  and  can  be  approximated  as  being  a  constant  fraction,  called  the 
compensation  ratio  of  No.  This  constant  fraction  is  typically  equal  to  0.6  ±  0.1  [77,78]  . 
Using  this  compensation  ratio  and  No  =10  cm" ,  equation  III-2  yields  a  295  K  carrier 
concentration  of  n  =  2.9xl016  and  2.5xl016  cm"3  for  x  =  0  and  0.3,  respectively.  In  both 
cases,  then,  the  donor  level  is  sufficiently  shallow  for  nearly  complete  ionization  at  room 
temperature.  At  77  K,  the  carrier  concentrations  are  3.2xl014  and  3.4xl013  cm'3  for  x  =  0 
and  0.3,  respectively.  Actual  silicon  ionization  energies  are  generally  shallower  (for  x  < 
0.3)  than  those  predicted  by  the  simple  theory  above  due,  in  part,  to  free  carrier 
screening.  For  x  =  0  and  Nd=  3x10  cm'  ,  it  is  found  experimentally  that  Eo  =  1 7  meV 
[41].  For  x  =  0.3  and  ND=  4xl017  cm'3,  ED  =  20  meV  [79].  The  effect  of  screening  upon 
donor  ionization  energy  can  be  approximated  by 
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III-3 


ED(ND)  =  ED0-aNDm, 

where  the  screening  factor,  a  ,  was  empirically  determined  to  be  2. 1  x  1  O  '  meV-cm  for  Si 
in  GaN  [80].  A  screening  factor  of  2.5xl0‘5  meV-cm  has  been  reported  for  Alo.23Gao.77N 
[81]. 

The  actual  measured  ionized  carrier  concentrations  (for  x  >  0)  proved  to  be 
somewhat  unpredictable.  One  reason  for  this  was  that  the  carrier  concentration  measured 
by  the  growers  for  quality  assurance  was  artificially  high  due  to  an  impurity  conduction 
layer  near  the  sapphire  substrate  interface.  However,  it  was  also  apparent  that  precise 
control  of  silicon  doping  concentration  was  fundamentally  more  difficult  for  x  >  0  than  it 
was  for  x  =  0.  The  aluminum  mole  fraction  of  the  wafers  was  reported  by  SVT 
Associates  based  on  CL  and  X-ray  diffraction  measurements.  The  CL  measurements 
performed  at  AFIT  showed  a  small  variation  in  the  radial  direction  such  that  aluminum 
mole  fractions  were  greatest  at  the  center  of  each  wafer.  A  2  pm  thick  GaAs  epitaxial 
layer  doped  at  6.0xl017  cm"'  with  silicon  was  also  included  in  this  study.  This  epitaxial 
layer  was  grown  by  MBE  on  a  GaAs  substrate  at  the  Air  Force  Research  Laboratory. 

The  purpose  for  including  this  wafer  was  to  compare  the  radiation  hardness  of  AlxGai_xN 
with  the  better  known  GaAs  grown  by  a  comparable  epitaxial  technique.  An 
Alo.27Gao.73N/GaN  single  heterostructure  grown  by  Cree,  Inc.  [82]  was  also  included  in 
this  study  to  compare  the  effects  of  radiation  doses  on  the  heterostructure  transport 
properties.  Table  III- 1  lists  the  alpha-numeric  designations  of  the  wafers  used  in  this 
study  as  well  as  their  measured  properties.  The  third  and  forth  columns  of  this  table 
show  the  apparent  carrier  concentration  and  mobility  detennined  by  room  temperature 
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Table  III-l.  Wafer  properties  and  designations. 


I.D. 

Material  Type 

Apparent 

Carrier 

Cone. 

(e'/cm3) 

Apparent 

Mobility 

(cm2/V-s) 

Corrected 

Carrier 

Cone. 

(e/cm3) 

Corrected 

Mobility 

(cm2/V -s) 

CL 

Peak 
Energy 
(eV)  ' 

A1 

AlxGai_xN  x  =  0 

0.65  xlO17 

139 

0.65xl017 

140 

3.49 

A2 

AlxGai_xN  x  =  0 

1.6  xlO17 

324 

1.6xl017 

324 

3.49 

B1 

AlxGai.xN  x  =  0.10 

2.3  xlO17 

78 

<  1.0  x  1017 

3.71 

B2 

AlxGai.xN  x  =  0.14 

2.7  xlO17 

56 

8.2  xlO16 

64 

3.87 

Cl 

AlxGai.xN  x  =  0.20 

6.5  x  1017 

44.3 

3.0  x  1017 

49.4 

4.00 

C2 

AlxGai.xN  x  =  0.20 

1.0  x  1017 

13.3 

0.8  x  1017 

14.6 

4.00 

D1 

AlxGai.xN  x  =  0.30 

2xl017 

25-35 

lxlO17 

30-45 

4.26 

D2 

AlxGai.xN  x  =  0.30 

8.6xl017 

15.7 

7-8  xlO17 

16.1 

4.26 

El 

Alo.27Gao.73N 

/GaN 

Sheet  cone. 
l.lxlO13  cm-2 

1410 

FI 

GaAs 

6.4xl017 

2980 

1.52 

Hall  measurement  of  a  representative  sample  from  the  wafer.  The  fifth  and  sixth 
columns  show  the  corresponding  carrier  concentration  and  mobility  corrected  for 
interfacial  layer  conduction  as  described  in  section  3.4.2  below. 

The  last  column  of  table  III- 1  lists  the  T  =  6  K  peak  CL  energies  of  the  samples 
used  in  this  work,  and  these  energies  are  plotted  as  data  points  in  figure  III- 1  to  compare 
with  the  predicted  band  gap  energy  assuming  different  bowing  parameters.  It  is  seen  that 
a  bowing  parameter,  b  =  0.35,  provides  the  best  fit  to  the  plot  of  band  gap  versus 
aluminum  mole  fraction  x  in  these  samples. 
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Figure  III-l.  Measured  CL  peak  energies  (black  squares)  of  AlxGai.xN  samples  versus 
x.  The  CL  measurements  are  performed  at  T  =  6  K. 


Most  sample  preparation  including  cleaning,  masking,  metal  deposition,  and 
contact  annealing  was  performed  in  the  AFIT  clean  room.  Each  2”  diameter  wafer 
arrived  with  a  thin  titanium  backing  that  had  to  be  etched  away  in  a  hydrofluoric  acid 
solution.  The  titanium  backing  had  been  necessary  to  distribute  heat  across  the  substrate 
during  MBE  growth.  The  wafer  was  then  cut  into  5x5  mm  squares  (die)  with  a  diamond 
blade  saw.  A  number  was  marked  on  the  back  of  each  sample  with  a  diamond  scribe  in 
order  to  keep  track  of  the  sample’s  original  location  on  the  wafer.  The  oxidized  layer 
expected  to  be  on  the  surface  of  each  sample  was  removed  by  immersion  in  boiling  aqua 
regia  (1:3  HNOsiHCl)  for  2-4  minutes.  Ohmic  contact  areas  were  then  defined  by 
photolithography  (discussed  below)  for  the  DLTS  samples  and  by  a  perforated  stainless 
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steel  mask  for  the  Hall  samples.  The  Ohmic  metal  was  deposited  by  electron  beam 
evaporation.  The  Ohmic  metallization  used  for  GaN  samples  was  Ti(400A)/Al(2000A), 
where  the  notation  indicates  that  400  A  of  titanium  is  deposited  first,  followed  by  2000  A 
of  aluminum.  The  standard  Ohmic  metallization  used  for  AlGaN  samples  was  Ti  (300 
A)/A1  (800  A)/Ti  (1200  A)/Au  (500  A).  Following  Ohmic  metal  deposition  (and 
photoresist  lift-off  in  the  case  of  the  DLTS  samples),  the  samples  were  annealed  at  900 
°C  for  45  s  in  a  nitrogen  gas  ambient.  This  step  was  perfonned  in  the  rapid  thermal 
annealing  furnace  located  in  the  AFIT  clean  room.  At  this  point,  the  Hall  samples  were 
ready  for  characterization.  The  DLTS  samples  required  the  additional  steps  involved  in 
making  Schottky  barrier  contacts.  As  with  the  Ohmic  contact  areas,  the  circular  Schottky 
contact  areas  (diameters  of  200,  300,  and  400  pm)  were  defined  via  photolithography.  In 
this  process,  1813  photoresist  was  spread  into  a  uniform  layer  on  the  surface  of  the 
sample  at  4000  rpm.  The  sample  was  then  placed  on  a  100  °C  hot  plate  for  5  minutes  in 
order  to  dry  the  photoresist.  The  sample  was  then  placed  in  contact  with  the 
photographic  mask  in  the  mask  aligner  and  exposed  to  UV  light  for  20  s.  The  exposed 
areas  of  the  sample  were  then  removed  by  immersing  the  sample  in  35 1  developer 
solution  for  30  s.  The  sample  was  then  rinsed  in  de-ionized  water  (DIW)  for  30  s 
followed  by  a  30  s  dip  in  dilute  hydrochloric  acid  solution  (1:10  HCkDIW)  and 
subsequent  DIW  rinse.  The  dilute  HC1  is  used  to  remove  surface  oxidation  that  grew 
since  the  previous  oxidation  removal  step.  Following  every  DIW  rinse,  the  sample  was 
blown  dry  with  nitrogen  gas.  At  this  point,  the  sample  had  a  photoresist-defined  pattern 
of  exposed  circular  areas  that  were  clean  and  de-oxidized.  The  sample  was  then  placed  in 
the  electron  beam  evaporator  such  that  the  entire  masked  surface  was  exposed  to  metal 
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deposition.  The  metallization  used  for  Schottky  barrier  contacts  was  Ni  (350  A)/Au 
(4000  A).  Following  metallization,  lift-off  of  the  photoresist  was  accomplished  by 
placing  the  sample  in  acetone  and  exposing  it  to  ultrasonic  agitation  for  10-20  minutes. 

At  this  point  the  sample  could  be  diced  and  packaged  for  DLTS  measurements.  The  5x5 
mm  samples  were  typically  diced  into  four  square  quarters  using  a  diamond-tip  scribing 
tool.  This  step  yielded  samples  small  enough  to  fit  onto  the  TO-5  cans  used  for  DLTS 
characterization.  A  manual  ultrasonic  gold  wire  bonder  was  used  to  make  a  gold  wire 
connection  between  one  of  the  circular  Schottky  contacts  on  the  sample  and  one  of  the 
gold  posts  on  the  can.  The  wire  bonder  was  also  used  to  make  a  gold  wire  connection 
between  the  Ohmic  area  on  the  sample  and  another  gold  post  on  the  can.  A  complication 
was  encountered  in  this  step,  as  it  was  found  that  the  wire  bond  would  not  adhere  to  the 
annealed  Ohmic  contact  area  directly.  Consequently,  the  procedure  was  altered  such  that 
Au  was  deposited  to  connect  at  least  one  of  the  Schottky  pads  and  the  Ohmic  contact 
area.  This  step  was  accomplished  by  masking  off  the  appropriate  area  with  rubber 
cement  or  the  sticky  part  of  a  Post-it™  note,  and  then  sputtering  gold  onto  the  surface  in  a 
table-top  sputterer.  After  this  step,  the  exposed  Schottky  pad  served  as  an  Ohmic  contact. 

3.2  Irradiation  Procedures 

The  electron  irradiations  were  perfonned  at  the  Wright  State  University  Van  de 
Graff  (VDG)  generator  facility.  The  available  electron  beam  energy  range  is 
approximately  0.40-1.8  MeV.  The  available  beam  current  depends  on  the  chosen  energy. 
Typical  beam  currents  that  allow  for  stable  operation  fall  in  the  range  of  1-30  pA.  At  the 
lowest  achievable  energies,  the  beam  current  drops  off  considerably,  and  the  beam  energy 
tends  to  be  unstable.  A  current  integrator  attached  to  the  cold  head  keeps  a  running  count 
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of  the  negative  charge  impinging  on  the  cold  head.  The  total  fluence  is  then  given  with 
explicit  dimensional  analysis  as 

0  =  Qtot  Coulomb  1 _ =  Qtot  e~  m4 

A  cm 2  1.602x10  19  Coulomb  1.602xlCTl9A  cm 1 

where  A  is  the  irradiated  area,  and  Qtot  is  the  total  integrated  charge.  Before  January  of 

2003,  this  area  was  approximately  equal  to  3.3  cm'.  In  January  of  2003,  however,  a  two 

inch  spacer  tube  was  inserted  between  the  aperture  and  the  cold  head.  The  beam  focal 

point  is  in  front  of  the  aperture.  Thus,  the  beam  diverges  after  this  point  in  the  beam  line. 

It  was  found,  by  irradiating  a  plastic  sheet  at  the  plane  of  the  cold  head,  that  the  new 

beam  area  was  approximately  7.9  cm",  and  this  is  the  number  used  in  calculating  the 

fluence  in  all  the  doses  of  table  III-2,  with  the  exception  of  the  dose  applied  on  30  Apr 

04.  At  that  time,  the  samples  were  mounted  on  a  different  cold  head,  and  this  other  cold 

head  was  further  from  the  aperture  than  the  normal  cold  head.  The  divergence  angle  was 

determined  from  the  previous  spot  size  measurement  to  be  approximately  equal  to  0.25 

radians.  This  angle  was  used  to  determine  that  the  beam  area  at  the  plane  of  this  second 

cold  head  was  18.6  cm  .  Because  it  took  several  weeks  to  complete  the  post-irradiation 

characterization,  the  samples  were  stored  in  liquid  nitrogen  following  irradiation.  This 

reduced  the  risk  of  obtaining  conflicting  results  due  to  long-term  room  temperature 

annealing.  All  of  the  samples  were  exposed  to  at  least  3  hours  of  room  temperature 

annealing  before  being  characterized.  Table  III-2  lists  the  parameters  pertaining  to 

irradiations  performed  in  the  course  of  this  work. 
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Table  III-2.  Irradiation  parameters. 


Irradiation 

Date/ 

Facility 

Energy/ 

Type 

Beam 

Current 

Irradiated 

Area 

Effective  Beam 
Flux 

Total  Fluence 

20  Oct  03 
WSU 

1.0  MeV 
electrons 

15  pA 

7.9  cm2 

lxlO13  ctrfV 

9x10 16  cm'2 

30  Apr  04 
WSU 

1.0  MeV 
electrons 

10  pA 

18.5  cm2 

6  xlO12  cm'V 

2x10 16  cm'2 

9  Jun  04 
WSU 

1.0  MeV 
electrons 

10  pA 

7.9  cm2 

8  xlO12  ctrfV 

9x1 016  cm'2 

31  Aug  04 
WSU 

0.62  MeV 
electrons 

1.5  pA 

7.9  cm2 

8  xlO12  ctrfV 

2.5xl016  cm'2 

3.3  Deep  Level  Transient  Spectroscopy  (DLTS) 

The  core  components  of  the  DLTS  system  used  in  this  dissertation  research 
(figure  III-2)  consist  of  an  evacuated  JANIS  Corp.  sample  chamber  with  closed  cycle 
refrigerated  helium  cold  head  and  coaxial  electrical  feedthroughs,  a  Lake  Shore  330 
temperature  controller,  a  SULA  Corp.  capacitance  meter,  which  supplies  a  100  mV  test 
waveform  at  1  MHz  to  measure  capacitance,  a  Lecroy  1  GHz  pulse  generator  for  system 
timing  signals,  and  a  PC  and  data  acquisition  board,  which  supply  all  pulsed  and  steady- 
state  bias  voltages  to  the  diode  and  control  the  entire  data  acquisition  process  in  an 
automated  manner.  The  temperature  range  of  the  sample  chamber  is  50  to  450  K,  and 
temperature  control  is  good  to  ±  0.1  K  in  most  temperature  ranges.  All  measurements  in 
this  study  are  made  on  a  100  pF  capacitance  range  setting.  The  auxiliary  equipment 
includes  a  Keithley  237  Source  Measurement  Unit  (SMU),  which  is  used  for  I-V 
measurements,  and  a  Lecroy  9410  oscilloscope,  which  is  used  to  monitor  diode  input  and 
output  voltage  waveforms.  Additionally,  a  low-pass  fdter  with  a  cutoff  frequency  of 
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Figure  III-2.  Schematic  diagram  of  DLTS  system  used  in  this  study.  The  Keithly  237 
SMU  is  used  for  I-V  measurements  (after  Scofield  [72]) 
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slightly  greater  than  1  MHz  was  inserted  in  the  signal  path  for  improved  signal-to-noise 
ratio. 

The  wire-bonded  SBD  is  placed  in  good  thennal  contact  with  the  cold  head  and 
fine  coaxial  wires  are  connected  to  the  Ohmic  and  Schottky  leads  that  protrude  from  the 
TO-5  can.  These  wires  run  to  a  coaxial  feedthrough  outside  of  which  BNC  cables  run  to 
the  capacitance  meter.  Before  beginning  the  DLTS  temperature  scan,  room  temperature 
capacitance  vs.  voltage  (C-V)  and  current  vs.  voltage  (I-V)  measurements  were 
performed  at  room  temperature.  These  measurements,  which  were  also  perfonned  at 
about  20  K  increments  throughout  the  50-450  K  temperature  range,  are  essential 
diagnostic  tools  for  determining  the  quality  of  the  SBD,  and  thus  the  suitability  for  DLTS 
measurements.  Additionally,  the  computer  control  program  that  makes  the  C-V 
measurement  also  determines  the  ionized  shallow  donor  concentration  No  from  the  slope 
of  the  1/C'  vs.  V  curve.  This  calculation  is  based  on  the  relationship  [83] 

N°  ~  e£sA2[d{^IC2)l dv]  m"5 

where  es  is  the  static  dielectric  constant  and  A  is  the  area  of  the  Schottky  contact.  These 
two  values  are  given  as  inputs  to  the  program. 

The  data  which  are  directly  obtained  by  the  DLTS  procedure  are  digitized 
capacitance  transients  recorded  at  2  K  temperature  intervals  across  a  certain  temperature 
range.  Capacitance  transients  are  depicted  in  figure  III-3.  The  bias  voltage  applied  to  the 
SBD  during  DLTS  measurement  is  depicted  in  the  upper  section  of  figure  III-3.  The 
SBD  is  initially  held  at  a  reverse  bias  voltage  of  magnitude  Vr .  The  capacitance 
corresponding  to  the  SBD  thus  biased,  C0  is  as  shown  in  the  lower  section  of  figure  III-3. 
This  static  situation  is  interrupted  by  a  forward  bias  voltage  pulse  Vfof  duration  tp.  This 
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Figure  III-3.  Applied  voltage  pulses  and  corresponding  capacitance  transients  in  the 
DLTS  measurement  procedure. 


Vf  pulse  causes  the  depletion  region  of  the  SBD  to  collapse  and  capacitance  C  to  spike 
upward  as  shown.  During  this  time,  the  introduction  of  majority  carriers  into  the 
previously  depleted  region  allows  previously  empty  majority  carrier  traps  to  be  fdled. 
When  the  Vf  pulse  is  removed,  the  depletion  region  returns  to  the  usual  width  W  plus  an 
extra  amount  A  IF  in  response  to  the  trapped  majority  carrier  charge.  This  extra  depletion 
width  corresponds  to  a  decrease  in  SBD  capacitance  indicated  by  AC  in  figure  III-3.  At 
reverse  bias,  the  majority  carrier  traps  emit  the  trapped  carriers  in  an  exponentially 
decaying  manner  with  time  constant  en  described  in  section  2.4.3.  The  emission 
corresponds  to  the  exponential  decay  of  A C{t)  to  the  quasi-equilibrium  value  of  C0.  This 
capacitance  transient  is  expressed  as 

AC(t)  =  AC(1  -  exp(-e„t)),  III-6 

where  t  =  0  at  the  end  of  the  filling  pulse.  In  order  to  increase  the  signal-to-noise  ratio  in 
the  measured  capacitance  transient,  500-800  transient  signals  are  collected  and  averaged 
at  each  temperature  step.  The  recorded  capacitance  transients  may  then  be  analyzed 
according  to  the  methods  of  section  4. 1 . 
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One  other  way  in  which  capacitance  transients  may  be  obtained  is  by  double-correlated 
DLTS  (DDLTS).  The  DDLTS  technique  simply  subtracts  the  capacitance  transient 
obtained  with  a  smaller  magnitude  of  forward  voltage  pulse  Vn  from  the  capacitance 
transient  obtained  with  a  larger  magnitude  of  forward  voltage  pulse  Vn .  The  resulting 
DDLTS  capacitance  transient  represents  the  carrier  trapping  and  emission  in  a  narrow 
width  of  the  depletion  region  defined  by  the  choices  of  Vn  and  Vn- 

3.4  Hall  Effect  Measurements 

3.4.1  Room  Temperature  Hall  Effect  Measurements 

The  Hall  effect  has  long  been  the  standard  method  for  measurement  of  free  carrier 
concentration  and  mobility  as  well  as  resistivity.  Figure  III-4  shows  a  diagram  of  the 
basic  components  of  the  Lake  Shore  HMS  7700  Hall  effect  measurement  system  used  in 
this  study.  As  usual  for  measurement  of  epitaxial  layers,  the  van  der  Pauw  technique  is 
used  to  measure  resistivity  and  Hall  voltage.  In  this  technique,  four  Ohmic  contacts  are 
placed  on  the  periphery  of  a  sample  as  illustrated  in  figure  III-5.  The  Ohmic 
metallization  layers  were  deposited  and  annealed  as  described  previously.  In  order  to 
make  electrical  connection  between  the  sample  and  the  Hall  measurement  system,  the 
sample  was  affixed  with  rubber  cement  to  a  removable  sample  mount,  indium  was 
soldered  to  the  Ohmic  contacts,  and  the  four  wires  from  the  sample  mount  were  then 
easily  soldered  to  the  indium. 

The  physical  concept  of  the  Hall  measurement  will  be  described  here,  briefly.  A 
fixed  current,  /,  is  sourced  between  two  opposing  (across  the  diagonal)  contacts  under  the 
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LakeShore  HMS  7700 


Figure  III-4.  Schematic  of  Temperature  Dependent  Hall  measurement  system. 
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Figure  III-5.  Illustration  of  Hall  effect  in  planar  sample  with  Van  der  Pauw  contact 
geometry  [15]. 

influence  of  a  perpendicularly  directed  magnetic  field,  B  .  In  this  situation,  the  carriers 
experience  the  Lorentz  force,  given  as 

F  =  q(vx  B),  III-7 

where  q  is  -e  for  electrons  and  +e  for  holes.  Because  electrons  and  holes  have  opposite 
charge  and  move  in  opposite  directions  under  the  applied  electric  field,  the  Lorentz  force 
causes  both  carriers  to  be  displaced  in  the  same  direction.  The  dominant  carrier  type 
causes  a  net  dipole  moment  to  result  from  this  displacement,  and  it  is  measured  as  a 
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voltage  between  the  opposing  contacts.  In  a  simplified  Hall  bar  geometry,  this  Hall 
voltage  would  be  given  by  the  expression  [15] 


/  Br 

VH  =  —  —  =  R-hsI  xBrH .  III-8 

ens 

2  2 

In  this  expression,  Rhs  is  the  sheet  Hall  coefficient  and  rn  (  =  <r  >/<r>  )  is  the  Hall 
scattering  factor,  defined  by  x,  which  is  the  mean  time  between  scattering  determined  by 
a  variety  of  scattering  mechanisms  [83].  Rhs  can  be  expressed  as 

RHs= HI-9 
IxBrH 


Assuming  n-type  material ,  the  sheet  carrier  concentration,  ns,  and  the  electron  mobility, 
/ 1„ ,  are  defined  in  terms  of  this  coefficient  as  [83] 


=  -- 


eR, 


III- 10 


and 


Pn 


rHps 


III- 11 


The  Hall  scattering  factor  can  theoretically  have  any  value  in  the  range  1.0  <  rn<  2.0, 
and  1.2  would  be  a  typical  values  for  many  semiconductor  samples.  Because  the 
individual  scattering  mechanisms  are  not  characterized,  th  is  assumed  to  be  equal  to  1.0. 
Thus,  the  quantities  obtained  from  the  Hall  measurement  are  understood  to  be  the  Hall 
carrier  concentration  and  Hall  mobility,  which  overestimates  the  actual  carrier 
concentration  and  underestimates  the  actual  mobility  by  the  factor  th.  In  the  van  der 
Pauw  geometry  of  figure  III-5,  ps  in  equation  III- 1 1  is  determined  from  a  total  of  eight 
different  current-voltage  pair  measurements,  and  then  four  different  current-voltage  pair 
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measurements  are  made  under  the  application  of  the  magnetic  field  to  obtain  Rhs  in 
equation  III- 1 0  and  III- 11.  In  the  implementation  here,  the  electromagnet  current  is 
reversed  and  Rh  is  also  obtained  under  the  reversed  magnetic  field.  The  average  of  the 
two  Rh  values  is  then  used  to  calculate  ns  and  jjn.  These  two  values  are  calculated  by  the 
measurement  system  software.  The  volume  carrier  concentration  is  then  obtained  from 
the  sheet  concentration  by 

n  =  III- 12 

d 

where  d  is  the  thickness  of  the  epitaxial  layer.  The  source  current  used  in  this  work 
ranged  from  10'  to  10'  A,  depending  on  the  resistivity  of  the  sample.  Specifically,  the 
current  was  chosen  in  order  to  keep  the  Hall  voltage  within  the  desired  range  of  10  to  500 
mV. 

3.4.2  Temperature-Dependent  Hall  Effect  Measurements 

In  the  case  of  temperature  dependent  Hall  effect  measurements,  the  above 
measurements  and  calculations  are  performed  at  different  temperatures  under  automated 
computer  control.  The  chamber  is  evacuated  and  an  air-cooled  helium  compressor 
enables  the  cold  head  refrigerator  to  achieve  20  K  sample  temperature.  The  automated 
temperature  control  system  then  raises  the  temperature  in  programmed  steps.  Steps  of 
2.5  K  were  used  when  20  K  <  T<  70  K,  and  then  5  K  steps  were  used  when  70  K  <  T< 
320  K.  A  wait  time  of  seven  minutes  was  used  at  each  temperature  step  to  give  the 
sample  time  to  equilibrate  with  the  sample  holder.  For  shielding  consistency  and  to  avoid 
persistent  photoconductivity  effects,  all  of  the  measurements  were  performed  in  the  dark, 
inside  the  vacuum  chamber. 
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Because  epitaxial  layers  of  GaN  and  AlGaN  are  often  observed  to  have  a 
degenerate  donor  conduction  layer  near  the  substrate  [45,  46],  low  temperature  Hall 
effect  measurements  were  necessary  in  order  to  deconvolve  the  conductivity  contribution 
of  the  epitaxial  layer  from  the  conductivity  contribution  of  the  degenerate  near-interfacial 
conduction  layer.  Following  the  deconvolution  method  outlined  by  Look  [33],  the 
measured  sheet  carrier  density  ns  and  mobility  fx  at  T  =  20  K  were  taken  to  be  equal  to 
that  of  the  sheet  carrier  density  nsi  and  mobility  /xi  of  the  near-interfacial  conduction 
layer.  This  step  was  based  on  the  assumption  that  the  conduction  in  the  near-interfacial 
layer  is  independent  of  temperature  and  that  the  epitaxial  layer  sheet  carrier  concentration 
nsi  and  mobility  jui  are  very  low  at  T  =  20  K.  With  these  assumptions,  a  two-layer 
conduction  analysis  [17]  may  be  applied  in  order  to  extract  the  volume  density  n\  and 
mobility /o  at  temperatures  higher  than  20  K.  This  extraction  is  accomplished  at  each 
temperature  step  by  the  equations 


Ai  = 


ju~n 


//  n 


t^2nS2 

-ju2nS2ld 


III- 13 


and 


(//  n  -  ju2ns2 1  d)2 
ju2n  -/u;nS2/d 


III- 14 


where  d  is  the  epitaxial  layer  thickness. 


3.5  Cathodoluminescence 

A  diagram  of  the  cathodoluminescence  experiment  is  shown  in  figure  III-6.  The 
CL  measurements  were  performed  at  6  K  under  vacuum.  The  excitation  source  was  a 
Kimball  Physics  EMG-12  electron  gun.  In  this  dissertation  research,  a  beam  energy  of  5 
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Figure  III-6.  Schematic  of  cathodoluminescence  system  [84]. 
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keV  was  chosen  with  a  beam  current  of  25  jli  A  and  focused  spot  size  of  approximately  2 
mm  diameter.  The  light  emitted  from  the  sample  was  dispersed  by  a  SPEX  500M 
spectrometer  with  a  0.2  A  resolution  and  directed  into  a  liquid  nitrogen-cooled 
photomultiplier  tube  (PMT)  biased  at  1500  kV.  The  PMT  signal  was  recorded  by  the 
control  computer  at  each  wavelength  step  of  the  spectrometer.  A  typical  scan  of  the 
spectrometer  started  at  2800  A  and  moved  to  8000  A  in  2  A  steps.  The  chosen 
integration  time  was  0.1s.  CL  characterization  was  perfonned  on  the  Hall  effect  samples, 
and  grounding  wires  were  soldered  to  one  Ohmic  contact  on  each  sample.  The  other  ends 
of  these  wires  were  tied  to  the  copper  sample  holder  to  prevent  excess  charge  build-up  on 
the  sample  surface  during  measurement. 
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IV.  Analytical  Procedures 


4.1  Deep  Level  Transient  Spectroscopy  (DLTS) 


4.1 .1  Capacitance  T ransient  Analysis 

In  this  study,  the  digitized  capacitance  (C)  transient  data  are  analyzed  according 
to  the  “boxcar”  method  (figure  IV- 1)  whereby  C  is  extracted  at  two  different  times,  t\  and 
E.  The  difference  SC  =  C(/|)  -  C(h)  is  a  negative  value  for  majority  carrier  trapping  and 
a  positive  value  for  minority  carrier  trapping  whether  the  material  is  n-type  or  p-type.  In 
the  remainder  of  this  discussion  the  behavior  of  a  majority  carrier  trap  (specifically, 
electrons  in  n-type  material)  will  be  assumed,  but  the  arguments  are  equally  valid  for 
minority  carrier  traps.  As  temperature  is  swept,  a  positive  peak  occurs  in  the  signal  - 
6C(T)  at  temperature  Tm  (figure  IV-1),  and  at  this  temperature  it  can  be  shown  [83]  that 


en(T,„) 


li aMi) 

t2  —  tx 


IV-1 


This  emission  rate  constant,  determined  by  the  choice  of  t\  and  V  is  referred  to  as  the 
“rate  window.”  Different  combinations  of  t\  and  h  are  chosen  that  will  yield  sufficiently 
different  en.  Taking  the  natural  logarithm  of  both  sides  of  equation  11-46  and  rearranging 
tenns  yields 


=  In 


/•cr„  (go) 


«V|  ,  (Etq  +EJ  1000 

kjj  1000k  T 


IV-2 


Equation  IV-2  shows  that  a  plot  of  1  nieJT1)  versus  1 000/ T  (an  “Arrhenius  plot”)  will 
produce  a  straight  line  having  a  slope  equal  to  (ET0  +  En)l  1 000k  and  a  T— >  oo  intercept 
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dependent  capacitance  transients  (after  Lang  [85]). 


equal  to  ln[}>  (gjg\)(a/k)  cr„(oo)].  Accordingly,  an  Arrhenius  plot  is  constructed  from  the 
measured  e„(Tm)  data  and  a  linear  fit  reveals  the  total  activation  energy,  ETOT  (=  ETo+Ea  ), 
and  the  apparent  capture  cross  section,  ona  (=  (, gjg\)(alk )  an(co)).  These  two  values  are 
often  referred  to  as  the  “DLTS  signature”  of  a  trap  level,  and  are  often  reported  in  the 
electrical  characterization  of  deep  trap  defects.  The  measured  activation  energy  and 
apparent  capture  cross  sections  can  be  used  with  the  general  expressions  for  emission 
from  a  single  trap  level  [83], 

A  C(T,  t )  =  AC(7\0)  exp[-  ent\  IV-3 


and 


e„  =  v(T)<jnaN  c  (T)  exp 


TOT 

kT 


IV-4 


in  order  to  model  the  capacitance  transient  due  to  that  level.  In  these  expressions,  A C(T, 
t )  is  the  capacitance  transient  at  time  t  following  the  end  of  the  trap  filling  pulse,  en  the 
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emission  rate,  v(T)  the  rms  thennal  velocity  and  NC(T)  the  conduction  band  density  of 
states  at  temperature  T.  Peak  broadening  may  be  introduced  into  the  above  model  by 
integrating  equation  IV-3  over  a  distribution  of  trap  energies  centered  at  Ejot- 

oo 

A C(t)  =  AC(0)J g(E)  exp[-e„  (E)t\lE ,  IV-5 

0 


where  g{E)  is  the  Gaussian  distribution  function  given  by 


g(E)  =  f-d=-]exp[-(£ -Emiy- /2S=] 
v  V  2  tc  S  J 


This  simple  model  was  applied  to  alloy-broadened  DLTS  data  by  Omling,  et  al.  [86]. 
The  parameter,  S,  is  a  measure  of  broadening.  Care  must  be  taken  when  applying  this 
method  to  the  case  of  overlapping  peaks,  as  the  choice  of  S  for  one  peak  can  alter  the 
positions  and  magnitudes  of  the  adjacent  peaks,  with  the  magnitudes  being  the  most 
susceptible  to  error. 


4.1.2  Trap  Concentration 

It  can  be  shown  [83]  that  when  the  trap  concentration,  Nt,  is  much  less  than  the 
net  shallow  donor  concentration  in  the  space  charge  region,  Noet  (=  No  -  Na), 

1  AC 

Nr=--2 - NnDe‘,  IV-7 

/,  Q 

where  it  is  also  assumed  that  AC  corresponds  to  complete  filling  of  all  the  traps  in  the 
region  of  interest  during  the  forward  voltage  pulse.  The  factor, yj,  accounts  for  the 
thickness  of  the  depletion  region  in  which  the  trap  level  is  already  filled  before 
application  of  the  forward  voltage  pulse.  The  factor^  is  given  as  [87] 
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and  the  thickness  X  is  given  as  [87] 


(w  -A)2  ~(wf  -  Ay 


w 


A  = 


2s(Ef  -Et) 


e 2  -Nr 


1/2 


IV-9 


where  is  the  depletion  depth  under  the  reverse  bias  voltage,  and  vvy  is  the  depletion 
depth  during  the  forward  bias  voltage  pulse  as  depicted  in  figure  IV-2.  Equation  IV-7 
can  be  applied  to  specific  rate  window  plot,  where  instead  of  A C/C0  (=  [C(0)  -  C(oo)]/CG), 
the  peak  magnitude  is  a  lesser  value  of  8C/C0  (=[C(C)  -  C(C)]/CG)  .  In  this  case,  equation 
IV-7  can  be  re-written  in  terms  of  the  ratio,  r  =  tXtu  as  [83] 


„  1  ^  SC 

Nt  = - 2 - 

fx  CQ 


exp{  r/(l  -r)  In  r 

(1  ~r) 


N 


D  • 


IV- 10 


For  r  =  3  and  r  =5,  two  commonly  used  ratios  in  this  dissertation  research,  the  value  in 
brackets  becomes  5.20  and  3.74,  respectively. 


4.1.3  Pulse  Width  Dependence 

During  application  of  the  forward  bias  voltage  pulse,  the  previously  empty  point 
defect  traps  of  a  given  species  capture  electrons  according  to  [28] 

N(t)  =  Vr[l-exp(-c,/)],  IV- 11 

where  N(t)  and  Nt  are,  respectively,  the  filled  and  total  concentrations  of  the  specific  trap 
species,  and  c„  is  the  electron  capture  rate  (introduced  in  equation  11-35)  for  that  trap.  By 
varying  the  filling  pulse  width  tp  and  recording  the  resulting  AC/C,  equations  IV-7  and 
IV- 1 1  can  be  used  to  directly  measure  cn  at  a  given  temperature.  By  inserting  typical 
values  of  n~10  cm'  ,<v>~10  cm/s,  and  an  ~  1 0'  cm“  into  equation  11-35,  it  is  found 
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that  cn~  10  s'  .  In  this  case,  the  traps  would  be  63%  filled  in  10  ns.  Thus,  a  pulse  width 
tp  on  the  order  of  10  ns  is  required  to  measure  trap  filling  versus  tp.  For  the  DLTS  system 
used  in  this  study,  the  lowest  available  tp  is  20  us,  so  such  a  measurement  can  not  be 
performed. 

Dislocation-related  defects  are  known  to  display  capture  kinetics  that  deviate 
significantly  from  the  form  of  equation  IV- 1 1.  It  has  been  observed  in  many  materials 
including  Si  [88],  GaAs  [89],  and  GaN  [90]  that  dislocation-related  traps  N(t)  do  not 
approach  a  saturation  value  except  when  very  long  filling  pulse  widths  (>100  ms)  are 
used.  The  increase  in  N(tp)  (or  A C(tp)/C)  is  linear  when  plotted  versus  log(C).  It  has  been 
observed  that  the  DLTS  signature  (Ea  and  r>na)  for  these  dislocation-related  traps  is  often 
the  same  as  for  known  point  defects  in  the  material.  It  was  also  observed  [88]  that  for 
short  pulse  widths  (when  only  a  small  fraction  of  Nt  are  filled)  normal  trap  filling  takes 
place  according  to  equation  IV- 1 1.  A  simple  model  [91]  was  proposed  to  describe  this 
dislocation-related  capture  behavior  in  terms  of  a  repulsive  Coulomb  potential,  (pd ,  that 
reflects  the  number  of  charges  captured  at  the  dislocations.  Figure  IV-3  shows  an  energy 
band  illustration  of  this  model.  With  the  addition  of  this  repulsive  potential,  the  electron 
capture  rate  valid  for  all  filling  pulse  widths  becomes 

^  =  (NT  ~  N)n  <V>Crn  exp[-  q<Pd  (tp  )/ kT  l  IV-  ■ 12 

where  (pAtp)  increases  with  tp  and  has  the  effect  of  reducing  the  total  number  of  electrons, 

ri{,  available  for  trapping  in  close  proximity  to  the  defect.  As  suggested  in  figure  IV-3, 
the  basic  physical  properties  of  the  defect-related  point  defect  remain  the  same  as  those  of 
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Figure  IV-3.  Energy  band  illustration  of  trapping  at  dislocation-related  point  defects. 
The  situation  for  small  occupation  numbers  (N(t)«Nj)  resembles  that  for  isolated  point 
defects  (after  Omling  et  al.  [88]). 
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the  isolated  point  defect. 


4.1.4  Field  Dependence 

In  all  DLTS  measurements,  carrier  emission  is  measured  under  the  influence  of  an 
electric  field.  If  a  trap  is  a  single  acceptor,  it  starts  out  neutral  and  becomes  negatively 
charged  upon  capture  of  an  electron.  If  a  trap  is  a  single  donor,  it  starts  out  positively 
charged  and  becomes  neutral  upon  capture  of  an  electron.  In  the  latter  case,  emission  of 
the  electron  takes  place  against  the  Coulombic  attraction  to  the  positively  charged  donor 
ion.  Expressed  as  a  function  of  distance  r  from  the  donor  ion,  the  potential  energy  is  U(r) 
=  -q~ler.  Under  the  application  of  the  electric  field  F  in  the  x-direction  the  potential 
energy  may  be  expressed  in  one  dimension  as  [72] 

e2 

U  (x)  = - +  eFx.  IV- 13 

ex 

Differentiating  with  respect  to  x  and  setting  the  expression  equal  to  zero  yields  the 
distance  xmin  at  which  the  Columbic  attraction  peaks.  Plugging  this  value  back  into 
equation  IV- 13  yields  the  new  Coulombic  energy  barrier  that  must  be  overcome,  and  this 
energy  barrier  is  seen  to  be  lower  than  the  zero-field  energy  barrier  by  an  amount  A  U 


given  as 


A  U  =  2e\ 


1/2 


IV- 14 


This  so-called  Poole-Frenkel  barrier  lowering  becomes  apparent  in  DLTS  measurements 
as  a  decreasing  linear  plot  of  Ea  or  ln(e„)  versus  F  .  Whenever  such  field-dependent 
behavior  occurs,  the  trap  may  be  assumed  to  be  donor-like.  The  absence  of  such 
behavior,  however,  is  not  necessarily  confirmation  that  the  trap  is  acceptor-like  [92].  A 
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classic  example  is  the  EL2  deep  donor  in  GaAs,  for  which  it  has  been  shown  that  the 


Poole-Frenkel  effect  is  greatly  suppressed.  Ganichev  et  al.  [92]  presented  a  model  that 
showed  how  the  presence  of  an  Et  can  cause  this  suppression.  It  was  also  shown  that 
classical  Poole-Frenkel  emission  can  only  be  expected  at  low  field  strengths,  whereas 
phonon  assisted  and  direct  tunneling  begin  to  dominate  at  higher  field  strengths  [92]. 
The  latter  two  processes  dominate  at  all  field  strengths  in  the  case  of  neutral  traps. 

The  field-enhanced  emission  processes  are  illustrated  in  figure  IV-4  for:  a)  a  donor-type 
electron  trap;  and  b)  an  acceptor-type  electron  trap. 


Figure  IV-4.  Potential  barrier  and  emission  processes  a)  for  a  charged  (donor-type) 
electron  trap  and  b)  for  a  neutral  (acceptor-type)  electron  trap  (after  Ganichev  [92]). 


4.2  Temperature  Dependent  Hall  Effect  Measurements 

The  corrected  n(T)  data  obtained  using  the  procedures  of  section  3.4.2  may  in 
principle  be  fitted  by  equation  II-6  in  order  to  model  the  given  donor  and  acceptor  levels 
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and  concentrations.  In  this  dissertation  research,  the  fit  was  performed  with  the 
assumptions  that  n(T)  is  detennined  by  two  dominant  donor  levels  located  within  several 
kT  of  the  Fermi  level  and  that  all  other  donor  and  acceptor  energies  are  more  than  several 
kT  below  the  Fenni  level  at  all  temperatures.  The  assumption  of  two  dominant  donors  is 
based  on  previous  results  by  Look  [33]  and  Green  [17]  for  epitaxial  layers  of  GaN  grown 
by  HVPE.  In  those  studies,  n(T)  fitting  indicated  that  a  shallow  donor  level  believed  to 
be  an  impurity  dominated  the  carrier  concentration  in  addition  to  a  deeper  donor  level 
believed  to  be  the  Vn  donor.  The  introduction  of  a  third  donor  with  any  chosen 
concentration  and  energy  level  provides  the  opportunity  to  obtain  a  better  fit  to  the  data  at 
the  expense  of  weakening  the  model  with  too  many  parameters.  The  second  assumption 
that  all  acceptors  are  located  more  than  a  few  kT  below  the  Fermi  level  and  thus  ionized 
for  20  K  <  T  <  300  K  allows  the  summation  of  k  acceptors  in  equation  II-6  to  be 
simplified  as 

'ZN#’*'ZN*=NT‘  IV-15 

k  k 


where  the  parameter  Nef  is  the  sum  of  all  acceptor  levels.  This  second  assumption  is 
based  on  the  findings  of  Look  [33]  and  Green  [17].  Using  equations  II-3,  II-4,  and  IV-15, 
equation  II-6  may  be  expressed  in  the  two  donor  model  as 
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where  Edi  and  Ed?  are  the  ionization  energies  of  the  first  and  second  donors.  The  fit  was 
accomplished  on  a  Mathcad  worksheet  by  first  plotting  the  carrier  concentration  data 
nj  T)  and  the  modeled  carrier  concentration  nm{  T)  from  equation  IV- 16  versus  1000/T  on 
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a  logarithmic  vertical  scale.  The  model  parameters  Ndi,  Tdi,  Nm,  Em,  and  Na  were  then 
individually  adjusted  until  the  model  curve  overlapped  the  data  points  over  the  widest 
temperature  range  possible.  When  the  model  could  not  match  the  data  over  the  entire 
temperature  range,  parameters  were  chosen  in  order  to  fit  to  the  higher  data  points.  Data 
in  higher  temperature  range  (T  >  70  K  )  are  less  susceptible  to  error  that  is  introduced 
from  the  near-interfacial  conductive  layer  subtraction. 
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V.  Experimental  Results  and  Analysis 


The  data  obtained  and  analyzed  according  to  the  procedures  in  chapters  III  and  IV 
are  plotted  and  discussed  below.  The  results  are  grouped  according  to  characterization 
technique,  with  DLTS  results  being  presented,  first;  Hall  effect  results,  second;  and 
cathodoluminescence  results,  third.  Analysis  is  presented  along  with  experimental 
observations,  with  an  additional  analysis  section  being  included  to  discuss  the  DLTS 
trends. 

5.1  Deep  Level  Transient  Spectroscopy  (DLTS) 

The  DLTS  rate  window  spectra  of  the  pre-irradiated  samples  were  found  to  vary 
significantly  with  aluminum  mole  fraction  and  silicon  doping  concentration. 

Experimental  results  from  each  wafer  are  thus  presented  separately,  with  comparisons 
being  discussed  in  section  5.1.9. 

5.1.1  Unirradiated  GaN  (wafer  A2) 

Figure  V-l  shows  the  DLTS  spectrum  of  unirradiated  GaN  (sample  A2-4b) 
recorded  with  the  pulse  parameters  V/=  0  V,  Vr  =  -4  V,  and  tp  =  0.1  s.  Six  separate  peaks 
are  labeled  in  the  figure.  An  additional  peak  (labeled  E)  was  revealed  with  DDLTS,  but 
the  DLTS  signatures  could  not  be  measured  for  a  reason  discussed  below.  Arrhenius 
plots  were  obtained  for  the  peaks  A,  B,  C,  and  D,  and  these  plots  are  displayed  in  the 
inset  of  figure  V-1.  The  measured  DLTS  parameters  corresponding  to  these  peaks  are 
listed  in  table  V-l .  It  is  seen  that  with  the  exception  of  trap  A,  the  measured  DLTS 
signatures  agree  well  with  those  previously  reported  in  the  literature. 
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Figure  V-l.  DLTS  spectrum  of  unirradiated  GaN  (sample  A2-4b).  Trap  labels  are 
assigned  based  on  similarity  to  previously  reported  traps.  The  Arrhenius  plots  of  peaks 
A,  B,  C,  and  D  are  shown  in  the  inset. 
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Table  V-l.  Measured  DLTS  trap  parameters  in  GaN  and  references  to  similar  reports 


Parameter 

E 

D 

c 

Cl 

B 

A 

W  (K) 

-110 

137 

207 

250-270 

301 

365 

Ej ot  (eV) 

0.22±0.02 

0.44±0.04 

- 

0.59±0.02 

0.82±0.04 

crOT  (10-16  cm2) 

2.8 

1300 

- 

56 

1500 

Similar  traps 

and  LTot 
(eV) 

ED2  0.11 

El90  0.19 

D2  0.16 

D  0.25 

E02  0.27 

El  0.26 

D,  0.20 
DLN,93  0.24 

C  0.41 

ES3  0.40 

C,  0.43-0.48 
ES4  0.45±  .1 

B  0.62 

E2  0.58 

D2  0.60 

A  0.67 

E05  0.61 

E394  0.67 

D395  0.67 

3  Peak  temperature  determined  with  47s  rate  window. 


The  peak  that  begins  to  develop  at  the  high  temperature  end  (~  450  K)  of  the 
spectrum  only  appears  for  filling  pulse  widths  larger  than  1  ms.  Limitations  on  the 
temperature  and  rate  window  do  not  allow  the  trap  parameters  to  be  extracted  for  this 
trap,  but  other  people  have  observed  a  trap  in  this  temperature  range  of  400-450K,  and  it 
has  been  labeled  Ai,  [49],  DLN4  [93],  and  E4  [94],  The  reported  activation  energy  of  Ai 
ranges  from  0.86  to  0.91  eV.  The  concentration  of  this  trap  has  been  shown  to  increase 
with  all  kinds  of  radiation  [8],  and  this  increase  shows  a  correlation  with  an  increase  in 
the  density  of  trap  E,  which  is  believed  to  be  the  isolated  nitrogen  vacancy.  In 
unirradiated  GaN,  however,  the  intensity  of  Ai  has  been  reported  to  be  anti-correlated 
with  trap  E.  As  pointed  out  by  Look  [8],  both  of  these  observations  are  consistent  with 
the  identification  of  this  trap  as  being  nitrogen  interstitial  (Ni)  related. 

The  next  peak,  labeled  A,  is  seen  at  365  K  using  a  47  Hz  rate  window.  Other 
people  have  observed  a  trap  in  this  temperature  range  which  has  been  labeled  A  [49], 
E05  [35],  E3  [94],  and  D3  [95].  The  Arrhenius  fit  indicates  an  activation  energy  E-tot  = 
0.82  ±  0.04  eV  and  cross  section  ona  =  1.5  x  10'  "  cm-.  This  activation  energy  and  cross 
section  are  both  significantly  greater  than  what  has  usually  been  reported  for  trap  A  [49] 
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(0.67  eV  and  ~10'  cm'),  but  as  will  be  seen  below,  the  peak  seen  here  appears  to 
actually  be  a  superposition  of  peaks  with  different  activation  energies.  When  this  is  the 
case,  the  Arrhenius  plot  from  the  combined  peak  sometimes  indicates  large  activation 
energy  and  cross  section.  The  next  peak,  labeled  B,  is  seen  at  301  K  using  a  47  Hz  rate 
window.  A  trap  has  been  reported  in  this  temperature  range  by  others,  and  it  has  been 
labeled  B  [49],  E2  [94],  and  D2  [95].  The  Arrhenius  fit  indicates  an  activation  energy 
ETOt  =  0.59  ±  0.02  eV  and  cross  section  atm  =  5.6  x  10'  cm  .  The  activation  energy  and 
cross  section  agree  with  previously  reported  values  for  this  trap  (0.58-0.62  eV  and  ~10'14 
cm2). 

The  next  slightly  perceptible  shoulder  peak  on  the  low  temperature  side  of  trap  B 
is  seen  in  the  250-270  K  temperature  range,  and  is  tentatively  assigned  to  another  trap 
labeled  Ci.  Other  people  observed  a  trap  in  the  same  temperature,  and  it  has  been  labeled 
Ci  [49]  and  ES4  [96],  Fang  et  al.  [49]  reported  that  Ci  was  peculiar  to  reactive  molecular 
beam  epitaxy  (RMBE)  GaN  layers.  They  also  reported  that  Ci  was  concentrated  close  to 
the  surface,  and  that  the  DLTS  peak  position  was  sensitive  to  electric  field  strength. 

Auret  et  al.  [96]  reported  that  ES4  is  created  in  the  process  of  sputtered  metal  deposition, 
and  that  it  was  concentrated  within  15  pm  of  the  surface. 

The  next  peak,  labeled  C,  is  seen  at  207  K  using  a  47  Hz  rate  window.  A  similar 
trap  reported  in  this  temperature  range  has  been  labeled  C  [49]  and  ES3  [96].  The 
Arrhenius  fit  indicates  an  activation  energy  Etot  =  0.44  ±  0.04  eV  and  cross  section  ana  = 
1.3  x  10'  ~  cm'.  The  activation  energy  and  cross  section  agree  with  the  two  previously 
reported  values  for  this  trap  (0.40-0.41  eV  and  10'14-10'13cm2).  Fang  et  al.  [49]  have 
reported  that  trap  C  is  associated  with  ion-beam  etching,  and  is  concentrated  within  27 
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jam  of  the  surface.  They  have  also  observed  this  trap  in  GaN  grown  by  plasma-assisted 
MBE.  Auret  et  al.  [96]  reported  that  as  with  ES4  (Ci),  ES3  was  created  during  sputtered 
metal  deposition,  and  that  it  was  concentrated  within  15  pm  of  the  surface. 

An  instance  in  which  trap  C  dominated  the  DLTS  spectrum  is  shown  in  figure  V- 
2.  Here,  it  is  seen  that  the  intensity  of  peak  C  increases  considerably  with  pulse  width. 
The  inset  shows  that  this  peak  height  increase  was  only  apparent  with  the  larger  rate 
windows.  As  shown  in  figure  V-3  the  increase  in  peak  C  was  not  observed  in  a 
subsequent  scan,  indicating  that  that  this  effect  can  be  suppressed  with  annealing  at  450 
K  under  an  applied  reverse  bias  voltage. 

The  next  peak  in  the  DLTS  spectrum  of  figure  V-l,  labeled  D,  is  seen  at  137  K 
using  a  47  s'1  rate  window.  A  similar  trap  reported  in  this  temperature  range  has  been 
labeled  D  [49],  E02  [35],  El  [94],  and  DLN1  [93].  The  Arrhenius  fit  indicates  an 
activation  energy  EYot  =  0.22  ±  0.02  eV  and  cross  section  a,w  =  2.8  x  10'16  cm2.  The 
activation  energy  and  cross  section  are  close  to  previously  reported  values  for  this  trap 
(0.24-0.27  eV  and  10'15-10'14cm2),  however,  DDLTS  measurements  reveal  that  this  peak 
has  a  strong  electric  field  dependence. 

Figure  V-4  shows  the  DDLTS  spectra  (16.5  s'1  rate  window)  of  the  same  GaN 
diode  discussed  above.  The  four  pulsing  conditions  isolate  four  sections  of  the  Vr  =  0  to 
-4  V  depletion  region.  It  is  seen  that  the  summation  of  these  four  different  signals  (solid 
dark  gray  line)  closely  matches  the  signal  obtained  with  the  pulsing  conditions  Vf  =  0,  Vr 
=  -4  V  (dashed  light  gray  line).  The  figure  suggests  that  peaks  A  and  D  are  strongly 
affected  by  electric  field,  and  thus  by  depth  in  the  depletion  region.  Figure  V-5  shows  the 
detailed  components  of  peak  D  under  the  indicated  pulsing  conditions.  Arrhenius  plots 
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Figure  V-2.  DLTS  spectrum  of  unirradiated  GaN  (sample  A2-4b)  at  different  filling 
pulse  widths.  The  trap  C  dominates  at  large  pulse  width  in  this  temperature  scan,  but  it 
greatly  diminished  in  a  subsequent  temperature  scan.  The  inset  shows  various  rate 
windows  for  tp  =10"2  s,  illustrating  that  trap  C  is  only  dominant  for  large  rate  windows. 
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Figure  V-4.  Double  correlated  DLTS  spectra  of  unirradiated  GaN  (sample  A2-4b).  The 
four  pulsing  conditions  isolate  four  sections  of  the  Vr=  0  to  -4  V  depletion  region.  The 
summation  of  these  four  different  signals  (solid  dark  gray  line)  closely  matches  the  signal 
obtained  with  the  pulsing  conditions  V/=  0  ,  V,  =  -4V  (dashed  light  gray  line). 
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corresponding  to  these  peaks  are  shown  in  figure  V-6.  The  peak  that  appears  under  the 
lowest  field  condition  has  activation  energy  £tot  =  0.23  eV  and  cross  section  ona  =  3.9  x 
10'16  cm2,  whereas  the  peak  that  appears  under  the  highest  field  condition  has  activation 
energy  Ejot:  =  0. 13  eV  and  cross  section  crna  =  1.1  x  10'  cm'.  The  two  low-field  peaks 
are  very  similar,  and  their  superposition  dominates  the  emission  properties  of  the 
combined  peak  D.  The  two  high-field  peaks  have  the  same  cross  section  but  different 
activation  energies  (0.13  and  0.19  eV).  This  fact  could  indicate  that  two  different 
electron  traps  are  being  revealed.  In  fact,  shoulders  on  the  low  temperature  side  of  trap  D 
have  often  been  reported  [64].  In  studies  on  GaN  grown  by  MOCVD  and  HVPE,  this  has 
been  labeled  as  trap  E|  [90],  and  has  similarities  to  trap  E  (combination  of  EDI  and  ED2 
[31])  introduced  by  electron  irradiation.  In  GaN  grown  by  RMBE,  deconvolution  of  trap 
D  by  fitting  has  revealed  two  traps,  Di  ( ETot=  0.25  eV  ;  on  =  1.4xl0'14  cm2)  and  D2  {Etot 
=  0. 16  eV  ;  an  =  2.1x10'  cm“),  where  Di  dominates  the  behavior  of  trap  D,  and  D2 
appears  as  a  shoulder  on  the  low  temperature  side  of  Di.  In  this  case,  however,  the  very 
low  capture  cross  section  may  indicate  that  phonon  assisted  tunneling  is  dominating,  and 
what  is  seen  in  figure  V-5  is  simply  trap  D  with  a  strongly  field-dependent  emission  rate. 

Figure  V-7  shows  the  DDLTS  spectra  (16.5  s'1  rate  window)  of  trap  A.  In  this 
case,  the  higher  rate  windows  were  too  noisy  for  Arrhenius  fitting  to  be  accomplished. 
Figure  V-8  shows  the  peak  heights  of  traps  A,  B,  C,  and  D  versus  filling  pulse  width.  It 
is  seen  that  traps  A,  B,  and  D  exhibit  a  linear  increase  with  log(tp)  up  through  tp  =  10  ms. 
Trap  C  increases  in  this  manner  without  saturation  up  through  tp  =  500  ms.  As  discussed 
in  section  IV,  this  behavior  indicates  that  all  four  of  these  defects  are  concentrated  in  the 
vicinity  of  dislocations. 
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Figure  V-6.  Arrhenius  plots  corresponding  to  the  low  temperature  peaks  shown  in  figure 
V-5.  Inset  are  the  activation  energies  plotted  against  the  square  root  of  the  mean  field 
strength.  The  error  bars  represent  the  standard  errors  detennined  from  data  variance. 
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Figure  V-7.  Double  correlated  DLTS  spectrum  of  peak  A  in  unirradiated  GaN  (sample 
A2-4b). 
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Figure  V-8.  Peak  heights  of  four  traps  versus  filling  pulse  width  in  GaN  (sample  A2 
4b).  The  47  s'1  rate  window  was  used  for  this  comparison. 


5.1.2  Irradiated  GaN  (Wafer  A2) 

Figure  V-9  shows  the  DLTS  spectrum  of  sample  A2-4a  irradiated  at  two 
consecutive  1.0  MeV  electron  doses.  It  is  apparent  that  the  peaks  A,  B,  and  C  do  not 
increase  with  radiation  damage.  A  small  decrease  in  some  of  these  may  be  attributed  to 
an  increased  in  the  Schottky  barrier  diode  series  resistance.  Only  one  radiation-induced 
peak,  labeled  E,  is  shown.  Arrhenius  fitting  could  not  be  accomplished  for  this  peak 
because  of  carrier  freeze-out.  In  this  situation,  radiation-induced  deep  acceptor  levels 
cause  the  carrier  concentration  to  drop  sharply  as  temperature  is  reduced.  This 
phenomenon  is  especially  severe  at  the  second  dose  as  shown  in  figure  V-9.  Figure  V-10 
shows  that  the  diode  capacitance  drops  off  considerably  at  low  temperatures  due  to 
carrier  freeze-out.  The  C-V  curves  taken  after  the  l.lxlO17  cm"2  dose  are  shown  at 
various  temperatures  in  figure  V-l  1.  In  this  figure,  it  is  seen  that  a  normal  C-V 
relationship,  showing  a  decrease  in  capacitance  with  increase  in  Vr,  is  not  obtained  for  T 
<  160  K.  A  comparison  with  a  nonnal  C-V  curve  at  130  K  is  shown  in  figure  V-l 2  for 
samples  before  and  after  electron  irradiation.  Figure  V-l 3  shows  double-correlated 
DLTS  spectra  of  GaN  (sample  A2-5a)  irradiated  with  2.5xl016  cm'2  of  0.62  MeV 
electrons.  The  rate  window  used  in  this  figure  is  610  s'1,  so  the  peaks  are  shifted  upward 
in  temperature  relative  to  the  121  s'1  rate  window  of  figure  V-9.  In  figure  V-l  3,  it  is  seen 
that  trap  E  is  rather  broad,  and  almost  certainly  is  a  combination  of  traps  ER1  and  ER2 
reported  in  [35].  Figure  V-14  shows  the  field  dependence  of  the  E  and  D  peaks,  where 
both  DDLTS  spectra  are  generated  in  the  same  spatial  region  defined  by  the  forward 
pulse  pairs  Vp  =  0  V  and  Vp=  -2  V.  This  figure  demonstrates  that  the  variations  in  the 
DDLTS  spectra  of  figure  V-l 3  are  caused  by  the  field  variation  in  the 
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Figure  V-9.  DLTS  spectra  of  GaN  (sample  A2-4a)  irradiated  at  two  consecutive  1.0 
MeV  electron  doses.  Unirradiated  GaN  (sample  A2-4b)  is  overlaid  for  comparison. 
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Figure  V-10.  Capacitance  versus  temperature  for  GaN  (sample  A2-4a)  before  and  after 
two  consecutive  1.0  MeV  electron  irradiations. 
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Figure  V-ll.  C-V  curves  at  various  temperatures  for  1.0  MeV  electron  irradiated  GaN 
(sample  A2-4a). 
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Figure  V-12.  C-V  curve  of  1.0  MeV  electron  irradiated  GaN  (sample  A2-4a),  showing 
the  effect  of  carrier  freeze-out  at  this  temperature.  The  normal  C-V  curve  from 
unirradiated  GaN  (sample  A2-4b)  is  overlaid  for  comparison. 
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Figure  V-14.  Double  correlated  DLTS  spectra  of  peaks  D  and  E  in  irradiated  GaN 
(sample  A2-5a)  showing  the  effect  of  field  strength  in  the  spatial  region  defined  by 
forward  pulse  pair  Vn  =  0  V  and  ffi  =  - 1  V.  The  rate  window  is  610  s'1. 
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DLTS  Signal:  AC/C  (arb.  scale) 
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Figure  V-15.  DLTS  spectra  of  peaks  D  and  E  in  irradiated  GaN  (sample  A2-5a)  before 
and  after  annealing  at  350  °C  for  15  minutes.  The  rate  window  is  610  s'1. 


depletion  region,  and  not  by  an  inhomogeneous  spatial  distribution  of  defects.  Figure  V- 
15  shows  that  trap  E  diminishes  greatly  after  annealing  at  350  °C.  The  extent  of  the 
reduction  caused  by  the  annealing  is  obscured  by  carrier  freeze-out,  however,  it  appears 
that  trap  E  anneals  out  considerably,  and  that  the  remaining  low  temperature  shoulder  of 
peak  D  may  be  similar  to  what  was  seen  in  the  unirradiated  GaN  samples. 

5.1.3  Unirradiated  Alo.i4Ga0.86N  (Wafer  B2) 

Figure  V-16  shows  a  DLTS  spectrum  of  unirradiated  Alo.14Gao.86N  typical  of 
samples  from  wafer  B2.  A  47  s'1  rate  window  is  used,  and  the  pulsing  parameters  are  V/= 
3  V,  Vr  =  -3  V,  and  tp  =  20  ps.  Two  electron  traps,  labeled  PI  and  P2  are  apparent  in  the 
spectrum. 

The  peak  of  figure  V-16  labeled  P2  is  noticeably  broader  than  the  peak  that  would 
be  expected  from  a  discrete  trap  level.  This  kind  of  broad  peak  can  result  from  the 
superposition  of  a  few  traps  with  different  closely  spaced  energy  levels.  With  the 
measurement  parameters  shown  in  figure  V-16,  the  central  value  of  activation  energy  for 
P2  is  Etot  =  0.98  ±  0. 10  eV  and  the  apparent  capture  cross  section  is  ona  =  4x10'  “  cm'. 
Using  the  0.1s  filling  pulse  width,  a  much  cleaner  Arrhenius  plot  is  obtained.  In  that 
case,  the  central  value  of  activation  energy  for  P2  is  Etot  =  0-88  ±  0.02  eV  and  the 
apparent  capture  cross  section  is  a,w  =  5±3  xlO'12  cm2. 

The  most  prominent  electron  trap  level,  PI,  peaks  at  255  K  in  the  47  s'1  rate 
window.  The  Arrhenius  fit,  shown  in  the  inset,  indicates  an  activation  energy  of  ETot  = 
0.54  ±  0.02  eV  and  cross  section  of  cr,ia=  3.9  ±  1.0  x  10'14  cm2.  This  level  was  also 
observed  by  Legodi  et  al.  [42]  in  Alo.12Gao.88N  grown  by  the  same  plasma-assisted  MBE 
method.  They  suggested  that  this  electron  trap  is  similar  to  the  one  observed  in  MOCVD 
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Figure  V-16.  DLTS  spectrum  of  pre-irradiated  Alo.14Gao.86N  (sample  B2-21d).  The 
Arrhenius  plots  of  the  two  peaks  are  shown  in  the  inset. 
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-grown  Alo.12Gao.88N  by  Gotz  et  al.  [97].  Gotz  et  al.  labeled  this  electron  trap  DLN2  and 
reported  an  activation  energy  of  Etot=  0.61  ±  0.02  eV.  The  authors  also  proposed  that 
this  DLN2  trap  in  Alo.14Gao.86N  corresponded  to  a  trap  in  GaN  having  ETot=  0.5  eV. 

They  also  reported  that  the  activation  energy  of  DLN2  decreased  linearly  with  the 
increasing  square  root  of  field  strength,  indicating  that  the  trap  was  donor-like.  In  order 
to  determine  if  PI  shares  this  donor-like  characteristic,  the  electric  field  dependence  of 
ETOt  is  measured  by  performing  DDLTS  measurements.  Figure  V-17  shows  Eror  versus 
the  square  root  of  the  average  electric  field  strength  in  a  narrow  spatial  segment  of  the 
depletion  region  defined  by  the  forward  pulse  pairs  Va  =  -1  V  and  Vn  =  -2  V.  The 
reverse  bias  was  increased  from  -2  to  -6  V  to  obtain  an  average  electric  field  strength 
ranging  from  1.3xl03  V/cm  to  9.6x1 03  V/cm.  It  is  apparent  that  there  is  no  clear 
decreasing  trend  in  ETot  as  the  average  electric  field  strength  is  increased  in  this  range. 
The  inset  of  figure  V-17  shows  the  47  s'1  rate  window  of  PI  at  electric  field  strengths 
1 ,3x  1 0  and  9.1x10  V/cm,  where  it  is  clear  that  the  increase  in  electric  field  strength 
causes  no  significant  shift  in  the  peak  position.  This  is  quite  different  than  the  situation 
reported  for  DLN2  by  Gotz  et  al.,  however,  in  that  work,  the  field  strengths  were  an  order 
of  magnitude  greater  due  to  a  higher  carrier  concentration  n  ~  5x10  cm"  .  Figure  V-18 
shows  the  increase  in  the  PI  and  P2  peak  heights  with  pulse  width.  Here,  it  is  seen  that  in 
the  pulse  width  range  of  20  ps  to  1  ms,  the  peak  heights  increase  linearly  with  the 
logarithm  of  pulse  width.  Again,  this  may  indicate  that  the  defects  corresponding  to  these 
peaks  are  associated  with  line  defects,  as  discussed  in  chapter  IV.  Because  increasing 
occupation  of  the  electron  traps  gives  rise  to  a  repulsive  potential,  it  appears  that  these 
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Figure  V-17.  Electric  field  dependence  of  electron  trap  PI  activation  energy  in 
Alo.14Gao.86N  (sample  B2-8b).  DDLTS  signals  are  measured  in  the  spatial  region  defined 
by  the  forward  voltage  pulse  pair  Vn  =  -1  and  Vn  =  -2  V,  with  field  increase  provided  by 
increasing  reverse  bias  Vr.  Error  bars  apply  to  random  sources  of  error  in  these 
comparative  measurements.  Inset  shows  47s'1  peak  positions  when  field  strength  is 
raised  nearly  one  order  of  magnitude.  No  peak  temperature  shift  is  observed. 
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DLTS  Signal:  AC/C 


Figure  V-18.  DLTS  spectrum  of  unirradiated  Alo.14Gao.86N  (sample  B2-21d)  at  various 
pulse  widths.  The  peak  heights  of  PI  and  P2  are  plotted  versus  filling  pulse  width  in  the 
inset. 
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Figure  V-19.  Normalized  peak  height  of  electron  trap  PI  versus  reverse  bias  voltage. 
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traps  both  behave  as  acceptors.  The  reverse  bias  voltage  dependence  of  PI  peak  height  is 
shown  in  figure  V-19.  The  resulting  plot  is  characteristic  of  a  trap  distributed  throughout 
the  thickness  of  the  epitaxial  layer  rather  than  being  concentrated  near  the  surface. 

5.1.4  Irradiated  Alo.i4Ga0.86N  (Wafer  B2) 

Figure  V-20  shows  the  DLTS  spectrum  of  the  irradiated  Alo.14Gao.86N  sample  B2- 
21b  measured  after  1.0  MeV  electron  irradiation  to  a  fluence  of  9x10  cm'  and  about  3 
hours  of  room  temperature  annealing.  Measurements  taken  with  two  different  filling 
pulse  widths  are  presented.  It  is  apparent  from  the  spectrum  corresponding  to  the  0.1s 
pulse  width  that  four  distinct  electron  traps  indicated  by  Rl,  R2,  R3,  and  R4  were 
generated  from  the  radiation  damage.  Arrhenius  analysis  was  most  readily  applied  to 
peaks  in  the  20  ps  pulse  width  spectrum.  The  apparent  activation  energy  and  the 
apparent  cross  section  for  the  electron  traps  of  R2,  R3,  and  R4,  determined  from  the  inset 
of  figure  V-21  are  0.21  eV  and  7xl0~17  cm2,  0.26  eV  and  1.4xl0'16  cm2,  and  0.33  eV  and 
1.3x10'  cm",  respectively.  An  Arrhenius  plot  for  Rl  could  not  be  made  due  to  lack  of 
resolution.  In  order  to  deconvolve  the  peaks  in  the  20  ps  pulse  width  spectrum,  the 
measured  parameters  from  the  Arrhenius  plot  were  used  with  equations  IV-4  to  IV-6  to 
generate  the  radiation  damage-induced  peaks  Rl,  R2,  R3,  and  R4  seen  as  the  dotted  lines 
in  figure  V-21.  In  the  application  of  these  equations  to  each  trap,  the  apparent  cross 
section,  ona,  was  allowed  to  vary,  the  activation  energy,  ETot,  was  held  fixed  at  the 
measured  value,  and  the  broadening  parameter,  S,  and  DLTS  signal  magnitude,  AC,  were 
chosen.  The  ona  that  produced  the  best  fit  was  within  the  error  bounds  of  the  measured 
ana  for  each  trap.  In  this  work,  the  S  corresponding  to  trap  PI  was  determined  from  the 
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DLTS  Signal:  AC/C 


Figure  V-20.  DLTS  spectrum  of  1.0  MeV  electron  irradiated  Alo.14Gao.86N  (sample  B2- 
21b)  at  short  (20  ps)  and  long  (0.1  s)  pulse  width. 
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Figure  V-21.  DLTS  spectrum  of  1.0  MeV  electron  irradiated  Alo.14Gao.86N  (sample  B2- 
21b)  at  20  ps  pulse  width.  Inset  is  the  Arrhenius  plot  constructed  from  visible  peak 
positions.  Open  circles  are  measured  data.  Dotted  lines  are  the  individual  peaks  modeled 
using  equations  IV-4  to  IV-6.  Dashed  line  is  spectrum  from  the  pre-irradiated  data.  The 
solid  line  is  the  sum  of  all  the  peaks. 
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pre-irradiated  signal.  It  was  then  found  that  this  value  had  to  be  scaled  down  for  the 
lower  energy  peaks  to  produce  a  good  fit  to  the  data. 

Traps  Rl,  R2,  and  R3  have  not  been  previously  reported,  as  such,  in  AlGaN.  One 
electron  trap,  EpRl  (0.187  eV),  was  previously  reported  in  proton-irradiated 
Alo.12Gao.14N  using  thermally  stimulated  capacitance  (TSCAP)  [42],  Though  the  authors 
drew  no  correlation  with  known  defects  in  GaN,  it  appears  likely  to  be  the  same  as  our 
observed  trap  R2,  or  perhaps  a  superposition  of  Rl  and  R2.  Traps  Rl,  R2,  and  R3  are 
similar  to  ER1,  ER2,  and  ER3  observed  in  GaN  [35].  As  was  discussed  in  chapter  II,  it 
has  been  reported  that  ER1  corresponds  to  the  isolated  nitrogen  vacancy,  and  that  ER2  is 
also  related  to  the  nitrogen  vacancy  with  some  difference  in  microscopic  configuration, 
causing  it  to  have  different  cross  sections  and  activation  energies  [31].  It  is  noteworthy 
that  R2  and  R3  are  50  and  60  meV  deeper  than  the  corresponding  traps  in  GaN.  Equation 
III-2  indicates  that  the  increase  in  electron  effective  mass  would  cause  a  shallow 
(hydrogenic)  donor  to  be  only  7  meV  deeper.  On  the  other  hand,  the  band  gap  energy 
difference  between  Alo.14Gao.g6N  and  GaN  is  approximately  270  meV,  so  it  is  clear  that 
the  energy  of  these  two  traps  is  not  fixed  with  respect  to  the  vacuum  level  but  is,  in  fact, 
partially  related  to  the  position  of  the  conduction  band  minimum. 

The  electron  trap  parameters  measured  in  unirradiated  and  irradiated  Alo.14Gao.86N 
are  summarized  in  table  V-2.  Although  an  Arrhenius  plot  was  not  obtained  for  trap  Rl,  a 
range  of  possible  values  for  its  activation  energy  can  be  proposed  by  assuming  a  capture 
cross  section  based  on  GaN  reports  in  the  literature.  Two  independent  reports  showed 
ER1  to  have  a  cross  section  5%  that  of  ER2  [35,  31].  In  another  report,  the  cross  section 
of  ER1  is  half  that  of  ER2  [37].  Using  these  values  as  lower  and  upper  bounds  of  the 
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Table  V-2.  Measured  defect  parameters  in  1.0-MeV  electron  irradiated  Alo.14Gao.86N  and 


references  to  simi 

lar  reports. 

Parameter 

R1 

R2 

R3 

R4 

PI 

P2 

Ge ak2  (K) 

~  110 

138 

166 

208 

253 

381 

Etot  (eV) 

0.21  ±0.02 

0.26  ±  0.02 

0.33  ±0.03 

0.54  ±0.02 

0.98±  0.10 

oM  10" 16  cm2) 

0.7 

1.0b 

1.4 

1.25b 

1.3 

1.25b 

400 

40000 

Similar  Traps  and 

Ej ox  (eV)  in  GaN 

ERl  0.13 

EDI  0.06b 

T1  0.12 
(0.08b) 

ER2  0.16 

ED2  0.1  lb 
T1A  0.16 

ER335  0.20 

Similar  Traps  and 
fiTOT(eV)  in 
AlxGai.xN 

EpRl  0. 1 9C 
x  =  0.12 

EOA2 

0.57 
x  =  0.12 

aPeak  temperature  determined  with  47s"1  rate  window. 
b  Fitted  value 

L  Determined  by  TSCAP.  Authors  in  Ref.  42  did  not  make  attribution. 


cross  section  of  R1  relative  to  R2,  activation  energy  of  0.15  ±  0.02  eV  is  obtained  from 
the  fitted  R1  peak  in  figure  V-21.  Trap  R4  does  not  have  an  obvious  analog  in  irradiated 
GaN,  suggesting  that  it  is  possibly  related  to  aluminum  displacement.  Assuming  binding 
energies  comparable  to  GaN,  the  1 .0  MeV  electron  energy  used  here  is  above  the 
threshold  energy  values  needed  for  displacing  the  three  constituents  of  Alo.14Gao.86N. 

Figure  V-22  shows  the  DLTS  spectrum  before  and  after  annealing  at  430  and  450 
K  sequentially.  After  cycling  the  diode  up  to  430  K,  the  peak  heights  of  both  traps  R2 
and  R3  decreased  significantly.  Some  further  reduction  was  observed  after  cycling  up  to 
450  K.  Trap  R4,  however,  was  thermally  stable  up  to  450  K.  Figure  V-23  shows  the 
DLTS  spectrum  of  Alo.14Gao.86N  after  0.62  MeV  electron  irradiation  with  a  fluence  of 
2.5x10  cm"  .  In  this  sample,  the  carrier  concentration  was  lower,  and  carrier  freeze-out 
was  a  problem  at  low  temperatures.  Nevertheless,  trap  R4  is  clearly  induced  by  the 
radiation,  and  the  introduction  rate  of  this  trap  with  0.62  MeV  electrons  is  comparable  to 
the  introduction  rate  with  1 .0  MeV  electrons.  In  this  figure,  it  is  also  clear  that  the 
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Figure  V-22.  DLTS  spectrum  of  1.0  MeV  electron  irradiated  Alo.14Gao.86N  (sample  B2- 
21b)  before  and  after  annealing.  The  dashed  and  dotted  lines  represent  the  spectrum  after 
cycling  up  to  430  and  450  K,  respectively. 
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Figure  V-23.  DLTS  spectrum  of  0.62  MeV  electron  irradiated  Alo.14Gao.86N  (sample 
B2-21d)  before  and  after  350  °C  annealing.  The  spectrum  of  unirradiated  Alo.14Gao.86N 
(sample  B2-21d)  is  overlaid  (dotted  line). 
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concentration  of  trap  R4,  which  was  seen  in  figure  V-22  to  be  stable  at  177  °C  (450  K), 
decreases  considerably  following  annealing  at  350  °C  (623  K)  for  15  minutes. 

5.1.5  Unirradiated  Al0.2o Ga08oN  (Wafer  Cl) 

Figure  V-24  shows  a  DLTS  spectrum  of  unirradiated  Alo.2oGa0.8oN  typical  of 
samples  from  wafer  Cl.  A  46  s'1  rate  window  is  shown,  and  the  pulsing  parameters  are 
Vf=  1  V,  Vr  =  -3  V,  and  tp  =  1  ms.  Three  electron  traps,  labeled  P0,  PI,  and  P2  are 
apparent  in  the  spectrum.  The  DLTS  signals  of  all  three  traps  are  seen  to  be  relatively 
small  because  of  the  larger  carrier  concentration  (~  3xl017  cm'3  at  295  K)  in  samples 
from  this  wafer. 

The  peak  that  appears  at  158  K  in  the  46  s'1  rate  window  of  figure  V-24,  labeled 
P0,  may  correspond  to  trap  D  in  GaN.  However,  resolution  of  the  DLTS  signal  is  too 
poor  for  an  Arrhenius  plot  to  be  constructed,  so  this  assignment  is  tentative.  The  peaks 
labeled  PI  and  P2  correspond  to  the  peaks  given  the  same  designation  in  Alo.14Gao.86N. 
The  Arrhenius  plots  corresponding  to  these  two  electron  traps,  shown  in  the  inset  of 
figure  V-24,  indicate  activation  energies  of  L’tot  =  0.67  ±  0.03  and  0.92  ±  0.06  eV,  and 
cross  sections  of  ona=  1.6xl0'13  and  9.0xl0'14  cm2  for  PI  and  P2,  respectively. 
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DLTS  Signal:  AC/C 


Figure  V-24.  DLTS  spectrum  of  unirradiated  Al0.20Ga0.80N  (sample  Cl-38c).  The 
Arrhenius  plots  of  the  two  peaks  are  shown  in  the  inset. 
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5.1.6  Irradiated  Al0.2oGao.8oN  (Wafer  Cl) 


Figure  V-25  shows  the  DLTS  spectrum  of  irradiated  Alo.20Gao.80N  sample  Cl -3  8c 
measured  after  1.0  MeV  electron  irradiation  to  a  fluence  of  l.lxlO17  cm'2  and  about  3 
hours  of  room  temperature  annealing.  The  rate  window  is  46  s'1,  and  the  pulsing 
parameters  are  V/=  1  V,  Vr  =  -3V,  and  tp  =  1  ms.  The  apparent  activation  energy  and 
cross  section  of  electron  trap  R4,  detennined  from  the  Arrhenius  plot  shown  in  the  inset 
of  figure  V-25,  are  E-mj  =  0.38  ±  0.02  eV  and  ona  =  2.7xl016  cm2.  Arrhenius  plots  for 
Rl,  R2,  and  R3  could  not  be  made  due  to  lack  of  resolution.  The  electron  trap  parameters 
measured  in  unirradiated  and  irradiated  Alo.20Gao.80N  are  summarized  in  table  V-3.  Also 
shown  in  figure  V-25  is  the  DLTS  spectrum  of  the  same  Alo.20Gao.80N  sample  following 
annealing  at  350  °C  for  15  minutes.  Figure  V-26  shows  a  similarly  irradiated 
Alo.20Gao.80N  sample,  and  in  this  case  the  sample  was  annealed  at  400  °C  for  15  minutes. 
Comparison  of  the  two  figures  indicates  that  the  radiation-induced  traps  anneal  out  nearly 
entirely  in  the  350  to  400  °C  temperature  range. 


Table  V-3.  Measured  defect  parameters  in  1.0-MeV  electron  irradiated  Alo.20Gao.80N  and 
references  to  similar  reports. _ 


Parameter 

R2 

R3 

P0 

R4 

PI 

P2 

TBe ak2  (K) 

-177 

-208 

158 

245 

296 

405 

Ftot  (eV) 

0.38±0.02 

0.67  ±0.03 

0.92±  0.06 

crOT  (1016  cm2) 

2.7 

1600 

900 

Similar  Traps  and 
Ftot  (eV)  in  GaN 

ER2  0.16 
ED2  0.1  lb 
T1A  0.16 

ER335 

0.20 

D  0.25 

E02  0.27 

El  0.26 
DLN,  0.24 

Similar  Traps  and 
FT0T(eV)  in 
AlxGai.xN 

EpRl  0. 1 9C 
x  =  0.12 

PI  0.16 
x=0.09 

EOA2 

0.57 
x  =  0.12 

aPeak  temperature  determined  with  47s'1  rate  window. 
b  Fitted  value 

c  Determined  by  TSCAP.  Authors  in  Ref.  42  did  not  make  attribution. 
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Figure  V-25.  DLTS  spectrum  of  1.0  MeV  electron  irradiated  Al0.2oGa0.8oN  (sample  Cl- 
38a)  before  and  after  350  °C  annealing.  The  spectrum  of  unirradiated  Alo.20Gao.80N 
(sample  Cl-38c)  is  overlaid  (dotted  line). 
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Figure  V-26.  DLTS  spectrum  of  1.0  MeV  electron  irradiated  Alo.20Gao.80N  (sample  Cl- 
40b)  before  and  after  400  °C  annealing. 
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5.1.7  Unirradiated  Al0.3oGa0.7oN  (Wafer  D1) 

Figure  V-27  shows  a  DLTS  spectrum  typical  of  Alo.30Gao.70N  samples  from  wafer 
Dl.  A  47  s'1  rate  window  is  shown,  and  the  pulsing  parameters  are  Vf=  3  V,  Vr  =  -5  V, 
and  tp  =  20  ps.  One  electron  trap,  labeled  P0,  dominates  the  low  temperature  end  of  the 
spectrum.  The  apparent  activation  energy  and  cross  section  of  electron  trap  P0, 
determined  from  the  Arrhenius  plot  in  the  inset  of  figure  V-27  are  Tstot  =  0-39  ±  0.02  eV 
and  ona  =  2.5x10'  cm-.  The  very  large  capture  cross  section  of  trap  P0  is  responsible  for 
this  trap  having  a  47  s'1  rate  window  peak  at  such  a  low  temperature  as  152  K  when  the 
activation  energy  is  so  large.  Under  the  pulsing  parameters  Vf=  4.5  V,  Vr  =  -1  V,  and  tp  = 
1  ms,  two  peaks  labeled  P 1  and  P2  are  better  resolved  in  the  high  temperature  range  of 
the  47  s'1  rate  window  spectrum  as  shown  in  the  figure.  The  apparent  activation  energy 
and  cross  section  of  electron  trap  PI,  detennined  from  the  Arrhenius  plot  in  the  inset  of 
figure  V-27,  are  ETOt  =  0.94  ±  0.06  eV  and  ona  =  3.0xl0'12  cm2.  Figure  V-28  shows  the 
DLTS  spectrum  with  reverse  bias  voltage  ranging  from  Vr  =  0  to  -5  V.  The  inset  of  this 
figure  shows  that  the  shape  and  location  of  peak  P0  does  not  change  appreciably  with 
reverse  bias  voltage.  Figure  V-29  shows  the  nonnalized  peak  P0  peak  height  versus 
reverse  bias  voltage.  The  shape  of  this  plot  is  characteristic  of  a  trap  distributed 
throughout  the  depth  of  the  epitaxial  layer.  The  electron  trap  parameters  measured  in 
unirradiated  Al0.30Ga0.70N  are  summarized  in  table  V-4. 
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DLTS  Signal:  AC/C 


Figure  V-27.  DLTS  spectrum  of  unirradiated  sample  Dl-38c  (Alo.30Gao.70N).  The 
Arrhenius  plot  of  peaks  PO  and  P 1  are  shown  in  the  inset. 
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Figure  V-28.  DLTS  spectra  of  unirradiated  Al0.30Ga0.70N  (sample  D 1-3  8c)  at  four 
different  reverse  bias  voltages.  Inset  is  the  nonnalized  peak  PO  for  each  of  the  traces, 
showing  that  the  peak  shape  is  unchanged  with  reverse  bias  voltage. 
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Figure  V-29.  Normalized  peak  height  of  electron  trap  PO  in  Al0.30Ga0.70N  versus  reverse 
bias  voltage. 
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Table  V-4.  Measured  defect  parameters  in  unirradiated  Alo.30Gao.70N  and  references  to 
similar  reports.  _ 


Parameter 

P0 

PI 

P2 

UeakVK) 

152 

369 

440 

Etot  (eV) 

0.94  ±0.06 

Coo  (1016  cm2) 

30000 

Similar  Traps  and 
Ejoj  (eV)  in  GaN 

D  0.25 

E02  0.27 

El  0.26 
DLN,93  0.24 

Similar  Traps  and 
^TOT(eV)  in 
AlxGai.xN 

PI98  0.16 
x=0.09 

EOA2 

0.57 
x  =  0. 12 

“Peak  temperature  determined  with  47s"1  rate  window. 


5.1.8  Irradiated  AI0.30Ga0.70N  (Wafer  D1) 

Figure  V-30  shows  the  DLTS  spectrum  of  irradiated  Alo.30Gao.70N  sample  Dl-38c 
measured  after  0.62  MeV  electron  irradiation  to  a  fluence  of  2.5x10  ctn'".  The  rate 
window  is  47  s'1,  and  the  pulsing  parameters  are  Vf=  3  V,  Vr  =  -5V,  and  tp  =  20  ps. 

Under  these  conditions,  no  radiation-induced  peaks  were  observed  in  the  120  to  290  K 
temperature  range.  It  should  be  noted  that  the  C-V  characteristics  of  samples  from  wafer 
D1  are  poor  due  to  high  series  resistance.  Thus,  the  DLTS  spectra  of  both  the 
unirradiated  and  irradiated  samples  were  suppressed  in  magnitude.  This  fact  makes  it 
difficult  to  determine  whether  or  not  the  2.5x10  cm'  fluence  of  0.62  MeV  electrons 
created  electron  traps  corresponding  to  those  seen  in  the  x  =0  to  0.20  AlxGai_xN  samples. 
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DLTS  Signal:  AC/C 
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Figure  V-30.  DLTS  spectra  of  Al0.30Ga0.70N  (sample  Dl-38c)  before  and  after  0.62  MeV 
electron  irradiation. 


5.1.9  DLTS  Analysis 


Figure  V-31  compares  the  46  s"1  rate  windows  of  unirradiated  AlxGai_xN  for 
x=0.0,  0.14,  0.20,  and  0.30.  Two  significant  electron  traps,  PI  and  P2,  appear  in 
unirradiated  AlGaN  at  all  aluminum  mole  fractions.  The  electron  trap,  P2,  appears  to  be 
a  superposition  of  traps,  as  yet,  unidentified.  However,  the  peak  trends  in  figure  V-3 1 
suggest  that  P2  could  be  a  superposition  of  traps  A  and  B.  These  two  traps  are  both 
thought  to  be  related  to  Voa-shallow  donor  complexes  [11],  and  are  present  in  GaN 
grown  by  all  kinds  of  methods.  Trap  PI,  which  has  been  previously  reported  in 
Alo.12Gao.88N  [42],  was  discussed  in  tenns  of  its  relationship  to  another  previously 
reported  trap,  DLN2.  This  DLN2  trap  was  reported  by  Gotz,  et  al.  [97]  as  being  donor- 
like.  In  figure  V-17,  it  is  shown  that  this  trap  does  not  exhibit  Poole-Frenkel  barrier 
lowering  at  low  electric  field  strength.  Additionally,  the  dislocation  related  capture 
kinetics  shown  in  the  inset  of  figure  V-18  suggests  that  trap  PI  is  acceptor-like.  It  may 
be  concluded  that  Gotz,  et  al.  were  either  observing  a  defect  that  was  different  from  PI  or 
that  they  mistakenly  attributed  Poole-Frenkel  barrier  lowering  to  what  was  actually 
phonon-assisted  tunneling  at  the  high  field  strength.  Figure  V-32  compares  DLTS 
spectra  of  the  radiation-induced  electron  traps  measured  in  AlxGai_xN  at  x=0.14  and  0.20. 
It  is  clear  that  the  radiation-induced  peaks  Rl,  R2,  R3,  and  R4  are  induced  at  each 
aluminum  mole  fraction,  and  that  the  peak  positions  of  these  traps  shift  upward  in 
temperature  with  aluminum  mole  fraction.  Arrhenius  analysis  of  radiation-induced  traps 
R2,  R3,  and  R4  indicates  that  this  peak  temperature  shift  is  caused  by  an  increase  in  F'tot. 
The  same  can  be  confirmed  for  as-grown  electron  traps  PI  and  P2.  Figure  V-3 3  plots  the 
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Figure  V-31.  DLTS  spectra  of  unirradiated  AlxGai_xN  for  x=0,  x=0.14,  x=0.20,  and 
x=0.30. 
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Figure  V-32.  DLTS  spectra  of  AlxGai_xN  (x=0.14  and  x=0.20)  before  and  after  1.0  MeV 
electron  irradiation. 
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measured  £tot  versus  x  for  several  different  electron  traps.  The  progression  of  £tot 
versus  x  for  electron  trap  P 1  is  found  to  be  equal  to  the  increase  in  band  gap  energy  as 
measured  by  CL  (table  III- 1).  This  trend  for  PI  is  shown  clearly  in  figure  V-34  where 
L'tot  is  plotted  versus  band  gap  energy.  In  this  figure,  a  straight  line  with  a  slope  of  1.0 
connects  the  three  P 1  data  points,  indicating  that  the  increase  in  £tot  is  equal  to  the 
increase  in  band  gap  energy  in  this  range.  In  figure  V-34,  the  dotted  line  represents  a 
possible  extension  of  AlGaN  trap  PI  to  GaN  trap  D.  The  L'tot  of  radiation-induced  traps 
R2,  R3,  and  R4  are  found  to  increase  significantly  with  increase  in  x,  but  by  a  lesser 
amount  than  the  increase  in  band  gap  energy.  The  slope  of  each  connecting  line  is 
labeled  in  figure  V-34.  It  is  apparent  in  this  figure  that  E-tot  of  the  shallowest  traps 
changes  the  least  with  increase  in  band  gap  energy  and  that  ETOt  of  the  deepest  traps 
changes  the  most  with  the  increase  in  band  gap  energy.  The  radiation-induced  defects  in 
GaN  corresponding  to  R1  and  R2  are  believed  to  be  VN  related.  Thus  R1  and  R2  in 
AlGaN  may  also  be  related  to  Vn-  While  these  defects  are  ubiquitous  in  irradiated  GaN, 
the  radiation-induced  trap  corresponding  to  R3  is  not  reported  by  all  groups.  However,  it 
has  been  observed  in  HVPE  GaN  by  Auret  et  al.  under  a  variety  of  different  radiation 
conditions  (Sr90  (3-decay  electron  source,  5.4  MeV  He-ions,  and  2  MeV  protons 
[96:262]).  They  reported  that  this  trap  did  not  exhibit  typical  Poole-Frenkel  electric  field 
enhanced  emission.  Given  this  fact,  and  that  the  capture  cross  section  was  relatively 
small  (a  =  8±4xl0‘  cm"),  the  authors  [96]  suggested  that  ER3  is  acceptor-like.  The 
radiation-induced  electron  trap  R4  observed  here  does  not  appear  to  correspond  to  the 
radiation-induced  trap  level  reported  previously  in  GaN.  This  raises  the  possibility  that 
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R4  is  related  to  aluminum  displacement  This  trap  is  introduced  with  1 .0  MeV  radiation  at 
a  greater  rate  than  the 


Aluminum  Mole  Fraction  (x) 


Figure  V-33.  Measured  electron  trap  energy  versus  aluminum  mole  fraction  for  traps 
labeled  R2,  R3,  R4,  D,  P0,  and  PI  in  this  study.  The  dotted  line  represents  a  possible 
extension  to  GaN  trap  D  with  decrease  in  aluminum  mole  fraction. 
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Figure  V-34.  Measured  electron  trap  energy  versus  band  gap  energy  for  traps  labeled 
R2,  R3,  R4,  D,  PO,  and  PI  in  this  study.  The  dotted  line  represents  a  possible  extensions 
to  GaN  trap  D  with  decrease  in  band  gap  energy. 


individual  introduction  rates  of  Rl,  R2,  and  R3.  Furthermore,  R4  appears  to  have  greater 
thermal  stability  than  the  other  three  traps.  The  existence  of  this  trap,  could,  therefore,  be 
a  very  important  contributor  to  radiation-induced  degradation  of  AlGaN  based  devices. 
Traps  A,  B,  C,  D,  PO,  PI,  and  P2  all  exhibit  dislocation-related  capture  kinetics  as 
discussed  in  chapter  IV.  According  to  the  theory  that  was  discussed  in  that  chapter,  these 
traps  are  expected  to  be  acceptor-like  defects  that  are  concentrated  in  the  vicinity  of  the 
multiple  dislocations  expected  in  thin,  epitaxial  layers  of  GaN  and  AlGaN.  All  of  the 
radiation-induced  traps  anneal  out  almost  entirely  following  annealing  at  400  °C  for  1 5 
minutes,  and  traps  Rl  and  R2  begin  annealing  (under  a  reverse  bias  electric  field)  below 
200  °C  in  Alo.14Gao.86N.  Comparing  figures  V-22,  V-23,  and  V-25,  the  radiation-induced 
traps,  R2,  R3,  and  R4  appear  to  have  greater  thennal  stability  in  Al0.2oGa0.8oN  than  those 
of  the  same  traps  in  Alo.14Gao.86N.  The  electron  trap  labeled  PO  in  Al0.30Ga0.70N  is  of 
unknown  origin,  and  it  has  an  unusually  large  capture  cross  section.  This  large  cross 
section  would  normally  indicate  a  donor-type  trap,  but  the  reverse  bias  voltage 
dependence  and  electron  capture  kinetics  do  not  appear  to  support  this  conclusion. 

Figure  V-35  shows  the  Arrhenius  plots  of  all  electron  traps  thus  characterized  in  this 
study.  The  corresponding  activation  energies  and  capture  cross  sections  can  be  found  in 
tables  V-l,  V-2,  V-3,  and  V-4 
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5.2  Hall  Effect  Measurements 


Hall  effect  measurements  of  AlxGai_xN  (x=0.0  to  0.30)  samples  were  made 
following  1.0  MeV  electron  irradiation  at  a  fluence  of  lxlO17  cm'2.  The  room 
temperature  carrier  removal  rate,  r),  was  calculated  as  the  change  in  carrier  concentration 
divided  by  the  fluence.  It  is  found  that  to  first  order,  the  carrier  removal  rate,  p,  is 
dominated  by  the  initial  carrier  concentration,  n0.  Figure  V-36  shows  that  the  carrier 
removal  rate,  r|,  has  an  approximately  linear  relationship  to  n0,  regardless  of  the 
aluminum  mole  fraction.  This  linear  fit  is  quantified  as  p  =  (3.96x10"  «0- 0.15)  cm"  for 
6.5x1016<  n0<  8.2xl017  cm"3.  Under  the  same  irradiation  and  characterization  conditions, 

17  3 

the  GaAs  sample  Fl-8  (n0  =  6.41x10  cm"  )  is  found  to  have  a  carrier  removal  rate  p  = 

5.1  cm"  ,  while  the  corresponding  p  in  AlxGai_xN  having  the  same  n„  is  2.4  cm"1.  Thus, 
the  removal  rate  due  to  1 .0  MeV  electron  irradiation  in  this  GaAs  sample  is  two  times 
greater  than  in  AlxGai_xN. 

Temperature  dependent  carrier  concentration  was  analyzed  in  order  to  detennine 
the  mechanisms  involved  in  this  carrier  removal  processes.  Figure  V-37  shows  the 
temperature  dependent  carrier  concentration  of  GaN  sample  A2-6  before  and  after  the 
l.lxlO17  cm"2  fluence  of  1.0  MeV  electrons.  The  fit  lines  are  produced  by  using  the  two 
donor  model  of  equation  IV- 16.  Before  irradiation,  the  fitting  parameters  indicate  that 
there  are  1.9x10  cm"  shallow  donors,  Ndi,  having  a  screened  thennal  activation  energy, 

16  3 

ED i,  equal  to  17  me V.  Also,  there  are  7.5x10  cm"  deeper  donors,  ND2,  having  a 
screened  thermal  activation  energy,  Ed2,  equal  to  59  meV.  There  are  also  2.3x10  cm" 
deep  (always  ionized)  acceptors,  NA.  Following  irradiation,  NDi  drops  to  1.64xl017  cm"3, 
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Figure  V-36.  Dose-averaged  carrier  removal  versus  initial  room  temperature  carrier 
concentration  in  AlxGai_xN. 
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Figure  V-37.  Temperature  dependent  Hall  data  for  GaN  (sample  A2-1)  before  and  after 
1.0  MeV  electron  irradiation.  Model  (lines)  are  for  a  two  donor  fit  with  (Ndi=1.9x1017, 
EDi=17  meV,  ND2=7.5xl016,  ED2=60  meV,  NA=2.3xl016)  and  (NDi=1.64x1017,  EDi=19 
meV,  Nd2=1.55x1017,  Ed2=59  meV,  NA=8.7xl016)  for  the  unirradiated  and  irradiated  data 
respectively. 
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Nd2  increases  to  1.55x10  cm"  and  Na  increases  to  8.7  xlO  cm'  .  The  unscreened 
thermal  activation  energy  of  the  shallow  donor  is  calculated  by  equation  III-3  to  be  Edi° 

=  17  +  (2.1x10'5)(1.9x1017)1/3  =  29  meV.  This  is  a  very  reasonable  value  for  Si  in  GaN. 
The  deeper  donor  has  an  unscreened  thermal  activation  energy  ED2°  =  59  +  (2.  lxl O' 
5)(7.5x1016)1/3  =  68  meV.  This  is  a  very  reasonable  value  for  VN  in  GaN,  agreeing  with 
that  reported  by  Look,  et  al.  [47].  With  confidence  in  the  assignment  of  the  two  donor 
levels,  physical  inferences  may  be  drawn  from  the  fitting  parameters.  Following 
irradiation,  the  total  increase  in  the  concentration  of  Vn  is  8.0xl016  cm'3.  This  is  nearly 
equal  to  the  sum  of  the  decrease  in  Si  donors,  2.6xl016  cm'3  and  the  increase  in  NA, 
6.4x10  cm"  .  This  situation  can  be  understood  as  5.4x10  cm'  nitrogen  interstitials 
acting  as  deep  single  acceptors  according  to  theory  [24,  29],  2.6xl016  cm"3  nitrogen 
interstitials  forming  Si-Ni  complexes  that  passivate  some  silicon  donors  [99],  and 
3.3x10  cm"  gallium  interstitials  acting  as  triple  acceptors  according  to  theory  [24,  29]. 
Assuming  no  recombination  of  nitrogen  Frenkel  pairs,  the  introduction  rate  of  Vn  with 
1.0  MeV  electrons  is  0.73  cm'1  according  to  the  above  analysis.  Referencing  figure  II- 10, 
it  is  seen  that  this  introduction  rate  would  correspond  to  a  nitrogen  displacement  energy 
of  approximately  25  eV,  which  falls  within  a  theoretically  reasonable  range  of  20  to  50 
eV  for  GaN.  Based  on  the  above  considerations,  however,  there  is  no  reason  that  the 
carrier  removal  rate  would  be  expected  to  increase  as  the  initial  carrier  concentration  is 
further  increased.  Normally,  the  increase  in  Si  donors  would  raise  the  odds  that  a  mobile 
Ni  forms  a  Si-Ni  complex  before  recombining  with  a  Vn.  In  order  for  q  to  increase  over 
the  whole  range  shown  in  figure  V-36,  based  on  this  mechanism  alone,  the  nitrogen 
displacement  introduction  rate  would  have  to  be  over  3  cm'1.  This  displacement 
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introduction  rate  is  3  times  larger  than  predicted  by  equations  11-30  through  11-32  and  the 
lowest  reasonable  value  of  displacement  energy,  Tj.  Thus,  other  mechanisms  of  carrier 
removal  such  as  ionization-induced  activation  of  acceptor  impurities  [7,  96]  may  be 
dominant  for  large  n0  in  figure  V-36. 

It  was  found  that  the  substrate  contribution  to  room  temperature  conductivity  was 
greater  than  or  equal  to  the  bulk  contribution  in  samples  from  AlGaN  wafers  Bl,  B2,  Cl, 
C2,  and  Dl.  In  these  cases,  the  In (n)  versus  1000/T  plot  was  very  sensitive  to  the  choice 
of  two-layer  correction  parameters  (low  temperature  values  of  ns  and  /iv).  An  example  of 
this  behavior  is  shown  for  unirradiated  Alo.14Gao.86N  sample  B2-16  in  figure  V-38  and  for 
unirradiated  Al0.30Ga0.70N  sample  D2-7  in  figure  V-39.  A  reliable  two-donor  fit  could  not 
be  obtained  for  any  of  the  AlGaN  samples,  whether  they  were  irradiated  or  unirradiated. 

It  may  be  that  for  these  samples,  the  two  layer  analysis  of  equations  III- 13  and  III- 14  is 
inadequate  to  yield  the  precise  bulk  carrier  concentration  and  mobility.  Despite  the 
apparent  inadequacy  of  the  two  layer  analysis  for  precise  temperature  dependent  Hall 
fitting  of  the  AlGaN,  it  is  very  necessary  and  helpful  for  determining  the  true  carrier 
concentration  and  mobility  at  room  temperature  and  is  applied  accordingly  in  this 
dissertation  research.  Figures  V-40,  V-41,  and  V-42  are  representative  of  the  temperature 
dependent  Hall  measurements,  corrected  for  substrate  conduction,  in  AlGaN  samples. 

The  Alo.27Gao.73N/GaN  single  heterostructure  sample  El -3  was  subjected  to  the 
9x10  cm'  electron  fluence  at  1.0  MeV  and  characterized  with  temperature  dependent 
Hall  measurements.  The  Ohmic  contacts  were  diffused  through  the  Alo.27Gao.73N  layer 
and  brought  into  contact  with  the  two-dimensional  electron  gas  (2-DEG)  by  annealing  at 
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Figure  V-38.  Carrier  concentration  versus  inverse  temperature  of  unirradiated 
Alo.14Gao.86N  (sample  B2-16)  before  and  after  correction  for  a  secondary  conduction 
layer.  The  corresponding  mobility  plots  are  shown  in  the  inset. 
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Figure  V-39.  Carrier  concentration  versus  inverse  temperature  of  unirradiated 
Al0.30Ga0.70N  (sample  D2-7)  before  and  after  correction  for  a  secondary  conduction  layer. 
The  corresponding  mobility  plots  are  shown  in  the  inset. 
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Figure  V-42.  Temperature  dependent  carrier  concentration  of  Alo.30Gao.70N  (sample  D2 
7)  before  and  after  1 .0  MeV  electron  irradiation. 


900  °C.  The  high  conductivity  of  the  2-DEG  region  causes  the  AlGaN  layer  to  be  shorted 
and  it  can  be  assumed  that  the  measured  sheet  carrier  concentration  and  mobility  reflect 
that  of  the  2-DEG  region.  Figures  V-43,  V-44,  and  V-45  show  the  temperature 
dependent  sheet  carrier  concentration,  mobility,  and  sheet  resistivity  of  the  sample  before 
and  after  irradiation.  It  is  found  that  sheet  carrier  concentration  and  mobility  are  reduced 
by  13  %  and  23%,  respectively,  at  room  temperature  (295  K).  This  causes  a  150  % 
increase  in  sheet  resistivity.  While  sheet  carrier  concentration  remains  relatively  constant 
with  temperature  before  and  after  irradiation,  the  mobility  is  greatly  affected  at  low 
temperatures,  where  ionized  defect  scattering  dominates.  Thus,  it  is  confirmed  that  the 
1.0  MeV  electron  doses  used  throughout  this  study  (~  lxlO17  cm"2)  and  the  accompanying 
defect  introduction  rates  are  relevant  to  HEMT  device  operation. 
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Figure  V-43.  Sheet  carrier  concentration  versus  temperature  in  Alo.27GaNo.73N/GaN 
single  heterostructure  (sample  El -3)  before  and  after  1.0  MeV  electron  irradiation. 
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Figure  V-44.  Hall  Mobility  versus  temperature  in  Alo.27GaNo.73N/GaN  single 
heterostructure  (sample  El -3)  before  and  after  1.0  MeV  electron  irradiation. 
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Figure  V-45.  Sheet  resistivity  versus  temperature  in  Alo.27GaNo.73N/GaN  single 
heterostructure  (sample  El -3)  before  and  after  1.0  MeV  electron  irradiation. 
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5.3  Cathodoluminescence 


The  CL  spectra  of  as-grown  AlxGai_xN  samples  are  shown  in  figure  V-46.  These 
and  all  CL  measurements  presented  here  were  made  at  T  =  6  K  with  a  25  pA  excitation 
current  of  5  keV  electrons.  In  the  GaN  spectrum,  a  donor  bound  exciton  (D°,X)  peak  is 
dominant  at  3.48  eV.  Lesser  features  are  typical  of  that  previously  reported  in  GaN.  The 
CL  spectrum  of  the  Alo.14Gao.86N  shows  a  (D°,X)  at  3.87  eV,  and  a  secondary  peak  at 
3.83  eV  that  may  be  a  neutral-donor-to-free  hole  transition  (D°,h)  [100].  The  CL 
spectrum  of  the  Alo.20Gao.80N  shows  the  (D°,X)  peak  at  4.00  eV,  and  the  secondary  peak 
at  3.92  eV.  The  CL  spectrum  of  the  Al0.30Ga0.70N  shows  the  (D°,X)  peak  at  4.26  eV,  and 
the  secondary  peak  at  about  4.19  eV.  Two  broad  luminescence  peaks  appear  at  lower 
energies  in  the  AlxGai_xN.  One  of  these  (as  designated  in  figure  V-46)  corresponds  to  the 
ubiquitous  yellow  luminescence  (YL)  band  in  GaN.  Interestingly,  another  broad 
luminescence,  labeled  UVL,  is  fixed  at  1.00  ±  0.05  eV  below  the  band  edge.  It  is 
probable  that  this  luminescence  is  caused  by  a  radiative  transition  from  a  defect  band 
centered  at  about  1.00  ±  0.05  eV  below  the  conduction  band  minimum  to  holes  in  the 
valence  band.  This  defect  band  may  correspond  to  the  observed  peak  P2  discussed  in  the 
DLTS  results  section  above.  There,  it  was  shown  that  the  apparent  activation  energy  of 
P2  was  0.98  ±  0.10  eV  in  Alo.14Gao.86N  and  0.92  ±  0.06  eV  in  Alo.20Gao.80N.  The  CL 
spectra  of  as-grown  and  irradiated  GaN  and  Alo.20Gao.80N  samples  are  shown  on  a  linear 
scale  in  figure  V-47.  The  CL  spectra  measured  from  the  as-grown  samples  are  depicted  in 
solid  lines  and  the  CL  spectra  from  the  irradiated  samples  are  in  dotted  lines.  Following 
1.0  MeV  electron  irradiation  at  a  fluence  of  l.lxlO17  cm"2,  the  peak  CL  luminescence 
intensities  of  both  GaN  and  Al0.2oGa0.8oN  decreased  by  about  50%  as  shown  in  the  figure. 
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Figure  V-46.  Normalized  CL  spectra  of  unirradiated  AlxGai_xN  for  x  =  0,  0.14,  0.20,  and 
0.30. 
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Figure  V-47.  Normalized  CL  spectra  of  AlxGai_xN  (x=0.0  and  x=0.2)  before  and  after 

17  2 

1 .0-MeV  electron  irradiation  with  a  fluence  of  1 . 1x10  cm'” 
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This  reduction  is  found  to  be  stable  at  room  temperature  and  is  attributed  to  the  creation 
of  non-radiative  recombination  centers.  Figure  V-48  shows  a  plot  of  the  normalized  CL 
peak  luminescence  intensity  of  five  different  AlxGai_xN  samples  as  a  function  of  the  1.0 
MeV  electron  fluence.  The  intensity  is  reduced  by  about  35  and  50%  after  electron 
irradiation  at  fluences  of  2xl016  and  l.lxlO17  cm"2,  respectively  from  the  value  before 
irradiation.  The  margin  of  error  can  be  rather  large  for  these  measurements,  and  it  is  not 
expected  that  any  trend  can  be  distinguished  based  on  the  compositional  differences  of 
the  samples.  Rather,  the  data  of  figure  V-48  are  averaged  and  taken  to  be  representative 
of  the  AlxGai_xN  samples  for  0  <  x  <  0.30,  and  this  average  is  shown  in  figure  V-49.  The 
lines  overlaid  in  figure  V-49  are  the  reported  electroluminescence  peak  intensities  from 
InGaN  LEDs  irradiated  with  0.94  MeV  (dotted  line)  and  1.14  MeV  (dashed  line) 
electrons  [53].  Interestingly,  the  exponential  character  of  all  the  curves  is  similar.  In  the 
referenced  work,  it  was  reported  that  the  damage  annealed  out  after  16  hours  at  room 
temperature.  Figures  V-50  to  V-55  show  the  spectra  of  several  AlxGai_xN  samples  (x  = 
0.0  to  0.30)  before  and  after  1.0  MeV  electron  irradiation.  These  are  presented  on  a 
logarithmic  scale  to  show  the  detail  of  the  broad  range  of  luminescence  peaks.  Generally, 
the  entire  spectrum  decreases  in  intensity  due  to  the  introduction  of  non-radiative 
recombination  centers.  Figure  V-56  shows  the  GaAs  sample  Fl-8  before  and  after  1.0 
MeV  electron  irradiation  to  a  fluence  of  9x10  cm'  .  The  peak  luminescence  intensity 
decreases  by  99.9  %  following  the  irradiation.  Thus,  luminescence  intensity  in  the 
AlxGai_xN  samples  is  much  more  tolerant  to  1.0  MeV  electron  radiation  than  in  the  GaAs 
sample. 
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Figure  V-48.  Normalized  luminescence  intensity  in  AlxGai_xN  irradiated  at  two 
consecutive  1.0  MeV  electron  doses  with  2xl016  and  l.lxlO17  e'cm"2. 
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Figure  V-49.  Normalized  luminescence  intensity  in  AlxGai.xN  irradiated  at  two 
consecutive  1.0  MeV  electron  doses  with  2xl016  and  l.lxlO17  e'cm"2.  Data  are  produced 
from  the  average  of  four  samples  (x=0.0  to  0.20).  The  reported  degradation  of  an 
electron  irradiated  InGaN  LED  [53]  is  overlaid  for  comparison. 
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Figure  V-50.  CL  spectrum  of  GaN  (sample  Al-16)  before  and  after  1.0  MeV  electron 
irradiation  with  a  fluence  of  l.lxl017cm"2. 
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Figure  V-51.  CL  spectrum  of  GaN  (sample  A2-1)  before  and  after  1.0  MeV  electron 
irradiation  with  a  fluence  of  1 . 1  x  1 017  cm  2. 
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Figure  V-52.  CL  spectrum  of  Alo.14Gao.86N  (sample  B2-39)  before  and  after  1.0  MeV 
electron  irradiation  with  a  fluence  of  l.lxlO17  cm  2. 
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Figure  V-53.  CL  spectrum  of  Alo.20Gao.80N  (sample  Cl -9)  before  and  after  1.0  MeV 
electron  irradiation  with  a  fluence  of  l.lxlO17  cm  2. 
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Figure  V-54.  CL  spectrum  of  Alo.20Gao.80N  (sample  C2-2)  before  and  after  1.0  MeV 
electron  irradiation  with  a  fluence  of  1.1x1 017  cm"'. 
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Figure  V-55.  CL  spectrum  of  Alo.2oGa0.8oN  (sample  D2-7)  before  and  after  1.0  MeV 

16  2 

electron  irradiation  with  a  fluence  of  9x10  cm"  . 
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Figure  V-56.  CL  spectrum  of  GaAs  (sample  Fl-8)  before  and  after  1.0  MeV  electron 
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irradiation  with  a  fluence  of  9x  1 0  cm'". 
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VI.  Conclusions  and  Recommendations 


6.1  Conclusions 

The  primary  objective  of  this  study  was  to  quantify  the  effects  of  electron 
radiation  upon  n-type  AlxGai_xN  (x  =  0.0  to  0.3)  grown  by  MBE  with  regard  to  carrier 
removal,  mobility  reduction,  and  luminescence  degradation,  and  to  electrically 
characterize  radiation-induced  defects  that  may  contribute  to  performance  degradation  of 
nitride-based  electronic  and  optoelectronic  devices.  A  secondary  objective  was  to 
characterize  the  electrically  active  defects  that  exist  in  the  unirradiated  AlxGai.xN 
material,  and  to  compare  these  defects  with  those  in  irradiated  material. 

The  research  results  from  this  dissertation  research  are  as  follows: 

1)  It  has  been  found  that  with  1.0  MeV  electrons  at  a  fluence  of  lxlO17  cm’2,  the  carrier 
removal  rate  depends  foremost  on  the  initial  carrier  concentration,  regardless  of 
aluminum  mole  fraction.  The  carrier  removal  rate,  q,  was  detennined  to  be  a  linear 
function  of  the  initial  carrier  concentration,  n0,  as  q  =  (3.96x10’  na  -  0.15)  cm’  .  This 
finding  provides  a  framework  for  understanding  seemingly  conflicting  reports  about 
carrier  removal  rates  in  GaN  and  AlGaN. 

2)  The  result  in  1)  supports  the  hypothesis  that  for  1.0  MeV  electron  irradiation,  a 
mechanism  of  carrier  removal  in  Si-doped  AlxGai_xN  is  donor  passivation  via  the 
formation  of  Si-N;  complexes.  This  Si-Ni  complex  hypothesis  is  also  supported  by  fitting 
temperature  dependent  carrier  concentration  to  a  general  two-donor  model.  It  was  found 
that  following  irradiation  of  GaN,  the  sum  of  passivated  shallow  donors  and  deep 
acceptors  is  approximately  equal  to  an  increase  in  the  number  of  presumed  nitrogen 
vacancies. 
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3)  A  comparison  of  carrier  removal  and  luminescence  degradation  in  AlxGai_xN  and 
GaAs  was  made  for  irradiated  samples,  and  they  were  characterized  under  identical 
conditions.  This  comparison  study  supports  the  claims  that  AlxGai_xN  is  intrinsically 
more  tolerant  to  electron  radiation  than  conventional  semiconductor  materials  such  as 
GaAs,  and  especially  so  with  regard  to  luminescence  degradation. 

4)  It  was  found  that  three  radiation-induced  deep  defect  levels  R1  (0.15  eV),  R2  (0.21 
eV),  and  R3  (0.26  eV)  exist  in  Alo.14Gao.86N  which  correspond  to  previously  reported 
radiation-induced  deep  defect  levels  ER1  (~0.1 1  eV),  ER2  (-0.16  eV),  and  ER3  (-0.20 
eV)  in  GaN.  The  corresponding  AlGaN  levels  were  shown  to  deepen  significantly  with 
increasing  aluminum  mole  fraction. 

5)  A  new  radiation-induced  defect  level  R4  (0.33  eV  for  x  =  0. 14,  and  0.38  eV  for  x  = 
0.20)  has  been  found  in  AlGaN,  and  it  does  not  appear  to  correspond  to  any  known 
radiation-induced  defect  level  in  GaN.  The  defect  was  observed  to  be  thermally  stable 
under  reverse  bias  voltage  at  450  K,  and  it  did  not  anneal  out  completely  until  after  600 
K.  Knowledge  of  this  and  the  other  three  radiation-induced  defect  levels  enables  more 
complete  theoretical  modeling  of  radiation-induced  degradation  in  the  gate  region  of 
AlGaN/GaN  HEMTs,  and  could  help  device  manufacturers  engineer  radiation  hardened 
designs. 

6)  Two  dominant  defect  levels,  PI  and  P2,  in  unirradiated  AlGaN  were  further 
characterized.  It  was  demonstrated  that  PI  is  unlikely  to  be  of  the  same  origin  as  a 
donor-type  electron  trap  reported  in  an  early  AlGaN  DLTS  study.  It  was  found  that  for 
PI,  the  increase  in  E'tot  with  aluminum  mole  fraction  was  equal  to  the  increase  in  band 
gap  energy.  Additionally,  CL  characterization  suggests  that  there  exists  a  radiative 
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transition  from  the  P2  level  to  the  valence  band.  The  energy  of  this  transition  also 
increases  with  aluminum  mole  fraction  according  to  the  increase  in  band  gap. 

6.2  Recommendations 

The  observation  of  dominant  radiation-induced  electron  traps  in  AlxGai_xN  and 
the  discovery  of  the  conditions  for  which  they  are  well-resolved  provides  impetus  for 
future  DLTS  characterization  of  these  traps  under  different  experimental  conditions.  The 
logical  next  step  in  electrical  characterization  of  these  traps  is  to  detennine  the  activation 
energy  capture  barrier,  Ea,  for  each  trap,  so  that  the  true  thermal  activation  energy,  Eti„ 
may  be  extracted  from  the  total  activation  energy,  Ejoi,  presented  here.  While  there  are 
several  ways  to  attempt  such  a  measurement,  the  most  straightforward  and  reliable 
investigation  would  be  to  measure  the  fdled  trap  concentration  X(C)  versus  tp,  for  10'9<  tp 
<  1CT6  s.  From  this,  the  true  capture  cross  section  on{  T)  and  thus  Ea  may  be  extracted. 
Another  way  in  which  a  CV-DLTS  system  could  be  extended  to  better  characterize  these 
traps  is  to  incorporate  optical  excitation.  When  performed  successfully,  a  direct 
measurement  of  the  Gibbs  free  energy  (Ej)  is  obtained,  but  more  importantly,  the  DX 
behavior  and  associated  persistent  photoconductivity  (PPC)  could  be  explored. 

Additional  traps  (both  electron  and  hole  traps)  could  be  probed  throughout  the  entire 
band  gap.  While  many  additional  techniques  could  be  applied  to  the  characterization  (and 
even  positive  identification)  of  these  new  radiation-induced  traps,  this  research  can  also 
be  extended  by  determining  how  the  observed  traps  affect  the  operation  of  AlGaN  based 
electronic  and  optoelectronic  devices.  Such  an  investigation  should  include  DC  and  RF 
characterization. 
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Appendix:  C-V  and  l-V  Measurements 


The  greatest  technical  challenge  encountered  in  this  dissertation  research  was  to 
obtain  good  quality  Schottky  and  Ohmic  contacts  that  permitted  reliable  electrical 
characterization.  The  Schottky  barrier  diodes  used  for  DLTS  characterization  were 
particularly  challenging  to  obtain  because  they  required  both  types  of  contacts  and  also 
were  subject  to  other  requirements  such  as  the  requirement  that  the  depletion  region  not 
extend  to  the  substrate  and  that  series  resistance  be  sufficiently  small  at  T  >  100  K. 
Adequate  SBDs  could  not  be  successfully  fabricated  from  the  AlGaN  wafers  Bl,  C2,  and 
D2.  The  DLTS  results  obtained  with  the  good  SBDs  were  presented  in  chapter  V,  and 
representative  C-V  and  I-V  plots  from  these  devices  are  included  in  this  appendix  for 
reference.  The  1/C2  vs.  V  plot  insets  show  the  calculated  carrier  concentration  versus 
depth,  and  these  insets  also  show  how  the  depletion  depth  increased  with  applied  reverse 
bias.  Aside  from  carrier  concentration  determination,  C-V  and  I-V  curves  were  only 
utilized  in  this  dissertation  research  as  a  qualitative  diagnostic  of  problems  encountered 
during  DLTS  measurements.  The  I-V  vs.  T  plots  give  an  indication  of  the  reverse  bias 
leakage  current  and  series  resistance  encountered  at  various  temperatures.  Following 
irradiation,  the  most  noticeable  changes  in  the  SBDs  resulted  from  the  increase  in  series 
resistance.  As  such,  leakage  current  usually  decreased  following  irradiation  and  the 
DLTS  peak  heights  were  reduced  slightly  compared  with  those  seen  before  irradiation. 
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Figure  A-l.  Plot  of  1/C  versus  applied  bias  for  unirradiated  GaN  (sample  A2-4b)  at 
T=290  K.  Inset  is  the  extrapolated  carrier  concentration  versus  depletion  depth. 
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Figure  A-3.  C-V  curves  at  various  temperatures  for  unirradiated  GaN  (sample  A2-4b). 
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Figure  A-4.  I-V  curves  at  various  temperatures  for  unirradiated  GaN  (sample  A2-4b). 
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Figure  A-5.  Plot  of  1/C  versus  applied  bias  for  irradiated  GaN  (sample  A2-5a)  at  T  = 
150  K.  Inset  is  the  extrapolated  carrier  concentration  versus  depletion  depth. 
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Figure  A-6.  Plot  of  1/C  versus  applied  bias  for  irradiated  GaN  (sample  A2-5a)  at  T  = 
1 10  K.  Inset  is  the  extrapolated  carrier  concentration  versus  depletion  depth. 
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Figure  A-7.  I-V  curves  at  various  temperatures  for  1.0  MeV  electron  irradiated  GaN 
(sample  A2-4a). 
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Figure  A-8.  Plot  of  1/C  versus  applied  bias  for  unirradiated  Alo.14Gao.86N  (sample  B2- 
2 Id)  at  T=294  K.  Inset  is  the  extrapolated  carrier  concentration  versus  depletion  depth. 
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Figure  A-9.  Plot  of  1/C"  versus  applied  bias  for  unirradiated  Alo.14Gao.86N  (sample  B2- 
2 Id)  at  T  =  240  K.  Inset  is  the  extrapolated  carrier  concentration  versus  depletion  depth. 
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Figure  A-10.  Capacitance  versus  temperature  in  Alo.14Gao.86N  (sample  B2-21b) 
following  1.0  MeV  electron  irradiation. 
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Figure  A-ll.  I-V  curves  of  irradiated  Alo.2oGa0.8oN  (sample  Cl-38a)  before  and  after 
annealing  at  350  °C  for  15  min. 
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